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Abstract 



The Sunyaev-Zel'dovich Array: 

Constraining a New Pressure Profile for Fitting SZE Observations of 

Galaxy Clusters 

Tony Mroczkowski 

The Sunyaev-Zel'dovich Array (SZA), an eight element interferometer designed 
to probe the Sunyaev-Zel'dovich effect (SZE) from galaxy clusters, which I helped 
construct and operate, is described here (Part I). I then use SZA observations to 
investigate the utility of a new, self-similar pressure profile for fitting SZE observations 
of galaxy clusters (Part II). 

The SZA 30-GHz receiver system probes angular scales ~ 1-5'. A model that can 
accurately describe a cluster's pressure profile over a correspondingly broad range 
of radii is therefore required. In the analysis presented here, I fit a 2-parameter, 
radial pressure profile, derived from simulations and detailed X-ray analysis of relaxed 
clusters, to SZA observations of three clusters with exceptionally high quality X-ray 
data. Prom the joint analysis of the SZE and X-ray data, I derive physical properties 
of the cluster, such as gas and total mass, gas fraction and the integrated Compton 
y-parameter. 

The parameters derived from the joint fit to SZE-|-X-ray data agree well with a 
detailed, independent. X-ray-only analysis of these same clusters. When combined 
with X-ray imaging data, this new pressure profile yields an independent estimate of 
the electron temperature profile that is in good agreement with spectroscopic X-ray 
determinations. In addition to yielding relationships between cluster observables and 
physical cluster properties, this model could prove to be a useful tool in helping to 



constrain the temperatures of high redshift clusters, for which X-ray spectroscopic 
data are difficult to obtain. 
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1.1 

1.2 



1.3 



1.4 



Image adapted from van Speybroeck ( 1999 ) 



SZE spectral dependence /(z^), plotted as a function of frequency (see 
Eq. 1.5). The upper panel shows /(z/) over a broad range of frequencies, 



while the lower panel shows a detail of the null in the SZE spectrum. The 
relativistically-corrected /(z/) for a cluster with temperature kT^. = keV 
(black, solid line) reduces precisely to the classical frequency dependence 



(Eq. 1.4), since electrons with no temperature are not moving at rela- 
tivistic random velocities. The other lines show how the relativistically- 
corrected /(z/) departs from the classical behavior for higher temperature 
electrons. The relativistic corrections to /(z/) shown here are computed 



out to fifth-order using the equations provided in Itoh et al. (1998). The 
classical SZE spectrum has a null at z/ 217.5 GHz, above which the 
SZE signal becomes an increment. Higher temperature electrons require 



relativistic corrections (Itoh et al. 1998) to the classical SZE frequency 



dependence, which shift the null to higher frequencies. A high tempera- 
ture cluster would have a non-negligible thermal SZ effect at the classical 
null (217.5 GHz) 



SZE frequency dependence below 100 GHz. See Fig. |1.2| for caption. Be- 
low the null (see lower panel of Fig. 1.2), we often use the term "SZE 



decrement" to refer to the strength of the thermal SZ effect, which is neg- 
ative at low frequencies. The lower panel shows the fractional deviation 
from the classical SZE for clusters at higher temperatures, due to the 
relativistic electron velocities in the ICM. Treating the thermal SZE at 
30 GHz from a massive 10 keV cluster as classical introduces a ~ —3.6% 
bias in quantities derived from the SZE fits (i.e. the line-of-sight electron 
pressure in Eq. |1.2| would be underestimated by this amount, since the 



strength of the SZE would be overestimated by the classical calculation). 
Example baseline formed by a pair of antennae. The dish separation is 
d, while the projected baseline as seen by the source is d'. The distance 
to the source (represented by a cloud) is labeled Rjf. Note the scale 
of the broken lines, representing this astronomical distance -R//, is not 
accurate; in reality, the lines of sight are nearly parallel and point at the 
same location on the source. The time delay T is equal to the difference 
in distances to the source, from each antenna, divided by the speed of 
light c. By delaying the signal measured by the antenna on the right by 
time T, we ensure that the same astronomical wavefront is used in the 
correlation of the two signals. Tracking an astronomical source requires 
both the physical pointing of each antenna toward the source in the far 
field, and the implementation of instrumental delays (T) that ensure each 
antenna measures the same wave front as the Earth rotates 
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1.5 



1.7 



1.9 



SZA n,^;-space coverage of CL1226.9+3332 (12:26:58.0, +33:32:45.0), a 
cluster that passes near zenith for the SZA (the SZA's latitude is ~ 37°N). 
Each blue point represents a data point's n,f-space location, using fre- 
quency independent units of kX (distance divided by wavelength, divided 
by 1000). Note that there are 2 groupings of points: those ~ 0.35-1.3 kA 
from the center, due to the short baselines of the compact inner array, 
and those ~ 3-7.5 kA from the center, due to the long baselines formed 
with the outer antennae (i.e. each baseline formed with Antennae 6 or 7; 
the antenna layout of the SZA is shown in Fig. 1.7). The SZA was de- 
signed to provide this broad, uniform coverage in -space on the shorter 
baselines, while simultaneously providing higher resolution coverage with 
longer baselines. The data locations in -space exhibit inversion sym- 
metry (i.e. {u,v) = —{u,v)) because the visibilities are the transform 
of real data in image space. The visibilities therefore exhibit Hermitian 
symmetry (i.e. V^{u,v) 



V* 



in Eq. 1.6). 



1.6 SZA u,v-space coverage of A1835 (14:01:02.03, +02:52:41.71), a low dec- 
lination cluster (the SZA's latitude is ~ 37°N). See Fig. 1.5 for caption. 
Note that the coverage provided by the inner baselines is nearly as com- 
plete as that in Fig. L5 The coverage on longer baselines, used to con- 
strain the fluxes of sources known to be point-like (using independent ra- 
dio surveys), is less complete. If a source is truly point-like, the it, f -space 
coverage provided by the long baselines does not need to be complete to 
remove it, since a point source has the same magnitude of flux over all 
of UjV-spsice (we generally do not use the SZA to determine whether a 
source is point-like, and need only constrain its flux. See ^6^ for more 

details about modeling unresolved radio sources.) 

SZA Antenna Locations. Antenna 2 is the reference antenna, and is there- 
fore located at the origin. Antennae 6 &; 7 provide 13 long baselines (be- 
tween each other and paired with each of the six inner antennae). These 
long baselines probe small scales, thus aiding point source subtraction. 
Figures L5 & L6 show the n^f -space coverage of two sources observed 

with this array configuration 

Radial distribution of scales probed in -space (y/u'^ + v"^), for coverage 
shown in Fig. 1.5 The distribution is extremely bimodal due to the short 
(~ 0.35-1.3 kA, plotted in green) and long (~ 3-7.5 kA, plotted in blue) 

baselines of the SZA (shown in Fig. 1.7) 

Radial distribution of scales probed in UjV-space for coverage shown in 
Fig. L6 See Fig. L8 for details. Because of the low declination of A1835, 
the projected baselines for this observation are shorter, on average, than 
those for CL1226 (Fig. 



1.8). 



1.10 SZA u,v-spa.ce coverage of CL1226. 9+3332, combining observations from 



both the 30-GHz (blue) and 90-GHz (magenta) instruments. See Fig. 1.5 



for caption. Observations with the SZA using the 90-GHz receivers (see 
^2.1), thus probe finer scales (larger u,v-scales) than the 30-GHz system 
for the same array configuration 
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1.11 Radial distribution of scales probed in u,v-space {\/v? + f 2), for cover- 



age shown in Fig. |1.10[ The ~ 1.3-3 kA gap in the coverage at 30 GHz 
(shown in blue) is filled in by performing complementary observations at 
90 GHz (show in magenta). For a given array configuration and observa- 
tion length, the 90-GHz -coverage can be obtained (approximately) by 
multiplying each u,v coordinate in the 30-GHz -coverage by ~ 3. . . 
1.12 The effective collecting area of the Chandra primary mirror, as a function 
of photon energy. This response dominates the energy response of the 
instrument. It also reduces the plasma emissivity that Chandra effectively 



sees from a high temperature plasma (Fig. 1.13) [19] 

1.13 The X-ray emissivity of a cluster plasma at redshift z = 0.25, redshifted 
to local photon energy range 0.7-7.0 keV, as measured by Chandra (us- 



ing the instrument's effective area shown in Figure 1.12). The plasma 



emissivity was computed using the plasma model of Raymond fc Smith 



(1977) for a range of cluster temperatures and metallicities (see ^5.5.1 ). 



The effective emissivity we measure is reduced by Chandra's efficiency. 



which declines for photon energies > 4 keV (see Fig. 1.12). For plasma 
temperatures < 2 keV, the cluster X-ray emission is dominated by lines 
produced by elements heavier than helium. Since Chandra's sensitivity 
peaks at energies between 1-2 keV, the effective emissivity of a plasma 
with temperature ~ 1 keV has a strong metallicity dependence (compare 
at ~ 1 keV the effective emissivity of a plasma with metallicity Z = 0.9, 
versus that with Z = 0.1) 

2.1 SZA System Overview. The SZA has eight 3.5 m antennae that com- 
municate via fiber optic connections to equipment in a general-purpose 
utility trailer, which is referred to as the "Correlator Trailer." This trailer 
houses the downconverter, correlator, control system computer, and other 
electronics common to the system. The observer commands the control 
system using an interface programmed by Erik Leitch that utilizes the 
SSH (secure shell) protocol 

2.2 Detail of a pair of antennae. The primary mirror diameter is D (solid 
line). The center-to-center antenna separation (long dashed line) is d, 
which when looking at zenith (as depicted here) is also the projected 
baseline length. Separation s (short dashed line) is a distance between 
two arbitrary points on the primaries separated by less than D. A cross- 
sectional representation of the main and secondary lobes of the antenna 
sensitivity pattern is shown above the left antenna, where the first side- 
lobe (secondary lobe) is greatly exaggerated in scale; It was, in reality, 
measured (by James Lamb) to be -25 dB (10^^'^ ~ 1/316*^) less than the 
sensitivity at the center of the primary beam (main lobe) 

2.3 Photo of the inner 6 telescopes (Antennae 0-5) of the SZA. Photo is taken 



from the northeast (roughly along Baseline 2-7; see Fig. 1.7). Antenna 5 



is behind Antennae 1 and 3, and Antenna 4 is behind Antenna 0. 
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2.4 Photo of the outer 2 telescopes of the SZA (Antennae 6 & 7, see 1.7), 
taken while standing near the inner array. Antenna 6 is shown in the left 
panel, and Antenna 7 is shown in the right panel [26] 

2.5 Overview of the SZA antenna optical design. The primary (light blue) and 
secondary (red) mirrors are on-axis reflectors. The tertiary mirror (silver) 
selects the 30 or 90-GHz receiver within the receiver cryostat (dark green). 
The receiver cryostat is located within the larger receiver enclosure box 
(not shown). Image adapted from David Woody [29] 

2.6 Flow chart of the closed-loop chiller system. Warm water is pumped 
through a cool chiller reservoir. The chilled water continues on to the 
sidecab, which is a weather-tight equipment rack containing the antenna 
computers, the motion control system and servos, and electronics power 
supplies, and to the receiver enclosure (described in the text). The warm 
air is blown across the fins of the heat-exchangers, through which the 
cool water is circulated, removing heat. The chilled air returned to the 
surrounding sidecab or receiver box, as appropriate, where it cools the 
electronics. The water, now warm, returns to the pump input, where it 

is recirculated through the system [20] 

2.7 Photo of Receiver Cryostat Electronics test setup. Items 1-10 are compo- 
nents of the 90-GHz receiver, and 11-16 are components of the 30-GHz 
receiver. See §2.3.3 for details on each RF component. The stages of 



the refrigerator (17 & 18) are discussed in ^2.3.4 The copper strapping 

attached to 17 is not the final version 

2.8 90-GHz receiver block diagram, courtesy of Amber Miller. See the photo 



in Fig. 2^, which shows these components: 1. 90-GHz Feedhorn. 2. 
Circular Polarizer & Circular-to-Rectangular Transition. 3. MMIC 
HEMT Amplifier. 4. Isolator. 5. 2'^'^ MMIC HEMT Amplifier, Isolator, 
& High Pass Filter (under bracket). 6. W-band Mixer. 7. Waveguide for 
the incoming signal from 90 GHz tunable LO, the bias-tuned Gunn. 8. 

IF Amplifier. 9. Bandpass Filter. 10. K-band Mixer 

2.9 30-GHz receiver block diagram, courtesy of Amber Miller. See the photo 



in Fig. 2.7, which shows these components: 11. 30-GHz Feedhorn. 12. 
Circular Polarizer & Circular-to-Rectangular Transition. 13. HEMT am- 
plifier (under a mount). 14. Isolator. 15. 27 GHz High Pass Filter. 
16. Mixer. The IF Amplifier in this figure was later moved outside the 
receiver 



viii 



2.10 Receiver Cryostat Thermal test setup. Test hard-soldered strapping (4) 
on the second stage of the cold head (7) is shown. The long copper shims 
(5), used to increase the thermal conductivity between the RF compo- 
nents and the cold head, were replaced with more flexible, nickel-plated, 
braided-copper straps. Also in the photo: 1 is the zotefoam window, 
within the window holder (2). The cryostat case is labeled 3 (the upper 
lid, which is not shown, mates here to make a vacuum seal). The cold 
plate, to which most of the RF components are mounted, is 6. The first 
stage of the refrigerator head is 8. The radiation shield is 9, and the 
exposed part of the refrigerator head is 10. Mylar blanketing can be seen 
filling the space between the cryostat case and the radiation shield. . . . [39] 

2.11 The load curve for the CTI-cryogenics Model 350 Cryodyne refrigerator. 
Figure from Brooks Automation, the company that now owns Helix Tech- 
nology CTI-Cryogenics HQ] 

2.12 Thermal conductivity of OFHC Copper, computed from tables available 
on the NIST website (NIST 2008). These values have been measured in 

the temperature range T = 4 — 300 K 132] 

2.13 Thermal conductivity of Aluminum Alloy 1100, computed from tables 
available on the NIST website (NIST 2008). These values have been 
measured in the temperature range T = 4 — 300 K |42] 

2.14 Thermal conductivity of G-10, computed from tables available on the 
NIST website (NIST 2008). These values have been measured in the 
temperature range T = 4 — 300 K 133] 

2.15 The window holder, as seen looking at the bed where the window is epoxied. 133 

2.16 Cross-section of the window holder, depicting where the first layer of epoxy 

is applied 133 

2.17 Photo of Ebox (lid removed), set up for 30-GHz observations. The Ebox 
houses the modules used to control and process signals from the receiver. 
See text in §2.4.1 for a description of the electronics. Note that the bias- 
tuned Gunn (BTG), the module to control the bias-tuned Gunn (the 



"BTG Mod" in Fig. 2.18), and the IF switch had not been installed at 



the time of this photo, as they were only required for 90 GHz observa- 
tions. They were later installed in locations 12 and 20, respectively. The 
fiber bundle and power cables enter the Ebox at location 16. The refrig- 
erator head (13) for the receiver (located behind the Ebox), the receiver 
enclosure's chiller line (14), the chiller heat exchanger (15), and the TEC 
cooling fans (17) can all be seen above and outside the Ebox. The inter- 
nal Ebox air circulation fans at label 18. A spring that assists in lifting 
the Ebox is labeled 19. The walls of the Ebox are lined with an open-cell 
PVC foam for insulation, which was painted white to prevent degradation 

and flaking due to weathering I3Z] 

2.18 Layout of Modules in the Ebox. Connectors are shown in green |38] 
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2.19 TEC Blower Conceptual Diagram, courtesy Marshall Joy and Georgia 
Richardson. Ambient air from the receiver enclosure is forced through the 
heatsink fins by the blower (see Fig. 2.17, #17). These copper heatsinks 
are thermally coupled to the hot sides of the TECs (which are > 20° C 
hotter than the worksurface side). The TECs diffuse more heat away from 
the worksurface than toward it, thus removing heat from the Ebox. . . . 

2.20 Block Diagram of the SZA Band Downconversion, showing how the sky 
observation frequencies 26.938-34.938 and 90.78-98.78 GHz are split into 
the sixteen digitized bands, each 500 MHz in bandwidth, that are the 
input to the correlator. See description in the text ( ^2.5.1 ). Bands shown 
in red are the side bands that are not used 

2.21 Detailed Diagram of the 90-GHz Receiver Downconversion, showing how 
the sky observation frequencies are brought down to the 1-9 GHz IF 
band (see Fig. 2.8 for details on the receiver RF components). This figure 



accompanies the broader downconversion scheme illustrated in Fig. 2.20 



In the current tuning scheme, sky frequencies 90.78-98.78 GHz are mixed 
((g)) with a 72.28 GHz LO (the bias-tuned Gunn), producing an 18.5- 
26.5 GHz USB product. Sky frequencies < 85 GHz are blocked by the high 
pass filter, so no LSB mixing product is produced (e.g. the band shown 
in red). The bandpass filter passes the 18.5-26.5 GHz USB product to a 
second mixer. The second LO, at 17.5 GHz, mixes with this 18.5-26.5 GHz 
USB product to place the 90-GHz receiver output in the 1-9 GHz IF band. 
This 1-9 GHz IF signal is the input to the back-end electronics common 
to both the 30 and 90-GHz systems. See ^ 2.5.1 for more details [55] 

2.22 Photo of the SZA Correlator and Downconverter, courtesy David Hawkins. [58] 

2.23 Plot of phases on a single baseline (Baseline 0-1) and the corresponding 
downconverter temperatures during the same period. The blue lines are 
the temperatures measured in the downconverters for Antennae (dashed 
line) and 1 (solid line). The black line is the phase of the raw data on 
baseline 0-1, taken while staring at a strong point source, which ideally has 
flat, constant phase (see 4.3.2 for details). Note that thermal variations 
in adjacent digitizers, due to cycling of the A/C, of 0.5°C (peak to peak) 
produced 6° phase variations. Modifications to the A/C removed these 
thermal oscillations on short time scales (see Fig 2.24) [HI] 

2.24 Plot of phases on a single baseline (Baseline 0-1) and the corresponding 
downconverter temperatures during the same period. The blue lines are 
the temperatures measured in the downconverters for Antennae (dashed 
line) and 1 (solid line). The black line is the phase of the raw, uncalibrated 
data on baseline 0-1, taken while staring at a strong point source, with the 
linear slope and mean value removed. Note that this is on a much shorted 
time scale than that of Fig. 2.23 but is long enough to have shown 3 full 
cycles of the A/C if they had persisted. Any residual phase error is both 
negligble and is removed by calibration (see ® [HI] 
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The rectangular window function and its transform. Left: The rectan- 
gular window function w{n) in the time domain. Right: Magnitude of 
the Fourier transform of w{n), \W{k)\. w{n) is unity for n = [0,31], cor- 
responding to 32 discrete samples in time (i.e. the digitized waveform), 
and zero outside this range. |Vr(A;)| is unity for k = 16 and k = —15, and 

zero at all other integer values of k. See discussion in text 

9.022 GHz birdie in the IF, downconverted to 1.022 GHz, sampled in time 
and FFT'd. The blue line is the amplitude of the bandpass that would 
be determined from the SZA spectral data. The birdie does not alias 
cleanly to the center of any channel, and therefore leaks into every other 
channel in the band. The DTFT of the birdie (red curve) is the sum of two 
sine functions, centered at 478 and 522 MHz, (or at ±478 MHz, since the 
transform of a the windowed sinusoid is W{±kQ)). Values of k = [0,16] 
from the DFT are plotted here, corresponding to the positive frequencies 



of Channels through 16 in the SZA correlator (see Table 2.2). The 
band wraps outside k = [0,31] (or k = [—15, 16]), which corresponds to 

frequencies of 0-1000 MHz (or -500-500 MHz) 

Retuned "9.03125 GHz" birdie (placing it at an integer k frequency), 
downconverted to 1.03125 GHz, then sampled and FFT'd. This birdie 



aliases cleanly to the center of the Channel 16 (see Table 2.2), which can 
be excised from the data. The other channels' centers sample the nulls of 
the rectangular window's transform (red curve), which is a summation of 



sines located at ±468.75 MHz. See Fig. 3.2 for more details 

Photo of the closest pair of antennae (Antennae 3 & 5, see Fig. 1.7), 
with the test Eccosorb® wrapped around the feedlegs. Ultimately, the 
Eccosorb® was replaced by crinkled aluminum foil 

System temperatures of a typical SZA cluster observation, computed for 
each band. One representative antenna is shown, since the other seven all 
have similar system temperatures. The y-axis shows Tgyg in K, and the 
X-axis is time in hours from the start of the observation. The dark blue 
point at the start of the observation is the Tgyg of the bandpass calibrator 
( ^4.3.1 ). The black points are the Tgys for the target data (usually a clus- 



ter), and the orange points are Tgyg for the phase calibrator ( ^4.3.2 ). The 



points differ due to the differing columns of atmosphere to each source. 
The minimum in each curve occurs when the source transits (reaches its 
highest elevation in the sky, and therefore the optical depth reaches its 

minimum for the observation) 

Amplitude vs frequency channel of the bandpass calibrator, measured for 
each band of Antenna 5, before calibration. The y-axis is the amplitude 



in Jy, and the x-axis is the channel number. Fig. 4.4 shows the calibrated 
amplitude of the bandpass 
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4.3 Phase vs frequency channel for the bandpass cahbrator, measured for each 
band of Antenna 5, before cahbration. The y-axis is the phase in degrees, 
and the x-axis is the channel number. Fig. 4^ shows the calibrated am- 
plitude of the bandpass [86] 

4.4 Amplitude (vs frequency channel) of the bandpass calibrator for each band 
of Antenna 5, after calibration. See Fig. 4.2, which shows these bandpass 
amplitudes pre-calibration, for caption [ST] 

4.5 Phase (vs frequency channel) of the bandpass calibrator for each band of 
Antenna 5, after calibration. See Fig. |4.3[ which shows these bandpass 
phases pre-calibration, for caption [88] 

4.6 Phase of the point source calibrator versus time, as measured by a short 
baseline (Baseline 0-1). The y-axis is the phase in degrees, and the x-axis 
is time in hours since the start of the track. Band was entirely flagged 
for this track, due to broken digitizers and a lack of spares at the time. 
The magenta points are those flagged due to the automatic flagging, and 
the green points are those flagged by a user-specified script's limits, which 
catches 30° outliers from the underlying, interpolated phase [HO] 

4.7 Phase of the point source calibrator versus time, as measured by a long 
baseline (Baseline 5-6). See Fig. 4.6 for further details. Note the slightly 
larger scatter in the phase on this baseline than on Baseline 0-1 (Fig. 4.6). 
This is due to atmospheric coherence being slightly poorer on long base- 
lines (the antennae are looking through different columns of air, and at- 
mospheric turbulence has a scale size on the order of tens of meters). The 
underlying slope in the phase calibrator is easily determined in this ob- 
servation, as the atmospheric coherence for the example track was typical 

of a clear spring day. [91] 

4.8 Amplitude of the point source calibrator versus time, as measured by a 
short baseline (Baseline 0-1). The y-axis is the amplitude in Jy, and the 
X-axis is time in hours since the start of the track. Since a point source 
is unresolved (the beam is larger than the point source), the amplitude 
in Jy equals that in Jy/beam. This is the amplitude of the calibrator 



corresponding to the phase shown in Fig. 4.6 



4.9 Amplitude of the point source calibrator versus time, as measured by 
a long baseline (Baseline 5-6). This is the amplitude of the calibrator 



corresponding to the phase shown in Fig. |4.7[ See Fig. 4^ for further 

details [93] 

4.10 Antenna-based gain of the point source calibrator versus time, for Antenna 
0. The |/-axis is the gain (ideally 1), and the x-axis is time in hours since 
the start of the track. The scatter seen here, which is simply due to noise 
in the measurement, is smaller than our absolute calibration uncertainty 
(^ 5%, seepuchovej et ah] (|2007|) M 
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4.14 



4.15 



4.16 



4.17 



4.11 Antenna-based phase of the point source cahbrator versus time. One 
antenna must be chosen as a reference in order to compute the antenna- 
based phase of the other antennae. In this plot, Antenna served as 
that reference for all but Bands 6 & 7, which had flagged calibration 
observations (see red points in Fig. 4.6) 

4.12 Orphaned data from Antenna 0. Automatically flagged data in Band are 
shown as green points, and were due to a digitizer hardware problem (and 
a temporary lack of replacement digitizer boards). Data orphaned by gaps 
in the calibrator data in Bands 1, 6, & 7 are shown in red (see Fig. 4.11). 
Useful data are represented as blue points. For simplicity, target data 
are plotted with zero phase (rather than plotting the noisy distribution of 
raw target data phases, since we are only trying to determine which data 
are not bracketed by calibrator observations). The antenna-based phases 
of bracketing calibration observations are shown as black X's (which are 
identically zero in bands where Antenna was the reference; see Fig. |4.1l" ). [971 

4.13 Orphaned data from Antenna 5. Data orphaned by gaps in the calibrator 
data in Band 2 are shown in red (see Fig. 4.7). See Figure 4.12 for details. 
Note that Antenna 5 was not the reference antenna for any band in the 
antenna-based phase calculation, so none of its phases are identically zero. 
rms of target data, taken on short Baseline 0-1. The theoretical pre- 
diction, based on the measured Tsys, is shown in black. Flagged data - 
already caught by other steps (see Figures 4.6 & 4.12 in particular) - are 
plotted in red and magenta. Newly flagged data are in green. Useful, 

unfiagged data are in blue 

rms of target data, taken on long Baseline 5-6. Flagged data - already 
caught by other steps (see Figures 4.7 & 4.13) - are plotted in red and 

magenta. See Figure 4.14 for more details IIOII 

Amplitude of visibility data (in Jy, see ^1.3[ ) for a short baseline (0-1), 
taken on the target source (typically a cluster). Blue points are unfiagged 
data. Red and magenta points are flagged data caught by previous steps 
in the data calibration (see Figures 4.6 & 4.12), while green points may 
indicate data flagged in this step as outliers or from previous steps (e.g. 
Bands 1, 2, & 14 contain some green points that are not outliers, flagged in 
Fig. 4.6). The SZE flux from the cluster and the fluxes of point sources in 

mJy level; they are therefore not noticeable 



the cluster field are on the 

in the raw data plotted here (which are noise dominated) 

Amplitude of visibility data in Jy (see ^1.3) for a long baseline (5-6) 



mi 



taken on the target source. See Fig. |4.16 Note that long baselines do 



not typically measure cluster scales, but the raw target data here and 
in Fig. 4.16 have similar amplitudes. Also note that some green, flagged 
points (e.g. in Band 13) were flagged in Fig. 4.7, not because they are 
out hers 
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4.18 Dirty maps of the short (upper panel, < 2 kA) and long (lower panel, 
> 2 kA). X and y axes are the map coordinates in degrees. The colors 
represent the signal in Jy/beam (the flux detected within the beam formed 
by each baseline). See text for more details 11061 

5.1 Geometry for the line of sight integral of a spherically-symmetric model. IllOl 

5.2 The above figure shows fits to the density profiles of 11 nearby, relaxed real 
clusters, as well as to the average density profiles of 16 clusters simulated 
using adiabatic (red) and cooling + star-formation (blue) physics. V06 



motivate the use of a nine-parameter density model (Eq. 5.12) to model 

these density profiles. Figure from N07 11161 

5.3 The above figure shows fits to the temperature profiles of 11 nearby, re- 
laxed real clusters (in green, magenta, and cyan), as well as to the average 
temperature profiles of 16 clusters simulated using adiabatic (red) and 
cooling + star-formation (blue) physics. All profiles are scaled to their 
best-fit rsoo values. V06 used the eight-parameter temperature model 



(Eq. 5.13) to capture the details of the real cluster temperature profiles. 
Figure from N07 EH] 

5.4 The above figure shows the X-ray-derived pressure profiles of 5 nearby, 
relaxed real clusters with Tx > 5 keV, as well as to the average 3-D, 
radially-averaged pressure profiles of 16 clusters simulated using adiabatic 
(red) and cooling -|- star-formation (blue) physics. The black lines are the 
best fit generalized NEW profiles for each type of cluster plotted. Figure 
from N07 HH 

5.5 This figure shows the density of the standard /3-model (red) and SVM 
(blue) using the steepest outer slope V06 allowed, e = 5. This was cho- 
sen to illustrate the largest disparity the radially-average outer cluster 
slope V06 considered "realistic." We chose typical parameters for the 
components common to the two profiles (namely, we chose Tc = 100" and 
P = 0.7, and set the normalization pq = 1 for simplicity). For the SVM, 
we show a typical massive cluster's scale radius for the slope to steepen, 

r, = 500" CMl 

5.6 The ratio of the thermal energy content within a given radius to the total 
as r — )■ oo, plotted for the fit to A1835 of the N07 pressure profile. The 
vertical line shows rsoo — 360" ~ 1.4 Mpc 11381 

6.1 X-ray image of A1835, showing it to be relaxed. The X-ray analysis of 
A1835 relies on a single Chandra ACIS-I exposure, with 85.7 ks of good 
time (unflagged exposure time). The pixels of the ACIS-I detector are 
binned to be 1.968" on a side. The X-ray image shown here is smoothed 
with a Gaussian that is 2 pixels in width (for display purposes only). See 



Table 6^ for more details on the X-ray observation. The inner 100 kpc 
(core) and all detected X-ray point sources were excluded from the X-ray 
surface brightness and spectroscopic analyses 11421 
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6.2 X-ray image of CL1226, showing this high redshift cluster to be approx- 
imately circular on the sky, and thus apparently relaxed. This image 
combines two Chandra exposures. The inner 100 kpc (core) and all de- 
tected X-ray point sources were excluded from the X-ray analysis. . . . 11431 

6.3 X-ray image of A1914, showing it to be disturbed. This image combines 
two Chandra exposures. A large subclump, measured to be hot by M08, 
is located to the east (left) of the cluster center. This subclump, the inner 
100 kpc (core), and all detected X-ray point sources were excluded from 

the X-ray analysis 11441 

6.4 Cleaned SZA 30-GHz map of A1835, showing the detection to have a 
high significance of ~ 22-cr at the peak. This image, made with Difmap 
using only the short baseline data (~ 0.35-1.5 kA), is for presentation 
purposes only; we fit our data directly in -space, and do not use any 
interferometric mapping to determine cluster properties. Note that for a 
given u,v-space coverage, any unresolved structure in an interferometric 
SZE map of a cluster essentially "looks" like the synthesized beam; this 
determines the effective resultion of the observation. The synthesized 
beam is depicted in gray in the lower left corner. Contours are overlaid 

at 2-0" intervals 11471 

6.5 Cleaned SZA 30-GHz map of CL1226, showing the detection to have a 
significance of ~ lO-a at the peak. See caption of Fig. 6A for further 
details EH] 

6.6 Cleaned SZA 90-GHz map of CL1226, showing the detection with this 
higher resolution instrument to have a significance of ~ 12-0" at the peak. 
This observation was included since 30-GHz data alone could not constrain 
the radial profile of this high-redshift, small angular extent cluster (see 
Figures 1.10 and 1.11). See Fig. 6.4 for further details, noting the short 
baselines of the 90-GHz instrument are ~ 1-4.5 kA 11491 

6.7 Cleaned SZA 30-GHz map of A1914, showing the detection to have a high 
significance of ~ 18-0" at the peak. See Fig. |6.4|for detail 11501 
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6.8 



Radially-averaged SZE model fits to A1835 in -space, from the jointly- 
fit N07-I-SVM and isothermal /3-model. The upper panels show the real 
(left) and imaginary (right) components of the visibilities, radially-averaged 
to be a function of (m, v) radius, and rescaled to units of intrinsic, line-of- 
sight integrated Compton y; I plot here the frequency- independent quan- 
tity Y{u,v)d\ = V{u,v) d\/g{x)Io, where each band is scaled appropri- 
ately before binning, and the angular diameter distance is computed using 
the assumed ACDM cosmology. The lower panels show the reduced of 
the fits to the data for the chosen binning. The black points with error 
bars (1-cr) are the binned Y{u, v) d\ data, with the point source models 
first subtracted from the cluster visibilities. The blue, solid line is a high 
likelihood N07 model fit, while the red, dashed line is a similarly-chosen 
fit of the /3-model. For the available data points, both SZE models fit 
equally well (see lower panels, which shows the for each model is indis- 
tinguishable). However, note that as the -radius approaches zero kA - 
where there are no data to constrain the models - the /3-model predicts 
a much higher integrated Compton y than the N07 model. For cluster 
data centered on the phase center of the observation (with both cluster 
and primary beam sharing this center), the mean imaginary component 
would equal zero. Since these model fits include the primary beam, which 
is not necessarily centered on the cluster, the small but non-zero imag- 
inary component in the upper right panel is expected (note the smaller 
units for the y-axis of the right hand plot) (See, e.g. Reese et al. 2002, for 



comparison) 11551 



6.9 Radially-averaged SZE model fits for CL1226 in n,f-space. See Figure 6.8 



for additional caption details. Note that the 30-GHz -space coverage 



alone does not probe a sufficient range of cluster scales (see Figures 1.5 



and 1.11 ) to determine where the SZE signal falls to zero. This results in 
poor constraints on the cluster's radial profile when using 30-GHz data 
alone, and is a result of this high-redshift source being relatively compact 
on the sky (compared to A1835 and A1914; note that r25oo is on the 
order of 1 arcminute, as shown in Table 6.4, ^6.7). I therefore included 



90-GHz SZA data (green) in the joint SZE+X-ray fit, since the 90-GHz 
instrument was designed to complement the w^f -coverage provided at 30- 
GHz (see Figures 1.10 & 1.11). Again, the small but non-zero imaginary 
component is expected, but note the smaller units for the plot of the 
imaginary component of the radially-binned Y{u,v) 



[T56] 



6.10 Radially-averaged SZE model fits for A1914 in fi^w-space. See Figure 6.8 



for details. Note that the slightly larger imaginary component in the upper 
right panel could be due to cluster asymmetry, since A1914 is disturbed 
and elliptical (see Fig. 6.7). Note that the reduced of this fit is no 
larger than those for the other two clusters (Figures 6.8 & 6.9) 
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6.11 



6.12 



6.13 



6.14 



6.15 



6.16 



X-ray surface brightness profile fits to the clusters. The vertical dashed 
line denotes the 100 kpc core cut. The blue, solid line is the surface 
brightness computed using a high-likelihood fit of the N07-I-SVM profiles 
(analogous to the SZE fits plotted in Figures 6.8-6.10), while red, dot- 
dashed line is the surface brightness fit of a /3-model. Both model lines 
include the X-ray background that was fit simultaneously with the cluster 
model (i.e. the plotted lines are the superpositions of each set of the 
background and cluster models, which were fit simultaneously to the X- 
ray imaging data). The black squares are the annularly-binned X-ray 
data, where the widths of the bins are denoted by horizontal bars. The 
vertical error bars are the 1-a errors on the binned measurements. Arrows 
indicate r25oo and t^qq derived from the N07+SVM profiles (see ^6.7). . 
Pe(r) for each set of models fit to each cluster. The pressure from the 
jointly-fit N07+SVM is plotted in blue with vertical hatching. Pressure 
constrained by the SZE fit of the isothermal /3-model is plotted using red, 
dot-dashed lines; note that the isothermal /3-model's shape is constrained 
by X-ray imaging data, and the only unique parameter to the SZE data in 
this fit is the central decrement (SZE normalization, see ^5.2.1). Pressure 
derived from the density and temperature fits of the V06 profiles in the 
independent X-ray analysis is shown in black with grey shaded regions. 
See text in §6.6| for details. The vertical, black dashed line shows r25oo 
derived from the N07+SVM fits, while the magenta dashed line is for t^qq 

(see Tables 6.4 and 6.5) 

Pgiisir) for each set of models fit to each cluster. Colors and hne styles 
are the same as in Fig 6.12 See text in ^6.6 for details. Note that the 



MB 



vertical, black and magenta dashed lines show r25oo and rsoo, respectively, 

derived from the N07+SVM fits 

Te(r) for each set of models fit to each cluster. Colors and hne styles are 

Note that the isothermal /S-model's constant 
constrained by the X-ray data, is plotted using red, dot- 
See text in 96.61 for details. Note 



USD 



the same as in Fig |6.12 

r) 



dashed lines and dark red shading, 
that the vertical, black and magenta dashed lines show r25oo and rsoo, 

respectively, derived from the N07-I-SVM fits 11631 

Yios computed within 6' (~ 1.4 Mpc, which is ~ v^qq for this cluster) for 
fits to SZA observation of A1835. The bold, black contours contain 68% 
and 95% of the accepted iterations to the jointly-fit Chandra + SZA data, 
while the thinner, blue contours are those for fits to SZA data alone. 
Note that the vertical, black and magenta dashed lines show r25oo and 

rsoo, respectively, derived from the N07+SVM fits 11661 

y{R) - Compton y (integrated along the line of sight, Eq. 1.2 ) as a function 
of sky radius R. Colors and line styles are the same as in Fig 6.12 
36.61 for details. 



See 

text in y6.6| for details. Note that the vertical, black and magenta dashed 
lines show r25oo and v^qq, respectively, derived from the N07+SVM fits. . 
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6.17 



6.18 



6.19 



6.20 



6.21 



Flos ^ Compton y integrated over sky radius R. Colors and line styles 
are the same as in Fig |6.12[ See text in §6.6| for details. Note that the 



vertical, black and magenta dashed lines show r25oo and rsoo, respectively, 

derived from the N07+SVM fits 

Mgas - the gas mass integrated within a spherical volume defined by cluster 
radius r. Colors and line styles are the same as in Fig 6.12 See text in 



Mtot - the total mass estimated assuming hydrostatic equilibrium at clus- 
ter radius r. Colors and line styles are the same as in Fig 6.12 See text 



7.1 Constraints on Mt, 



and /, 



gasv 



[T69] 



for details 1170] 



in §6.6| for details 11711 

1-D histograms of the N07+SVM jointfit estimates, for A1835, of Mtot and 
Flos normalized by their respective median values, Mtotis the cyan region 
with a dashed outline, and Y\os is the vertically hatched region with a solid 
black outline. Both M^t and Yjos are computed within a fixed radius of 
9 = 360". The derived Fios, which scales with integrated SZE flux, has a 
more tightly constrained and centrally peaked distribution than that of 
Mtot, as Mtot is sensitive to the change in slope of the pressure profile 



(see Eqs. 5.16 & 5.32). Since the N07 model was developed primarily to 
recover Yint from SZE observations, it is useful to note how well this model 

performs in this capacity 11711 

/gas - the gas mass fraction, computed as Mgas/Mtot at cluster radius r 
for each accepted iteration in the Markov chain. Colors and line styles 



are the same as in Fig 6.12 See text in ^6.6 for details 11731 



using the N07+SVM combined with 



the spectroscopically-measured temperature, using Tx = T^h to fit A1835. 
The spectroscopically-measured Tx used here is the same as that used in 
the isothermal /3-model analysis presented in Chapter |6| The N07+SVM 
results that include X-ray spectroscopy are plotted using red, dot-dashed 
lines and red shading. The N07+SVM results without spectroscopic con- 
straints (i.e. the resulted detailed in Chapter [6]) are plotted in blue with 
vertical blue hatching. Results derived from the density and temperature 
fits of the V06 profiles, from the independent X-ray analysis, are shown 
in black with grey, shaded regions. Both panels show that including spec- 
troscopic information in fits of N07-t-SVM profile tightens constraints and 
improves the already remarkable agreement between it and the indepen- 
dent X-ray analysis 11771 
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Chapter 1 
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Introduction 

"What's so amazing that keeps us stargazing? What do we think we might see?" 
- Kermit the Frog ( The Rainbow Connection, written by Paul Wilhams and Kenneth 
Ascher) 

1.1 Clusters of Galaxies 

The detailed expansion history of the Universe and the growth of large scale structure 
are two of the most important topics in cosmology. Clusters of galaxies are the 
largest gravitationally-bound systems in the Universe, and thus provide a unique 
handle on cosmic expansion and structure formation. Measurements of the growth 
of structure, traceable using large cluster surveys, provide critical clues to the nature 
and abundance of dark matter and dark energy]^ 

At the time of this writing, a low density, cold dark matter (CDM) cosmology, 
dominated by a dark energy or cosmological constant (A) component, is heavily fa- 
vored]^ This is called ACDM cosmology, where 'A' refers to dark energy's contri- 
bution, expressed as Q^, to the total energy density of the Universe, Q. Cold dark 
matter is a form of non-baryonic matter travelling much slower than the speed of 
light, which - along with the baryonic component fib ^ provides matter's energy den- 
sity contribution ^m- In general, any quantity Qx is the ratio of x's energy density 
to the critical energy density it would take to close the Universe. 



^ For recent reviews of how cluster studies can be used to constrain cosmology, see the Dark 
Energy Task Force (DETF) report ( jAlbrecht et al.||2006[ ) and |Rapetti fc Allen] ( [2007 1 . 

^It is unknown what this dark energy component is, whether it truly acts as a "cosmological 
constant," or whether it evolves over time. 
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ACDM cosmology started to become the favored cosmology in the 1990's, when it 
displaced the then-favored - and now inappropriately-named - Standard Cold Dark 
Matter (SCDM) picture of the Universe. SCDM holds that = ^tot = I {^a = 
0). ACDM cosmology reconciles two results that cannot be resolved in the SCDM 
paradigm: the flatness of the Universe combined with the low value for the universal 
matter density f^M- 

In the mid-1990's, X-ray gas mass measurements began to demonstrate that Qm < 
l|^ Because galaxy clusters collapse out of representatively large, comoving volumes 
of the Universe (~ 10 Mpc), we assume their baryonic/dark matter ratio approaches 
the universal value. Using X-ray measurements of galaxy clusters, one can obtain 
both the hot gas mass (Mgas) and an estimate of the cluster total mass Mtot, and 
use these to compute the gas mass fraction /gas = Mgas/Mtotj^ Provided that the 
bulk of a cluster's baryons are in the hot gas probed by X-ray observations, we can 
approximately equate the gas mass fraction /gas ~ ^gas/^M ~ ^b/^M (within a factor 
of ~ 2, this is well-supported by the simulations). 

The combination of /gas measurements with constraints from Cosmic Microwave 
Background (CMB) measurements and Big Bang Nucleosynthesis (BBN) predictions 
yields an estimate for ^m- Since CMB measurements, coupled with BBN theory, 
independently constrain Qx^h"^ (where h is the local Hubble constant divided by 
100 km s~^ Mpc~^), one can use cluster gas fraction measurements to approximately 
constrain Qu- By accounting for the composition of the Coma cluster - adding up 



the observable mass in the hot gas and optically luminous stars - White et al. ( 1993 ) 



argued that the measurements implied a low value for in this way. Later, David 



et al. (1995) applied this method to many more clusters observed by ROSAT, and 



Optical measurements of the stellar and total masses of galaxies and galaxy clusters also showed 
this, though I focus here on the dominant baryonic component of the largest collapsed structures in 
the Universe - the X-ray emitting gas in galaxy clusters. 

'^How /gas is obtained from X-ray and X-ray-|-SZE observations is discussed in J 5.5.5 
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used this to estimate that ~ 0.1-0.2. 

Furthermore, X-ray cluster surveys - such as those provided by the ROSAT All- 
Sky Survey (RASS) and the Wide Angle ROSAT Pointed Survey (WARPS) - pro- 
vided independent indications that we live in a low density universe with Qm ~ 
0.2-0.3. This low density was inferred from the lack of strong cluster number evolu- 
tion, which would have been seen if = 1 (see, e.g. [Mushotzky fc Scharf 1997), as 



the matter density strongly affects when clusters form (see, e.g. Holder et al. 2000 



Haiman et al. 2001, and references therein). 



Whilst cluster measurements consistently implied we live in a low density universe, 
direct evidence for dark energy was provided by those who set out to measure the 
deceleration of the Universe's expansion, using type la supernovae (SNIa) as "stan- 
dard candles" (objects with a known luminosity). Surprisingly, their measurements 



indicated that the expansion of Universe is accelerating (Riess et al. 1998), imply- 
ing it is dominated by some form of dark energy. Soon after, measurements of the 
CMB made by the BOOMERanG, TOCO, and MAXIMA strongly constrained the 
Universe to be spatially flat, a result that has been confirmed by WMAP and other 
CMB measurements since. With no curvature component [Q]^ = 0), "flatness" means 
that the sum of the angles within a triangle is 180°, and that Qm + = 1- 

With strong evidence now in place for ACDM, and many more recent results 
confirming this, a "concordance cosmology" with = 0.7, = 0.3, and = 
can be defined. This cosmology is assumed for most of the results presented in this 
thesis, and is in agreement with the parameters published in the first year results 



of the Wilkinson Microwave Anisotropy Probe (Spergel et al. 2003, WMAP), which 
jointly fit data from a number of independent tests (e.g. their own and finer scale 
probes of the CMB, weak lensing, cluster counts, galactic velocity field, and HST's 



Key Project) described in Spergel et al. (2003). Defining a "concordance cosmology" 
provides a convenient framework for both observers and theorists. As cosmological 
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constraints are refined, any results assuming this cosmology can be scaled to reflect 
the updated parameters]^ 

Recently, a number of studies have used clusters alone to constrain the dark energy 



component of the Universe (see Allen et al.| 2004 LaRoque et al. 2006 Allen et al. 



2007 Vikhlinin et al. 2008, and references therein). This is done primarily using 
X-ray measurements of the intracluster medium (ICM), the hot (> 10"^ K) gas that 
comprises the majority of the baryons in a galaxy cluster. Assuming a theoretically- 
motivated functional form for the gas fraction, one can solve for the cosmological 
parameters that force the data to fit the expected fgas{z) evolution for redshift z. 

In this thesis, I explore the joint constraints provided by two independent probes 
of the ICM, using these to measure /gas, Mgas, Mtot, and other cluster astrophysical 
properties. In the next section, I describe the primary tool used here to constrain the 
properties of the ICM - the Sunyaev-Zel'dovich effect. 

1.2 The Sunyaev-Zel'dovich Effect 

The Sunyaev-Zel'dovich effect (SZE) is a unique probe of the hot gas in clusters, 
since it does not depend on emission processes, unlike X-ray emission or the broad 
majority of luminous astrophysical processes routinely measured by astronomers. The 
SZE arises by inverse Compton scattering of CMB photons off the hot electrons in 



the ICM (depicted schematically in Fig. 1.1). This leaves a spectral signature on the 



CMB that is independent of redshift, since at higher redshifts the CMB is both denser 
in photon number and less redshifted in energy]^ 

There are two separable components to the Sunyaev-Zel'dovich effect: the kinetic 

^This assumes, of course, that new results remain consistent with ACDM. 

^ For a given state of the electrons in the ICM at any redshift, the same number fraction of CMB 
photons are inverse Compton scattered by the same fraction of the photon energy. Therefore, the 
(1 + z)'* dependence in how the energy of the CMB is redshifted, which also maintains its blackbody 
spectrum, ensures the SZE spectral signature on the primary CMB is constant relative to the CMB. 
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Figure 1.1 Image adapted from van Speybroeck (1999). 



SZ effect, or KSZ, due to the line-of-sight proper motions of clusters, and the thermal 
SZ effect (TSZ, simply called 'SZE' in the chapters that follow). I do not discuss the 
KSZ here, since we could not measure it with the Sunyaev-Zel'dovich Array (SZA, the 
instrument presented in this thesis; see Chapter [2]). The thermal SZE is measurable 
as a spectral distortion of the CMB, computed 



AT 



fix) y, 



XI) 



CMB 



where AT is the temperature distortion of the CMB, which has temperature Tcmb, 
y is called the "Compton y parameter" and is integrated along the line of sight, and 



f{x) (given in Eqs. 1.4 and 1.5) contains the SZE frequency dependence. 



The line-of-sight Compton y parameter is computed 



^ ln,{e)T,{i)di, 



'1.21 



where ks is Boltzmann's constant, Ut is the Thomson scattering cross-section of the 
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electron, rrie is the mass of an electron, c is the speed of light, and ne{i) and Te{i) are 
the electron density and temperature along sight line i. The Compton y parameter 
has a linear dependence on electron pressure Pg when the ideal gas law, Pg = ksUeTe, 
is assumed: 

y=^[Peii)di. (1.3) 
The classical, non-relativistic frequency dependence /(x) of the SZE is 

f{x)=x(^^]-A (1.4) 



where the dimensionless frequency x is 

x = hu/kBTcMB- (1-5) 
Here h is Planck's constant, and u is the frequency of the observation. The classical 



spectral dependence f{x) given in Eq. 1.4, and relativistic corrections to it - necessary 



due to the high thermal velocities of electrons in the ICM - are plotted in the upper 



panel of Fig. 1.2 over a broad range of frequencies. Here I used the calculations of 



Itoh et al. (1998), which correct the SZE f(x) out to fifth-order. The SZE decrement 



< 218 GHz becomes an increment at > 218 GHz. A detail of this crossover is 



shown in the lower panel of Fig. 1.2, which shows that the relativistic velocities of 



the electrons, due to the cluster's temperature, shift the precise location of SZE null. 



Finally, Fig. 1.3 shows a detail of the SZE decrement at the frequencies the SZA can 
probe. 

1.3 Inter ferometry Overview 



Monochromatic light from an arbitrarily-shaped, spatially-incoherent aperture propa- 
gates via Fraunhofer diffraction to become the spatial Fourier transform of the light's 
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Figure 1.2 SZE spectral dependence /(z^), plotted as a function of frequency (see 
Eq. 1.5). The upper panel shows /(z/) over a broad range of frequencies, while the 
lower panel shows a detail of the null in the SZE spectrum. The relativistically- 
corrected /(z/) for a cluster with temperature kTe = keV (black, solid line) reduces 
precisely to the classical frequency dependence (Eq. 1.4[ ), since electrons with no 
temperature are not moving at relativistic random velocities. The other lines show 
how the relativistically-corrected /(z/) departs from the classical behavior for higher 
temperature electrons. The relativistic corrections to /(z/) shown here are computed 
out to fifth-order using the equations provided in Itoh et al. (1998). The classical 
SZE spectrum has a null at z^ ~ 217.5 GHz, above which the SZE signal becomes an 



increment. Higher temperature electrons require relativistic corrections (Itoh et al. 



1998) to the classical SZE frequency dependence, which shift the null to higher fre- 
quencies. A high temperature cluster would have a non-negligible thermal SZ effect 
at the classical null (217.5 GHz). 
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Figure 1.3 SZE frequency dependence below 100 GHz. See Fig. L2 for caption. 
Below the null (see lower panel of Fig. 1.2), we often use the term "SZE decrement" 
to refer to the strength of the thermal SZ effect, which is negative at low frequencies. 
The lower panel shows the fractional deviation from the classical SZE for clusters at 
higher temperatures, due to the relativistic electron velocities in the ICM. Treating 
the thermal SZE at 30 GHz from a massive 10 keV cluster as classical introduces a 
~ —3.6% bias in quantities derived from the SZE fits (i.e. the line-of-sight electron 
pressure in Eq. [1.2 would be underestimated by this amount, since the strength of 
the SZE would be overestimated by the classical calculation). 



intensity pattern, after traveling a distance many times its wavelength. Radio astro- 
nomical interferometry takes advantage of this simple fact, where the astronomical 
source serves as the arbitrary aperture. An interferometric array measures the Fourier 
transform of the source's spatial intensity distribution, probing angular scales \/d' 



(in radians) for wavelength A and projected baseline separation d' (see Fig. 



1.4). 
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Figure 1.4 Example baseline formed by a pair of antennae. The dish separation is c?, 
while the projected baseline as seen by the source is c?'. The distance to the source 
(represented by a cloud) is labeled i?//. Note the scale of the broken lines, representing 
this astronomical distance -R//, is not accurate; in reality, the lines of sight are nearly 
parallel and point at the same location on the source. The time delay T is equal to the 
difference in distances to the source, from each antenna, divided by the speed of light 
c. By delaying the signal measured by the antenna on the right by time T, we ensure 
that the same astronomical wavefront is used in the correlation of the two signals. 
Tracking an astronomical source requires both the physical pointing of each antenna 
toward the source in the far field, and the implementation of instrumental delays (T) 
that ensure each antenna measures the same wave front as the Earth rotates. 

1.3.1 Criteria and Assumptions for Interferometry 

The conditions necessary to take advantage of astronomical interferometry constitute 



the "van Cittert-Zernike theorem" (see Thompson et al. 2001, for a derivation and 



many further details). I summarize the necessary criteria here: 

• The source must be in the far field, meaning that its distance i?// ^ {d'Y jX 



(illustrated in Fig. 1.4), where d' is the longest projected baseline in the array 



for a given observation, and A is the wavelength at which the observation is per- 
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formed. A large distance is essential for statistically independent emission (or 
scattering events, in the case of the SZE) from the source to become coherent 
plane waves. For the Sunyaev-Zel'dovich Array, this means we must observe 
sources more than 360 km away, a condition met by any astronomically inter- 
esting source]^ In meeting the far field condition, for a sufficiently narrow band, 
we can treat incoming light from the source as a series of plane waves. 

• As a corollary to the above, we assume there are no sources in the near field, 
so we are only observing sources in the far field. 

• The source must be spatially incoherent, meaning that two arbitrary points 
within the source must not have statistically correlated emission. Since two 
points within a source are physically separated and entirely independent, emis- 
sion from these two points will not in general be coherently emitted. 

We also take advantage of one more approximation: the small angle approxima- 
tion. While this is not necessary for the van Cittert-Zernike theorem, it simplifies 
the analysis of interferometric observations, as it allow us to treat the data as a 2-D 



Fourier transform (discussed in the next section, ^ 1.3.2). 

As I discuss later in this section, the largest radial scales probed by the SZA 

are on the order of ~ 5' (1.45 x 10^^ radians). The small angle approximation is 

therefore well-justified for single, pointed observations made with the SZA (as opposed 

to mosaicked observations). We therefore treat the spatial intensity pattern of the 

source as if it were truly in a plane perpendicular to the line of sight. We call this the 

"image plane," and justify this by pointing out that the sources we observe are far 

enough away that the distance to any part of the source is negligibly different from 

the radial distance to the source's center. 

^ The Moon is ^ 3.8 x 10^ km away, and is much closer than the closest source we have observed, 
Mars. 
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At a given instant in time, all antennae must measure a single, monochromatic 
wavefront for the astronomical signals to be correlated. The signals are measured 
within spectral channels (a "channel" is a smaller range of frequencies within some 
band) that are small compared to the central frequency of each band, assuring each 
signal is approximately monochromaticj^ To ensure the condition that we are measur- 
ing a single plane wave front, the differing distances from each antenna to the source 
are corrected by adding (computationally) an adjustable time delay to each antenna's 
signal. This adjustable delay is the difference in path lengths from the astronomical 



source to each antenna, divided by the speed of light (labeled T in Fig. 1.4). The path 
length cT changes throughout the course of an observation, as the source traverses 
the sky. 

In radio interferometry, computing and applying the proper delays to the signal 
from each antenna, as the Earth moves, is called "fringe tracking." The Cartesian 
location (x, y) = (0, 0) (typically north-south and east- west angular offsets) is assigned 
to the point in the sky called the "pointing center." Mechanical tracking keeps this 
point (e.g. a cluster's center) in the center of each antenna's primary beam; the 
pointing center is also the "phase center" for the observation, since the adjustable 
delays T are computed in order to precisely compensate for the different distances to 
this point. These differing distances are based simply on the geometry illustrated in 



Fig. 1.4 



1.3.2 Probing Sources in Fourier Space {u,v-space) 

In the context of astronomical interferometry, Fourier space is often referred to as 
"ttjV-space," since u and v are the Fourier conjugates of image space coordinates x and 



^The SZA observes at sky frequencies ~ 30 and 90 GHz, and breaks each of its sixteen ~ 500 MHz 



bands into fifteen usable 31.25 MHz channels (for more details, see { 2.5.1 1. Each channel is therefore 
^ 1/1000 the sky frequency. 
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SZA M,v-coverage of CL 1226+3332 (at 30 GHz) 






u (kX) 



Figure 1.5 SZA u,v-space coverage of CL1226.9+3332 (12:26:58.0, +33:32:45.0), a 
cluster that passes near zenith for the SZA (the SZA's latitude is ~ 37°N). Each 
blue point represents a data point's u,v-space location, using frequency independent 
units of kX (distance divided by wavelength, divided by 1000). Note that there are 2 
groupings of points: those ~ 0.35-1.3 kA from the center, due to the short baselines 
of the compact inner array, and those ~ 3-7.5 kA from the center, due to the long 
baselines formed with the outer antennae (i.e. each baseline formed with Antennae 



6 or 7; the antenna layout of the SZA is shown in Fig. 1.7). The SZA was designed 



to provide this broad, uniform coverage in u,v-space on the shorter baselines, while 
simultaneously providing higher resolution coverage with longer baselines. The data 
locations in M,f-space exhibit inversion symmetry (i.e. {u,v) = —{u,v)) because the 
visibilities are the transform of real data in image space. The visibilities therefore 
exhibit Hermitian symmetry (i.e. Vi,{u,v) 



in Eq. 1.6). 
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SZA M,v-coverage of A 1835 (at 30 GHz) 
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Figure 1.6 SZA u,v-space coverage of A1835 (14:01:02.03, +02:52:41.71), a low dec- 
lination cluster (the SZA's latitude is ~ 37°N). See Fig. 1.5 for caption. Note that 



the coverage provided by the inner baselines is nearly as complete as that in Fig. 1.5 



The coverage on longer baselines, used to constrain the fluxes of sources known to 
be point-like (using independent radio surveys), is less complete. If a source is truly 
point-like, the -space coverage provided by the long baselines does not need to 
be complete to remove it, since a point source has the same magnitude of flux over 
all of u,v-spa.ce (we generally do not use the SZA to determine whether a source is 



point-like, and need only constrain its flux. See ^6.2 for more details about modeling 
unresolved radio sources.). 
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SZA Antenna Locations 
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Figure 1.7 SZA Antenna Locations. Antenna 2 is the reference antenna, and is there- 
fore located at the origin. Antennae 6 & 7 provide 13 long baselines (between each 
other and paired with each of the six inner antennae). These long baselines probe 



small scales, thus aiding point source subtraction. Figures |1.5| & |1.6| show the u,v- 
space coverage of two sources observed with this array configuration. 



yjj When calibrated against an astronomical source of known flux, an interferometric 
array's output can be expressed as the flux (in Janskiee^°) at each point probed in 
u,v-space, for each frequency probed by the instrument. For a given integration time 
(typically ~ 20 sec for the SZA), each baseline and band of the radio interferometric 
array produces one binned point in u,v-space, called a "visibility."!^ The visibilities 
behave collectively as a Fourier transform of the spatial intensity pattern /^(x, y) (flux 



^ Since we arc performing interferomctry on sources that are far (compared to the observation's 
wavelength), we can ignore the z spatial component along the line of sight, as well as its transform 
w. This is equivalent to stating that our sources lie in the image plane. 

lOln S.I. units, 1 Jy = 10^26 w • m^^ . Hz-i. 

-•^^For the 28 baselines and 16 bands of the SZA, four hours of useful, on-source observation time 
produces « 322, 560 independently measured visibilities. This is after the channels of each band are 
binned into one visibility per unit time. 
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^ Radial M,v-Space Coverage of CL1226 (at 30 GHz) 
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Figure 1.8 Radial distribution of scales probed in u,v-space (a/m^ + v'^), for coverage 



shown in Fig. 1.5 The distribution is extremely bimodal due to the short (~ 0.35- 



1.3 kA, plotted in green) and long (~ 3-7.5 kA, plotted in blue) baselines of the SZA 



(shown in Fig. 1.7). 



^ jq4 Radial M,v-Space Coverage of A1835 (at 30 GHz) 
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Figure 1.9 Radial distribution of scales probed in -space for coverage shown in 



Fig. L6 See Fig. L8 for details. Because of the low declination of A1835, the 
projected baselines for this observation are shorter, on average, than those for CL1226 
(Fig. pi). 
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per unit solid angle) at frequency z/, and is calculated (Thompson et al. 2001): 



-j2iT{ux+vy) 



a/1 - ,t2 - 



: dx dy. 



;i.6) 



Here A^{x,y) is the spatial sensitivity of each antenna, the central lobe of which is 



called the primary beam (illustrated in Fig. 2.2). Each antenna's beam was mapped 



by Ryan Hennessy and Mike Loh, and was found to be well-approximated by a sin- 
gle, circularly-symmetric Gaussian with a half-power radius of ~ 4.9'. Because our 
observations are on the scales of arcminutes, we can make the approximation that 
a/1 — x^ — y^ ~ 1, where x and y in are radiansj^ Doing this, Eq. 1.6 simplifies to 
the 2-D Fourier transform of the spatial intensity pattern multiplied by the beam: 



A,{x,y)h{x,y) e-^'^^'^^+^yUx dy. 



;i-7) 



Figures 1.5 & 1.6 show the u,v-space coverage of the SZA interferometer during 
two typical cluster observations. This coverage was provided using the array configu- 



ration shown in Fig. 1.7, performing observations with the 30-GHz receivers. Figures 



L8]fc 1.9 respectively show the distributions of u,v-Tadu probed by the «,f-spacings 



plotted in Figs. L5 & L6 The incompleteness of this coverage necessarily means 
information on those scales not probed is missing; an interferometer filters informa- 
tion on scales it cannot access. In fact, the SZA was designed specifically to probe 
cluster (arcminute) scales. Since the scales to which we are sensitive are proportional 
to X/d in radians, long baselines of a 30-GHz SZA observation (~ 3-7.5 kA) probe 
angular scales 6 oc X/d = 0.46-1.15', while short baselines of that same observation 
(~ 0.35-1.3 kA) probe angular scales 6 (x X/d = 2.6-9.8'. The SZA 30-GHz system 
is essentially not sensitive to radial scales larger than ~ 5' (the 90-GHz system, dis- 



^For example, observations that probe a radial distance on the sky of 5' implies y'^x^ 



,2 < 



0.0015 radians, and therefore the term -^/l — — > 0.9985. This justifies our use of the small 
angle approximation for single, targeted cluster observations used to probe arcminute scales. 
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Combined 30+90 GHz SZA M,v-coverage of CL 1226+3332 
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Figure 1.10 SZA u,v-space coverage of CL1226. 9+3332, combining observations from 



both the 30-GHz (blue) and 90-GHz (magenta) instruments. See Fig. 1.5 for caption 



Observations with the SZA using the 90-GHz receivers (see ^2.1), thus probe finer 



scales (larger -scales) than the 30-GHz system for the same array configuration. 

cussed below, probes scales one-third these sizes). In this way, the SZA isolates the 
small (AT in Eg. 1 1.1 [ typically peaks ~ 10 mK) cluster signal from the relatively large 



background (e.g. the atmospheric and instrumental noise discussed in ^2.3.1, as well 
as the 2.73 K primary CMS). 

By including SZA observations performed using the 90-GHz receivers, using the 



same array configuration shown in Fig. |1.7[ the gap at ~ 1.3-3 kA in the f -coverage 
at 30 GHz can be filled. Since the wavelenth is ~ 3 times shorter at 90 GHz than it 
is at 30 GHz, each baseline at 90 GHz is effectively ~ 3 times longer in A (number 



of wavelengths). The combined 30+90 GHz coverage is shown in Fig. 1.10, while the 
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^ Combined Radial w,v-Space Coverage of CL1226 
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Figure 1.11 Radial distribution of scales probed in u,v-space (a/m^ + v'^), for coverage 



shown in Fig. 1.10 The ~ 1.3-3 kA gap in the coverage at 30 GHz (shown in blue) 



is filled in by performing complementary observations at 90 GHz (show in magenta). 
For a given array configuration and observation length, the 90-GHz -coverage can 
be obtained (approximately) by multiplying each u,v coordinate in the 30-GHz u,v- 
coverage by ~ 3. 



corresponding histogram of radial -scales probed by combining the observations is 



shown in Fig. 1.11 The short baselines of the 90-GHz system, ranging ~ 1.2-3.8 kA, 
constrain cluster signals in the intermediate angular scales 9 oc X/d = 0.9-2.9'. 



1.4 X-ray imaging of Galaxy clusters 

In this thesis, I use X-ray imaging data taken with the Chandra X-ray Observatory}^ 



extending the work of Reese et al. (2002), LaRoque et al. (2006), Bonamente et al. 



(2006), and others to combine X-ray imaging data with new SZE data taken with the 



SZA. 



At X-ray wavelengths, emission from a cluster is predominantly due to thermal 



^ ^http : / / cxc . harvard . edu/ 
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Chandra Effective Area (0.7-7.0 keV) 
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Figure 1.12 The effective cohecting area of the Chandra primary mirror, as a function 
of photon energy. This response dominates the energy response of the instrument. It 
also reduces the plasma emissivity that Chandra effectively sees from a high temper- 



ature plasma (Fig. 1.13). 



processes that occur in the same gas that produces the SZE. The two thermal pro- 
cesses by which the ICM emits X-rays are bremsstrahlung and line emission (see, 



e.g. Longair 1998 Sarazin 1988), both of which depend on collisions between pairs 



of particles within the gas (i.e. two particles are involved in the emission process). 
The X-ray emission therefore scales as number density-squared. In contrast, the SZE 



depends linearly upon electron pressure (Eq. 1.2). These two ways of probing the 
ICM therefore complement each other. 



X-ray imaging data are sensitive to the surface brightness Sx (in cts arcmin ^ s 



2 



s 



X 



where ne(i) and Te{i) are the electron density and temperature along sight line 
^"^Free-free emission due to electron collisions. 
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Broadband X-ray Emissivity (0.7-7.0 keV, z=0.25) 
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Figure 1.13 The X-ray emissivity of a cluster plasma at redshift z = 0.25, redshifted 
to local photon energy range 0.7-7.0 keV, as measured by Chandra (using the instru- 
ment's effective area shown in Figure 1.12). The plasma emissivity was computed 



using the plasma model of Raymond & Smith (1977) for a range of cluster tempera- 



tures and metallicities (see |: 5.5.1). The effective emissivity we measure is reduced by 



Chandra^s efficiency, which declines for photon energies > 4 keV (see Fig. 1.12). For 
plasma temperatures < 2 keV, the cluster X-ray emission is dominated by lines pro- 
duced by elements heavier than helium. Since Chandra^s sensitivity peaks at energies 
between 1-2 keV, the effective emissivity of a plasma with temperature ~ 1 keV has 
a strong metallicity dependence (compare at ~ 1 keV the effective emissivity of a 
plasma with metallicity Z = 0.9, versus that with Z = 0.1). 

Aee(Te, Z) (iu cts cm^ s~^) is the X-ray emissivity measured by the instrument within 
the energy band used for the observation, z is the cluster's redshift, and Z is the 
plasma's metallicity. Metallicity accounts for elements heavier than helium, where 
Zq = 1 is defined to be the elemental abundance measured in the solar atmosphere, 
and the abundances in a cluster are measured relative to the solar abundance (see 



5.5.1 for a more detailed discussion of elemental abundances and metallicity). For 



T > 3 keV, the X-ray emissivity of the plasma depends weakly on temperature. 
Before accounting for the instrument's response (discussed below), Aee(Te) oc T^/^ 



since the X-ray emission is bremsstrahlung-dominated at these energies. 
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In addition to the already-weak temperature dependence of the plasma emissivity, 
the "effective emissivity" of the plasma as measured by Chandra is reduced by the 
instrument efficiency, which declines for photon energies > 4 keV. Chandra's effi- 
ciency is dominated by the "effective area" of the primary mirror for a photon of a 
given wavelength; the typical energy response of the Chandra ACIS-I CCD is plotted 



in Fig. 1.12 The effective emissivity measured by Chandra in a given energy band 



accounts for the instrument's efficiency at that energy, and is plotted in Figure [1.13 
for a range of plasma temperatures and metallicities for a plasma at z = 0.25j^ To 
obtain the total number of counts per pixel in an exposure, the result of Eq. |1.8|must 



be multiplied by the exposure time and the field of view (in arcmin^ ) of a pixel0 

1.5 Structure of the Thesis 

I describe the SZA instrumentation in Chapter |2| and discuss how we addressed 
challenges encountered during SZA commissioning observations in Chapter [3j In 
Chapter |4| I discuss SZA data calibration and the final data product used in cluster 
analyses. 

In Chapter |5| I discuss the joint modeling of SZE-I- X-ray data, and present the 
new models used. I also discuss our data fitting routine and how we derive cluster 
parameters of interest from the models. In Chapter [HI I apply the models to real 



observations and discuss the results. Finally, in Chapter 6.8, I present my conclusions 



and some ideas for extensions to these modeling techniques. 

The effective area Aeff{E) = e{E) x A for the efficiency e{E) for photon energy E (local to 
the instrument) of an aperture with physical area A. The "effective emissivity" could be computed, 
then, as Aee,eff = J^^ Aee,E ^eff{E) dE , for an X-ray exposure in the local energy range E1-E2, 
given a physical plasma emissivity of Aee,_B redshifted to local photon energy E. 

-•^^There are additional complications. First, the effective area changes over the field of view of an 
X-ray instrument. We use an exposure map, which gives the effective area each pixel sees at 1 keV 
(near the peak in Chandra's effective area) to correct for off-axis effects. We also correct the X-ray 



image for the quantum efficiency of the CCD. See Bonamente et al. (20041 for more details on the 
X-ray analysis. 
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Chapter 2 

The Sunyaev-Zel'dovich Array 

2.1 Overview of the Sunyaev-Zel'dovich Array 

The Sunyaev-Zel'dovich Array (SZA), an interferometer composed of eight 3.5 meter 
telescopes, is located at the Owens Valley Radio Observatory, in Big Pine, Califor- 
nia. Its coordinates are 37°14"02" latitude and 118°16"56" longitude, and it is at an 
altitude of 1222 meters above sea level. The Owens Valley is a desert, and therefore 
provides suitable atmospheric conditions for performing centimeter and millimeter- 
wave observations for most of the year. 



The SZA has a digital correlator with 8 GHz of bandwidth (see Hawkins et al. 2004 
for details on the correlator), and each antenna is equipped with two wideband receiver 
systems, capable of observing from 27-36 GHz (in the Ka band, hereafter referred to 
as the "30-GHz" band) and from 85-115 GHz (in the W band, referred to as the "90- 



GHz" band). See Fig. 2.1 for a broad overview of the SZA system. The large, 8-GHz 
receiver and correlator bandwidth provides the SZA with the high sensitivity required 
to detect rapidly cluster signals. This wide bandwidth also provides the ability to 
probe a wide range of u,v-space (Fourier conjugate of image space) simultaneously, 
meaning it is sensitive to a wider range of angular scales than a comparable instrument 
with a smaller bandwidth. 

The SZA receivers contain high electron mobility transistors (HEMT) for high 
gain, low noise amplification of incoming signals. The 30-GHz receivers contain 
the same HEMT amplifiers that were in the Degree Arcminute-Scale Interferometer 



(DASI) (see [Leitch et al.|2005D , in the same configuration used for the OVRO/BIMA 
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Figure 2.1 SZA System Overview. The SZA has eight 3.5 m antennae that commu- 
nicate via fiber optic connections to equipment in a general-purpose utihty trailer, 
which is referred to as the "Correlator Trailer." This trailer houses the downcon- 
verter, correlator, control system computer, and other electronics common to the 
system. The observer commands the control system using an interface programmed 
by Erik Leitch that utihzes the SSH (secure shell) protocol. 



SZE receivers (see Carlstrom et al. 1998 Carlstrom et al. 2000). The 90-GHz re- 



ceivers utilize HEMTs implemented through monolithic microwave integrated circuit 
(MMIC) technology, which integrates four HEMTs into each MMIC block. 

The SZA was designed to detect and probe clusters at intermediate and higher 
{z > 0.2) redshifts. The relatively small (3.5 m) primary mirrors and their short 



focal lengths (see ^2.2) provide two advantages for SZE observations of clusters: a 



large primary beam, which scales as X/D for observational wavelength A and primary 
mirror size D, and the ability to closely pack the antennae without resulting in inter- 
antenna collisions. The full width half maximum (FWHM) size of the truncated 
Gaussian primary beam, where the sensitivity falls to -3 dB of the peak, is ~ 10.7' 
at the center of the 30 GHz band. This is necessary so that objects on arcminute 
scales, such as clusters are intermediate redshifts, are relatively unattenuated by 
the primary beam. The scales probed by an interferometer are determined by the 
projected separation distance d (the distance between the centers of two antennae as 
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d 



Figure 2.2 Detail of a pair of antennae. The primary mirror diameter is D (solid line). 
The center-to-center antenna separation (long dashed line) is d, which when looking 
at zenith (as depicted here) is also the projected baseline length. Separation s (short 
dashed line) is a distance between two arbitrary points on the primaries separated 
by less than D. A cross-sectional representation of the main and secondary lobes 
of the antenna sensitivity pattern is shown above the left antenna, where the first 
sidelobe (secondary lobe) is greatly exaggerated in scale; It was, in reality, measured 
(by James Lamb) to be -25 dB (10~^'^ ~ 1/316*'^) less than the sensitivity at the 
center of the primary beam (main lobe). 

seen from the source), and scale as \/d. Because the SZA routinely observes clusters 
all the way to the antenna shadowing limit of 3.7 meters (i.e. the projected antenna 
spacing is nearly as small as the dish size), it can in principle probe radial angular 
scales larger than ~ 4.8' at the central frequency of 30.938 GHz. This is because, for 
a projected baseline of length d, there are points on the primary mirrors separated 
by distances s G [d — D,d + D] that are not entirely attenuated by the primary beam 



(see Fig. 



2.2). 
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Figure 2.3 Photo of the inner 6 telescopes (Antennae 0-5) of the SZA. Photo is taken 



from the northeast (roughly along Baseline 2-7; see Fig. 1.7). Antenna 5 is behind 
Antennae 1 and 3, and Antenna 4 is behind Antenna 0. 




Figure 2.4 Photo of the outer 2 telescopes of the SZA (Antennae 6 & 7, see 1.7), 
taken while standing near the inner array. Antenna 6 is shown in the left panel, and 
Antenna 7 is shown in the right panel. 
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Six of the eight SZA antennae, forming fifteen basehnesQ comprise the closely- 
packed "inner array," which is sensitive to ~ 1-5' scales. The inner array is shown 



in Fig. 2.3 Two outer antennae (Fig. 2.4), which are identical to the inner six, form 



thirteen more baselines (i.e. one baseline between the two outer antennae, and six 



baselines formed between each outer antenna with each of the inner six; see Fig. 1.7). 
These outer antennae provide the ability to fit simultaneously any unresolved (< 20") 
radio sources (hereafter "point sources") in the cluster field, which could otherwise 
mask the SZE decrement at 30 GHz. 



As discussed in [1.3, when observing an astronomical source, in addition to the 



physical pointing of each antenna that keeps the source in the main lobe of the 



antenna beam pattern (see Fig. 2.2), the delays necessary for fringe tracking must 
be computed. Fringe tracking corrects the phases for the instantaneous projected 
baseline changes due to geometry, as the source moves through the sky, relative to 



the array (see Fig. 1.4). Phase changes can also be due to effects that are more difficult 
to calculate. For example, variations in the lengths of the fiber optics and cables that 
carry the signals can change the phases of the signals from each antenna. Phase shifts 
due to properties of the electronics themselves can also act as path length differences 
(at a particular frequency). These instrumental phase effects are particularly sensitive 
to temperature fluctuations. 

In this chapter, I provide an overview of the components that comprise the SZA. 
Emphasis is placed on thermal and mechanical aspects of the instrument, as this was 
the focus of the instrumentation component of this thesis. Thermal stability and the 
consequences of poor thermal regulation are discussed as appropriate. 



^Thc number of baselines iVbi is computed: iVbi = N(N — l)/2 
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2.2 Telescope Optics 

The antennae of the SZA are designed as on-axis, altitude-azimuth telescopes with 



small primary mirrors that have short focal lengths (see Fig. 2.5). Their design allows 
for a compact array configuration with the antennae spaced as close as 1.2D, where D 
is the diameter of the primary, without ever colliding. The panels of the primary mir- 
rors were machined by Jerry Forcier of Forcier Mechanics using a conventional, com- 
puter numeric controlled (CNC) mill, programmable using standard computer-aided 
design/machining (CAD/CAM) code. Since the mill operates in standard Cartesian 
coordinates, the machining process provided the mirrors with a scalloped surface. 
The tool and toolpath used to machine the primaries were chosen so that the scal- 
loping scatters higher frequency light (e.g. IR and optical), which would otherwise 
damage the instrument if, for example, the sun were focused on the receiver]^ The 
surface was measured by Marshall Joy to have an rms roughness of 25 /im, meeting 
the design requirements necessary for the 1.3 mm upgrade the SZA may see as part 



of CARMA (see e.g. Woody et al. 2004 Scott fc Pound 2006). 



The secondary mirrors are convex paraboloids which slow the focu^ before it 
enters the receiver enclosure. The receiver enclosure is simply a large, weather-proof, 
thermally-regulated box that holds the receiver cryostats, the back-end electronics 
(which are further enclosed in a smaller, thermally-regulated electronics box), the 
tertiary mirror, an ambient calibrator load, and various support electronics. In the 
fall and spring, the receiver enclosure is cooled by air circulating through closed-cycle 
air/air heat-exchangers. In the summer, the enclosure is cooled by a refrigerated 
water /antifreeze solution, which is pumped through a closed liquid/air heat-exchanger 
(similar to a car radiator, containing refrigerated fluid) between the air in the box 



^Not all telescopes can point directly at the sun without damaging the receivers; the SZA, 
however, can. 



^See Fig. 



2.5 



which shows that the rays are being quickly focused by the primary; the secondary 



mirror extends the focal length. 
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Secondary Mirror 



■primary Mirror 




Figure 2.5 Overview of the SZA antenna optical design. The primary (hght blue) and 
secondary (red) mirrors are on-axis reflectors. The tertiary mirror (silver) selects the 
30 or 90-GHz receiver within the receiver cryostat (dark green). The receiver cryostat 
is located within the larger receiver enclosure box (not shown). Image adapted from 
David Woody. 



and the chilled fluid (see Figure 2.6). In the winter, the box is heated by ceramic, AC- 
powered heaters. In addition to protecting equipment from the elements (e.g. dust, 
harsh sunlight, and rain), the receiver enclosure serves as the first step in achieving 
thermal stability, and its stability helps reduce expansion and contraction of the 
section of fiber optics closest to the receiver. 

The signal enters the receiver enclosure through a microwave transparent (has a 
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Water Chiller 



cool water 



Pump 



warm water 



Cool air 



Receiver Enclosure 
Heat Exchanger 



Receiver Enclosure 




Sidecab 
Heat Exchanger 



SideCab 




Figure 2.6 Flow chart of the closed-loop chiller system. Warm water is pumped 
through a cool chiller reservoir. The chilled water continues on to the sidecab, which 
is a weather-tight equipment rack containing the antenna computers, the motion con- 
trol system and servos, and electronics power supplies, and to the receiver enclosure 
(described in the text). The warm air is blown across the fins of the heat-exchangers, 
through which the cool water is circulated, removing heat. The chilled air returned to 
the surrounding sidecab or receiver box, as appropriate, where it cools the electronics. 
The water, now warm, returns to the pump input, where it is recirculated through 
the system. 

very small absorption/reflection coefficient - less than a few percent), weather-proof 
window made of 75 fim mylar and 1-inch thick zotefoam (a nitrogen-filled, polyethy- 
lene foam). Zotefoam, a commercially-available foam used in housing construction, 
was chosen because it is highly transparent at microwave frequencies. It is a closed-cell 
foam that does not permit air or water diffusion, and is very robust both chemically 
and mechanically. This mylar layer in front of the zotefoam of the receiver enclosure's 
window provides an additional - and easily replaceable - layer of protection from dust 
and water, which cling to the porous surface of the zotefoamj^ Thin mylar, such as 



*A porous surface develops as the zotefoam weathers and degrades in sunUght. 
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that used for the windows, is also microwave-transparent. 

The next mirror in the optical path - directly after the window to the receiver 
enclosure - is an adjustable tertiary mirror: a concave, off-axis paraboloid that bends 
the optical path by 90° into the receiver. This programmable tertiary selects the 
receiver into which the signal is focused. 

2.3 Receivers 

2.3.1 Receiver Noise Considerations 

The primary goals of receiver design are the attainment of high gain and low noise. We 
characterize noise in terms of the "receiver noise temperature" Trx, which is derived 
in this section. 

The power emitted in the Rayleigh- Jeans tail of a blackbody is approximately 
linearly proportional to its temperature. For two thermal sources at temperatures Thot 
and Tcoid, we measure output powers Phot and Pcow, which are linearly proportional to 
the source temperatures as long as the gain of the receiver system also remains linear 
over that range. Our measurement includes an additional, constant amount of poweij^ 
due to the noise in the receiver, which can be attributed to a hypothetical thermal 
source at temperature Tj.^. This receiver noise temperature therefore characterizes 
the output power that we would measure in the absence of any input power. We 
define the receiver y factor (the ratio of the measured output powers): 



^This is the zero-intercept in a hypothetical graph of measured output power versus input black- 
body temperature. 



y = 




Thot + ^r: 
-^cold ~)~ 



rx 



rx 
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Solving for Tj-x, we have: 



rp ^hot yTcold /2 2\ 



Contributions to Tj-x can include noise from the receiver cryostat windows, amplifiers, 
mixers, waveguides, and the feed horns. These contributions are all affected by the 
physical temperatures of the components. 

The total contribution Tj-x from the receiver to the noise is calculated in terms of 
each component's gain G (which can be less than unity) for n stages: 

T2 T, 
Trx = Ti + — + „ „ + . . . + — (2.3) 



In Eq. |2.3[ the noise due to each successive stage is reduced by the combination 
of gains before it {Gi ■ G2 ■ ■ ■ ■ ■ Gn-i)- The first components in a receiver - the 
window and the feed horn - are chosen to have low loss and be well-matched to 
the incoming signal, respectively, in order to reduce their noise contributions. In 
addition to being well-matched to the incoming signal, the feed horns are kept at 
low temperature, which reduces their noise contribution. The first stage amplifier 
typically dominates the overall receiver noise temperature Tj.^. We therefore optimize 
the system to have a high-gain, low-noise amplifier at the first (and at every, if 
possible) stage of amplification. Higher current, higher noise amplifier stages follow 
this stage, providing most of the overall amplification of each receiver system. 

For the overall noise contribution to a signal from a telescope - including contri- 
butions from the back-end electronics, the atmosphere, optical elements in the signal 
path to the receiver, and of course the noise from the receiver itself - we define the 



system temperature Tgys, where the y factor in Eq. 2.1 is instead measured at the 
output of the entire system, and T^ys is scaled to above the Earth's atmosphere (i.e. 
we include the noise contribution of the atmosphere with the instrumental noise in 



one term; see ^ 4.3.4 ). The dominant contributions - in order of importance - to the 
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total system temperature Tgyg are therefore the receiver noise T^, the sky (~ 10-20 K 
at zenith), noise scattered from the grounc^ (which is called "spillover," and is ~ 4 K 
at very low elevations, which we avoid when observing), loss due to the windows in 
the optical path, poor coupling with the feed horn, and the noise contribution from 
first stage of the HEMT or MMIC amphfier. 



2.3.2 Measuring T^^ 

In the laboratory, we used a piece of Eccosorb^Has a "hot load" (a thermal load with 



T = Thot in Eq. 2.2). The "cold load" (a thermal load with T = T^om in Eq- 2.2) 
was another piece of Eccosorb® at liquid nitrogen's (abbreviated LN2) boiling point 
(~ 77 K at 1 atmosphere of pressure). This is achieved by keeping the blackbody 
material in a bath of LN2, and removing it for the measurement (for only a few 
seconds, so it retains evaporating LN2 throughout the measurement). 

On the SZA telescopes, the hot load is a mechanized version of the same ambient 
blackbody used in the lab, mounted to the front of the receiver cryostat. However, a 
true cold load with an LN2 bath would be impractical in our system. Instead, careful 
tracking of atmospheric temperature and humidity allow a model to be developed for 
the power contribution from the sky, and therefore the sky serves as the "cold load" 
for calibration purposes. 

Under ideal observational conditions the atmospheric contribution to the noise 
scales as the optical depth. The optical depth scales roughly as sec (6*) for angle 
6 measured from zenith. This sec (6') dependence only holds if assuming a plane- 
parallel atmosphere; this is a gross simplification, and the true calculation used in 

^Dirt is essentially a blackbody at ^ 300 K. Power emitted by the ground can scatter off the 
feed legs holding the secondary mirrors. This small amount of noise was measured to be ^ 4 K at 
elevations < 12°. We typically limit our observations to sources above a 30° horizon. 

^Eccosorb® is a microwave blackbody made by Emerson & Cuming Microwave Products. De- 
tails about this product can be found here: |http : / /www . eccosorb . com/ europe/english/ page/] 
|63/eccosor b 
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data calibration does not make this assumption. To reduce the noise contribution 
from the atmosphere, we typically observe sources with elevations > 30° above the 
horizon. 



2.3.3 Receiver RF Components 

Within the cryostat there are two complete, independent receiver systems (see photo 



in Fig. 2.7). These receivers are capable of observing at sky frequencies 27-36 GHz 



(the "30-GHz receiver") and 85-115 GHz (the "90-GHz receiver"). Each receiver (see 



Figures 2.9 and 2.8) begins with a corrugated feed horn (^1 and #11 in Fig. 2.7) 



that couples the incoming signal to a circular polarizer (#2 and #12 in Fig. 2.7). A 
circular-to-rectangular transition, which unfortunately rejects one linear polarization, 
couples the circular polarizer to a rectangular (linearly polarized) waveguide. The 
components that follow - starting with the HEMT and MMIC amplifier^ (#3 and 



#13 in Fig. 2.7) - are specific to each receiver system. I present the 30-GHz receiver 



first, but note that it is similar in design to the 90-GHz system. 



30-GHz Receiver RF Components 

In the 30-GHz receiver, a single HEMT amplifier with four stages (i.e. four discrete 



HEMTs) amplifies the signal. An isolator (#14 in Fig. 2.7), which reduces reflections 



coming back from the filter, follows the HEMT amplifier. Next, a high-pass filter 



#15 in Fig. 2.7) eliminates unwanted signals below 27 GHz, and helps to determine 



the range of the intermediate frequency (IF) signal, which is the sky signal mixed 



down to 1-9 GHz (see [2.5.1 for a description of the SZA downconversion scheme). 



We refer to the 1-9 GHz band as the "IF band," and provide the corresponding sky 
frequencies of each band within the IF band (see Table 



2.1 



^2.5.1). The IF band 



®Their definitions are repeated here, for convenience. HEMT is "high electronic mobiUty tran- 
sistor," and MMIC is "monoUthic microwave integrated circuit." 
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Figure 2.7 Photo of Receiver Cryostat Electronics test setup. Items 1-10 are compo- 
nents of the 90-GHz receiver, and 11-16 are components of the 30-GHz receiver. See 
^2.3.3 for details on each RF component. The stages of the refrigerator (17 & 18) 



are discussed in ^2.3.4 The copper strapping attached to 17 is not the final version. 
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85-115 GHz Receiver Block Diagram 



IF Amplifier 
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Bandpass fitter 
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Figure 2.8 90-GHz receiver block diagram, courtesy of Amber Miller. See the photo in 
Fig. 2/7, which shows these components: 1. 90-GHz Feedhorn. 2. Circular Polarizer 
& Circular-to-Rectangular Transition. 3. 1*^* MMIC HEMT Amplifier. 4. Isolator. 
5. 2^^^^ MMIC HEMT Amplifier, Isolator, & High Pass Fiher (under bracket). 6. 
W-band Mixer. 7. Waveguide for the incoming signal from 90 GHz tunable LO, the 
bias-tuned Gunn. 8. IF Amplifier. 9. Bandpass Filter. 10. K-band Mixer. 



26 - 36 GHz Receiver Block Diagram 



Mixer 

27 - 36 GHz RF 
Zt GHz LO 
1 -9 GHz IF 



IF Amplifier 
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27 GHz High Isolator HEMT circ to 
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Figure 2.9 30-GHz receiver block diagram, courtesy of Amber Miller. See the photo 



in Fig. 2/7, which shows these components: 11. 30-GHz Feedhorn. 12. Circular Po- 
larizer & Circular-to-Rectangular Transition. 13. HEMT amplifier (under a mount). 
14. Isolator. 15. 27 GHz High Pass Filter. 16. Mixer. The IF Amplifier in this 
figure was later moved outside the receiver. 
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is the band in which all the back-end electronics (^2.4) common to both the 30 and 



90-GHz receivers work. A mixer (#16 in Fig. 2.7), using a 35.938 GHz local oscillator 



(LO) as its reference, mixes sky frequencies 26.938-34.938 GHz down to the 1-9 GHz 



IF band (the LO is discussed in ^2.4.1, while the full downconversion from sky to IF 



bands is discussed in ^2.5.1). 



The high pass filter determines the upper end of the IF band by attenuating signals 
below 27 GHzj^ The high filter attenuates unwanted harmonics related to the LO, 
since the LO is actually the 4th harmonic of an 8.972 GHz oscillator (discussed in 



^2.4) mixed with a 50 MHz phase reference frequency, and other harmonics of the 



8.972 GHz oscillator can be present in the system. I discuss RF contamination - and 



how we eliminated it - in J: 3.1 



90-GHz Receiver RF Components 



The 90-GHz receiver (see Fig. 2.8) is similar to the 30-GHz receiver, with a few 



exceptions. Notably, it has two MMIC amplifiers (#3 and #5 in Fig. 2.7), each of 



which contains four HEMTs (i.e. four stages), since these higher-frequency MMIC 
HEMTs have less gain than their 30-GHz counterparts. The MMICs are separated by 



an isolator (7^4 in Fig. 2.7), which reduces coupling between the gain stages. Feedback 
between the stages could otherwise lead to oscillations, since coupling from reflections 
is a form of positive feedback (which is intrinsically unstable). 



There are two mixers in the 90-GHz system. The first (7^6 in Fig. 2.7) mixes the 
sky signal with a tunable LO that, along with an 18.5-26.5 GHz bandpass filter (#9 



in Fig. 2.7), selects the range of observable sky frequencies for the 90-GHz receiver. 



The second mixer (#10 in Fig. 2.7), using a fixed-frequency, 17.5 GHz dielectronic 



resonant oscillator (DRO), brings the signal down to the 1-9 GHz IF band. An 



additional IF amplifier (#8 in Fig. 2.7) is present in the 90-GHz system for two 

noting that none of the 



^Note that the channel at 26.938 GHz is not used; See Table 
channels at the band edges are used. 



2.2 



2.5.2 



38 



reasons - the MMIC HEMTs have less overall gain than their 30-GHz counterparts, 
and mixers are lossy components that split the input signal power among the output 
mixing products. Since the output of the 90-GHz receiver is handled by the same 
electronics setup as that which handles the 30-GHz receiver's output, the two systems 
need to have comparable output power levels. 



2.3.4 Thermal Considerations for the Receivers 



The signal enters the receiver cryostat through a zotefoam and mylar window. The 



RF components in the receiver are sensitive to temperature, and therefore must be 
kept cool (11-20 K) to reduce noise. This necessitates the placement of the amplifiers 
and other noise-sensitive components in a vacuum-sealed cryostat. 

The cryostat is cooled by a CTI-cryogenics Model 350 Cryodyne refrigerator. 



which offers two stages of cooling: the first (#8 in Fig. 2.10) is a higher power 



warmer stage that runs at ~ 55 K under a load of about 10 W, and the second (#7 
in Fig. 2.10[ ) is a lower power, cooler stage that runs at ~ 12 K when the load is 



about 1 W (see Fig. 2.1l[ ) p] The cryostat is therefore built in two stages as well: 
A warmer radiation shield, made of aluminum (> 90% reflectivity in the microwave 
and far infrared, where the thermal emission at < 100 K peaks), and a cooler stage 
within the shield, where the most sensitive components of the system are mounted. 
The cold head of the refrigerator interfaces with the two stages using copper 



strapping. Copper was chosen for its high thermal conductivity (see Fig. 2.12). For 
the coldest stage, it was necessary to machine and hard solder (using silver solder) a 
piece that ensured good thermal coupling on each end without mechanically stressing 



the refrigerator as the system is warmed and cooled (labeled #4 in Fig. 2.10). The 



^"The charge carriers in a semiconductor are less likely to be in the conduction band due to 
thermal motion as T — > K, which means most of the current in the amplifier is due to the applied 
field from the in coming signal. 
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Figure from 



http : //www . brooks . com/documents . cf in?docnmentID=4932 
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Figure 2.10 Receiver Cryostat Thermal test setup. Test hard-soldered strapping (4) 
on the second stage of the cold head (7) is shown. The long copper shims (5), used 
to increase the thermal conductivity between the RF components and the cold head, 
were replaced with more flexible, nickel-plated, braided-copper straps. Also in the 
photo: 1 is the zotefoam window, within the window holder (2). The cryostat case 
is labeled 3 (the upper lid, which is not shown, mates here to make a vacuum seal). 
The cold plate, to which most of the RF components are mounted, is 6. The first 
stage of the refrigerator head is 8. The radiation shield is 9, and the exposed part of 
the refrigerator head is 10. Mylar blanketing can be seen filling the space between 
the cryostat case and the radiation shield. 
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Model 350 Two-Stage Cryodyne ^^frigerator 

typical Perfortnance (60 Hz] 




FIRST STATE TEMPERATURE (K) 



Figure 2.11 The load curve for the CTI-cryogenics Model 350 Cryodyne refrigerator. 
Figure from Brooks Automation, the company that now owns Helix Technology CTI- 
Cryogenics. 

thermal link from the second stage of the refrigerator cold head to the cold plate is 
made by four thin copper shims, each of which has a small kink (zig-zag) that helps 
provide strain relief. The thermal power P conducted across a thermal link between 
temperatures Ti and T2 can be calculated, 

P=y[ \{T)dT (2.4) 
L Jti 

where k is the thermal conductivity in (Wm~^K~^). 

Indium film (0.01" thick) provides a soft, conductive layer between the nickel- 
plated copper refrigerator head and the copper strapping. Indium is also used in the 
interface between the strapping and the cold plate. The four oxygen-free high-purity 
(OFHC) copper shims that conduct heat from the cold plate to the second stage of 
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the refrigerator head are each 0.05" thick, 2" wide, and about 2" long, yielding an 



area-to-length ratio A/ L = 4 x 0.05" ~ 0.005 m (see Eq. 2.4). At 12 K, OFHC copper 



has a high thermal conductivit}- of ~ 7500 (Wm~^K~^), so a load of 1 W produces 
a temperature differential of only ~ 0.03 K along the length of the strapping to the 
cold plate. The measured differential across this copper strapping agreed with this 



calculation to within the resolution of the Lakeshore temperature sensors we used. 
Aluminum alloy 1100, the common alloy of which the receiver parts are machined - 
including the cold plate where the sensitive receiver electronics are mounted - has 
a far lower conductivity at low temperatures than copper (aluminum alloy 1100 has 



a conductivity k ^ 10 - 100 (Wm^^K-^); see Woodcraft (2005); NIST (2008) and 



Fig. 2.13). Even with a thickness of 0.25", the aluminum cold plate exhibits a ~ 3 K 
thermal differential. Braided copper strapping was therefore used to link the most 
sensitive components - the 30-GHz HEMT and two 90-GHz MMIC HEMT blocks - 
directly to the refrigerator head. 

The stages of the cryostat are separated by stand-offs made of a machinable fiber- 
glass/epoxy composite known as "G-10" or "garolite." It has a thermal conductivity 
~ 0.1-0.2 (Wm^^K-i) at temperatures from 10-100 K, and ~ 0.6 (Wm-^K-^) at 



room temperature, making it a good insulator (see Fig. 2.14). The six stand-offs for 
the radiation shield, each 0.813 in (2.06 cm) long and with a cross-section of 0.062 in^ 
(0.40 cm^), contribute about 0.8 W to the total thermal load on the radiation shield. 
Five stand-offs are used to separate the cold plate from the radiation shield. These 
stand-offs are each 0.5 in (1.25 cm) long and have the same cross-sectional RrGci cis 
the first set of stand-offs; they contribute roughly 50 mW to the cold plate. 

A multi-layered blanket of aluminized mylar fills the space between the inner wall 
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See NIST website, 



http : //www. cryogenics .nist .gov/MPropsMAY/material 



properties , 



htm 



^^The temperature resolution of the readout was 0.125 K, and the time resolution was ^ 1 s. 
Given this temperature resolution and the fact that the refrigerator cycles ±0.5 K in ~ 2 s, the 
temperature at each end of the strapping changes too rapidly to measure precisely this differential. 
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Figure 2.12 Thermal conductivity of OFHC Copper, computed from tables available 



on the NIST website (NIST 2008). These values have been measured in the temper- 



ature range T = 4 — 300 K. 
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Figure 2.13 Thermal conductivity of Aluminum Alloy 1100, computed from tables 



available on the NIST website (NIST 2008). These values have been measured in the 
temperature range T = 4 — 300 K. 
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Figure 2.14 Thermal conductivity of G-10, computed from tables available on the 



NIST website (NIST 2008). These values have been measured in the temperature 
range T = 4 - 300 K. 

of the cryostat outer case (which is, of course, ~ 300 K) and the radiation shield (see 



Fig. 2.10). The radiation shield equilibrates at ~ 80 K, even though the first stage of 
the refrigerator is at ~ 55 K, due to the poor conductivity of the aluminum of which 
the shield is made. This blanket reduces thermal loading due to emission, as each 
layer is highly reflective and will equilibrate to some temperature between 90-300 K. 
Thermal emission can be computed using the Stefan-Boltzmann Law (P = bcxT'^A), 
where b is the emissivity coefficient of the material {b ~ 0.10 for aluminum), a is 
the Stefan-Boltzmann constant, T is temperature (in K), and A is the area of the 
absorbing/emitting surface. The strong temperature dependence encapsulated in the 
Stefan-Boltzmann Law implies that the innermost blanketing layer, which equilibrates 
at ~ 80 K, radiates ~ 200 times less than the 300 K cryostat inner wall. 



In Fig. 2.7, the item labeled #19 is called a "getter." The getter contains activated 
charcoal, which acts as a trap for contaminants in the cryostat. The charcoal, which 
is cold, provides a large surface area onto which contaminants are adsorbed. Note 
that the final receiver configuration had the getter mounted to the cold plate. 



44 



Zotefoam was chosen as the material out of which to build the receiver windows, 
as it also can hold vacuum (in addition to the properties discussed in the previous 
section). However, zotefoam is very elastic, and a window made of zotefoam alone 
bows under vacuum enough to be compressed against the corrugated feed horns. 
Since these feed horns couple directly via waveguide to the HEMT amplifiers, and 
are mounted on the electronics plate, they are ideally kept at the same temperature 
as the amplifiers. Even given the poor conductivity of zotefoam, this strong contact 
between it and the feed horns is enough to dominate the thermal load on the cold 
stage. 

The bowing problem was solved by introducing a layer of mylar as a backing for 
the zotefoam. Thin, microwave transparent mylar, like that in the receiver windows, 
is strong enough to support the zotefoam, but is not robust to sharp impacts. The 
zotefoam layer acts as a mechanical protection for the mylar, and makes the vacuum 
seal. 

The windows are held by a custom-designed piece machined from a solid alu- 



minum block (see Fig. 2.15). Two holes, each about 1.2 times larger than the outer 
diameter of the feed horns, allow the signal to enter. The bed of the holder provides 
a large, fiat surface which was sandblasted to further increase the surface area. This 
surface was coated with a thin layer of epoxy (specifically, a clear Emerson & Cuming 
Eccobond® 24, two-part epoxy with high viscosity, which avoids getting epoxy in un- 
wanted places), and a single layer of mylar is joined here, carefully removing bubbles 



by hand from the interface (see Fig. 2.16 for an illustration of this interface). After 
the first layer sets, the zotefoam is joined here by applying another layer of epoxy 
directly above the first layer. This is done to avoid having epoxy in the optical path. 
After this layer of epoxy sets, one last layer of epoxy is used to fill in the ~ I/IO*'^ 
inch gap around the zotefoam's edge, and a protective plate with a large elliptical 
gap for the signal is mounted on the front. The final step in this process is to add a 
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epoxy surface 




Figure 2.16 Cross-section of the window holder, depicting where the first layer of 
epoxy is applied. 



46 



pinhole in the mylar (on the vacuum side), allowing any air/moisture to escape from 
between the mylar/zotefoam interface when vacuum is drawn. The final product is a 
zotefoam window that bows gently into the mylar without making contact with any 
of the components inside the cryostat. 



2.4 Back-end Electronics 



2.4.1 Overview of the back-end electronics 



All support electronics for the receiver are located in a large, thermally-regulated. 



insulated box, called the "Ebox" (see Fig. 2.17). The Ebox is placed within the 
receiver enclosure, just outside the cryostat. The primary components in this box 



(shown schematically in Fig. 2.18) are: 



An IF switch, used to select the receiver output appropriate for 30 or 90-GHz 



observations. The IF switch is mounted in location ^20 in Fig. 2.17 



A pre-amplifier module (PAM in Fig. 2.18 #1 in Fig. 2.17), which is primarily 
an IF amplifier, is common to both receivers as it follows the IF switch. The 
PAM also includes a set of precisely-calibrated attenuators (in 0.5 dB steps, up 
to a total attenuation of -30 dB), both on the input and output of the module. 
These attenuators are useful in avoiding compression effects when the power 
becomes high enough to saturate electronics further down the signal chain. 
Since the hot load calibrator (a thermal blackbody at ambient temperature) 
increases the receiver output power by a factor ~ (300 K)/Tsys ~ 4 — 6, where 
Tgys ~ 40-50 K, these attenuators are also inserted during calibrations to prevent 
compression of the signal (for more information, see Q. 



A power supply for the receiver (labeled "rxmod" in Fig. 2.18 #2 in Fig. 2.17). 
This powers both the 30 and 90-GHz receivers. 
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Figure 2.17 Photo of Ebox (lid removed), set up for 30-GHz observations. The Ebox 
houses the modules used to control and process signals from the receiver. See text in 
^2.4.1 for a description of the electronics. Note that the bias-tuned Gunn (BTG), the 
module to control the bias-tuned Gunn (the "BTG Mod" in Fig. 2.18), and the IF 
switch had not been installed at the time of this photo, as they were only required for 
90 GHz observations. They were later installed in locations 12 and 20, respectively. 
The fiber bundle and power cables enter the Ebox at location 16. The refrigerator 
head (13) for the receiver (located behind the Ebox), the receiver enclosure's chiller 
line (14), the chiller heat exchanger (15), and the TEC cooling fans (17) can all be 
seen above and outside the Ebox. The internal Ebox air circulation fans at label 18. 
A spring that assists in lifting the Ebox is labeled 19. The walls of the Ebox are 
lined with an open-cell PVC foam for insulation, which was painted white to prevent 
degradation and flaking due to weathering. 
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Figure 2.18 Layout of Modules in the Ebox. Connectors are shown in green. 

• An high-frequency oscillator, based on a Yttrium Indium Gallium (YIG in 



Fig. 2.18, #10 in Fig. 2.17) diode in a cavity, that operates at 8.972 GHz. 



The module that monitors the YIG's status is labeled #3. 



The varactor, a "variable capacitance" component (#11 in Fig. 2.17) that is 
used to produce harmonics of the YIG frequency as well as lock the phase of 
the LO. A phase-lock loop (PLL), controlled by electronics in the correlator 
trailer, locks the phases of all antennae to a single reference. The module that 
controls the varactor is #3. 

A bias-tuned Gunn (BTG) oscillator, which is the 90-GHz system's equivalent 
of the YIG. This uses a Gunn diode forced into oscillation within an adjustable 



waveguide cavity. This was later installed in location #12 in Fig. 2.17, while 
the module to control it was installed between #3 and #4 ("BTG Mod" in 



Fig. 2.18). 



A dielectric resonance oscillator (DRO), which supplies the fixed 17.5 GHz sec- 
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ondary LO used by the 90-GHz system's secondary mixer. This is mounted on 



the same plate as the BTG, so it was also installed in location #12 in Fig. 2.17 



Various CANBUS (a protocol like TCP/IP that allows daisy-chaining) con- 
trolled electronics, used to control the various oscillators and reference frequen- 
cies of the PLL and lobe rotator (which accounts for the part of fringe tracking 
due to movement of the source across the sky), regulate the temperature, etc. 



Among these is the interface module ("Int Mod", #5 in Fig. 2.17), which mon- 
itors the phase reference and lobe rotator. These references are received by the 



10 and 50 MHz receivers (#7 and #8 in Fig. 2.17). 



The local oscillator terminator ("LO Term" in Fig. 2.18, 7^6 in Fig. 2.17), which 



locks the YIG to the ninth harmonic of a central, synthesized reference frequency 
with the 10 MHz signal mentioned above. This phase lock loop is separate from, 
but analogous to, the one which locks the varactor to the YIG using the 50 MHz 
reference. 



An optical IF band transmitter (called the "OTX"; labeled "IF Tx" in Fig. 2.18 



#9 in Fig. 2.17), which takes as its input the output of the PAM. This sends the 



signal over fiber optics (^^16 in Fig. 2.17) to the downconverter and correlator 
(which are discussed in ^2.5.1 & 2.5.2, respectively). 



2.4.2 Thermal considerations for the electronics box 

Thermal stability is the primary motivation for the implementation of this insulated 
electronics box. The thermal sensitivity of the amplifiers, local oscillators, and fiber 
optics (whose lengths change with temperature) leads to phase variations during peri- 
ods of poor thermal regulation. While these effects, discovered during commissioning 
operations, varied from antenna to antenna, one extreme example exhibited ±45° 
phase oscillations as the temperature varied sinusoidally about its nominal regulation 
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temperature by ±2.5°C (which could happen even with the final thermal regulation 
in place). Since these oscillations had periods on the order of 15-20 minutes, the same 
frequency as that of our calibration observations, the sampling of these phase changes 
is not frequent enough to remove the effects through data calibration. This problem 
was eventually solved with better regulation algorithms and the occasional manual 
intervention by the observer. 

Thermal regulation of the electronics box, utilizing thermo-electric coolers (TECs, 
or Peltier junctions), is set up so that the temperature at a single point on the 
electronics mounting plate is kept within ±0.125°C of the 32°C (28°C) set point in the 
summer (winter). The thermal differentials within this box are reduced by fans, and 
by the large thermal mass and high thermal conductivitj- of the 1 / 4-inch aluminum 
plate to which all the electronics are mounted. A small amount of silicone grease 
(heatsink compound) is used to ensure decent conductivity between the electronics 
modules and the mounting plate. 



A conceptual overview of the TEC setup for the electronics box is show in Fig. 2.19 



Ambient air is forced, by a large impeller fan (#17 in Fig. 2.17), through the fins of the 
heatsink attached to the hot side of each TEC, thereby carrying away the heat. Good 
thermal contact between the heatsink and the TEC is ensured by the pressure plate, 
which is held down by five springy Belleville washers, each providing 20-30 pounds 
of force. This setup also forces the cold side of the TECs onto a copper plate that 
makes thermal contact with the electronics mounting plate, labeled "worksurface" in 



Fig 2.19 Since the TECs are ceramic, their porous surfaces are coated with a layer of 
heatsink compound, which fills the pores and provides more surface area for thermal 
contact. Open-cell foam insulation isolates the temperature-regulated copper plate 
and worksurface from the warmer outside air. 
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Note, as in Figs. 



2.13 


and 


2.12 


copper, 


3ut is 



that - at room temperature - aluminum is only a factor of '-^ 2.5 



consideration, as the Ebox has to be lifted from the receiver enclosure to access the receivers. 
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TEC Blower Placement - Top view 

Air flows from the blower to the fins through an air duct. 




Air Duct 



Air Intake 



TEC Blower Placement - Side View 
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Figure 2.19 TEC Blower Conceptual Diagram, courtesy Marshall Joy and Georgia 
Richardson. Ambient air from the receiver enclosure is forced through the heatsink 



fins by the blower (see Fig. 2.17, #17). These copper heatsinks are thermally coupled 
to the hot sides of the TECs (which are > 20°C hotter than the worksurface side). 
The TECs diffuse more heat away from the worksurface than toward it, thus removing 
heat from the Ebox. 

The TECs can be operated at a maximum power of Pmax = 300 W (10 A at 30 V, 
DC), removing Pmax.remove = ^-Pmax for TEC efficiency ^. With cooling efficiencies 
of ^ ~ 30%, they can keep a thermal load of up to 90 W regulated at a desired set 
temperature. The electronics within the Ebox draw ~3 A at 24 V (DC), and therefore 
dissipate ~ 72 W of heat. There is a remaining thermal load that is due to both the 
finite thermal conductivity of the G-10 supports and foam insulation of the Ebox, 
and to the imperfect thermal isolation provided by the Ebox. 
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The Ebox insulation is a rigid PVC foam, chosen because it is the best insulating 
foam offered by McMaster-Carip^ The conductivity of the foam is 90 mW ■ in/ft^/K 
(0.17 btu ■ in/ftV°F). With ~ 11 ft^ of 0.375 in foam insulation, 2.6 W of heat flow 
through this insulation for each degree of differential between the outside and inside of 



the box (see Eq. 2.4), In practice the electronics box cannot maintain a temperature 
more than ~ 1°C less than the ambient temperature of the receiver enclosure, due to 
leaks that allow some air exchange between the inside and outside of the box. 

2.5 The Correlator Trailer 

The signal from each telescope is sent over fiber optic cables, located in conduits 



underground (as illustrated in Fig. 2.1), to an air-conditioned trailer - referred to as 
the "correlator trailer" - containing the instrumentation necessary to downconvert 
([ |2.5.1 ) and correlate {[2.5.2) the signals. The correlator trailer also houses the 



computer that controls the array and collects all data taken during an observation. 
2.5.1 Downconvert er 

The 1-9 GHz IF is processed by an array of mixers, amplifiers, and bandpass filters, 
known collectively as the "downconverter," which splits the IF into sixteen 500 MHz 



bands running from 0.5-1 GHz (see Fig. 2.20). The reason for using this range is that 



electronics that can handle 0.5 GHz and 1 GHz are similar, while those that could 
handle a signal from 0-0.5 GHz; for instance, DC signals cannot be capacitively 
coupled. 

Here I describe the full signal downconversion chain, starting from the sky fre- 
quencies 26.938-34.938 and 90.78-98.78 GHz, and ending with the sixteen digitized 



500 MHz bands that constitute the input to the correlator. Figure 2.20 illustrates 
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Figure 2.20 Block Diagram of the SZA Band Downconversion, showing how the sky 
observation frequencies 26.938-34.938 and 90.78-98.78 GHz are spht into the sixteen 
digitized bands, each 500 MHz in bandwidth, that are the input to the correlator. 
See description in the text ( ^2.5.1 ). Bands shown in red are the side bands that are 
not used. 



54 



Band 


IF Band Range 


30-GHz Sky Frequencies 


90-GHz Sky Frequencies 


(#) 


(GHz) 


(GHz) 


(GHz) 





1.0-1.5 


34.688 


91.030 


1 


1.5-2.0 


34.188 


91.530 


2 


2.0-2.5 


33.688 


92.030 


3 


2.5-3.0 


33.188 


92.530 


4 


3.0-3.5 


32.688 


93.030 


5 


3.5-4.0 


32.188 


93.530 


6 


4.0-4.5 


31.688 


94.030 


7 


4.5-5.0 


31.188 


94.530 


8 


5.0-5.5 


30.688 


95.030 


9 


5.5-6.0 


30.188 


95.530 


10 


6.0-6.5 


29.688 


96.030 


11 


6.5-7.0 


29.188 


96.530 


12 


7.0-7.5 


28.688 


97.030 


13 


7.5-8.0 


28.188 


97.530 


14 


8.0-8.5 


27.688 


98.030 


15 


8.5-9.0 


27.188 


98.530 



Table 2.1 The 1-9 GHz IF band is separated into 16 bands, as outlined in Fig. 2.20 



These are the central sky frequencies of both the 30 and 90-GHz SZA systems that 



correspond to the downconverted IF bands, discussed in §2.5. 1[ Note that higher 
band number corresponds to a lower sky frequency for the 30-GHz system because it 
uses the lower side band. 

this signal chain. 

The tertiary mirror of each telescope selects either the 30 or the 90-GHz receiver. 
In the 30-GHz receiver system, the lower side band (LSB) mixing product of the sky 
with a 35.938 GHz local oscillator (LO) places sky frequencies 26.938-34.938 GHz in 
the 1-9 GHz IF band. In the 90-GHz receiver, a 72.28 GHz LO (in the current tuning 
scheme) mixes with sky frequencies 90.78-98.78 GHz, producing a 18.5-26.5 GHz 
upper side band (USB) product; a second LO at 17.5 GHz mixes with this to place 



the 90-GHz receiver output (also USB) in the 1-9 GHz IF band (see Fig. |2.21| for 
more details on the downconversion in the 90-GHz receiver). Note that the 30-GHz 
receivers have a full bandwidth of 10 GHz, while the 90-GHz receivers have a much 
broader bandwidth that allows observations from 85-115 GHz. For both systems, the 
back-end electronics determine what range of sky frequencies are observed. 
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Figure 2.21 Detailed Diagram of the 90-GHz Receiver Downconversion, showing how 



the sky observation frequencies are brought down to the 1-9 GHz IF band (see Fig. 2.8 



for details on the receiver RF components). This figure accompanies the broader 
downconversion scheme illustrated in Fig. |2.20 In the current tuning scheme, sky 



frequencies 90.78-98.78 GHz are mixed ((8>) with a 72.28 GHz LO (the bias-tuned 
Gunn), producing an 18.5-26.5 GHz USB product. Sky frequencies < 85 GHz are 
blocked by the high pass filter, so no LSB mixing product is produced (e.g. the band 
shown in red). The bandpass filter passes the 18.5-26.5 GHz USB product to a second 
mixer. The second LO, at 17.5 GHz, mixes with this 18.5-26.5 GHz USB product to 
place the 90-GHz receiver output in the 1-9 GHz IF band. This 1-9 GHz IF signal 
is the input to the back-end electronics common to both the 30 and 90-GHz systems. 
See 92.5. 11 for more details. 
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Back-end electronics in the electronics box select the proper receiver output and 



send the signal over fibers to the correlator trailer (see Fig. 2.1). The IF is then split 
into 1-5 and 5-9 GHz bands by the block downconverter. The lower half of the IF, 
from 1-5 GHz, is passed directly to the second stage downconverter. The 5-9 GHz 
band is mixed with a 10 GHz LO, also in the block downconverter, and the LSB 
mixing product (from 1-5 GHz) continues to the second stage downconverter. I use 
the ® symbol to represent the mixer, which is located within the block downconverter, 
to clearly denote how the 1-9 GHz IF band is first split into 1-5 and 5-9 GHz. 

The second stage downconverters mix and filter each 1-5 GHz band into sixteen 



separate 0.5-1.0 GHz bands, listed in Table 2.1 Each 0.5-1.0 GHz band is then 



digitized at a sampling frequency fs = 1 GHz. This sampling operation digitally 



aliases each 0.5-1.0 GHz band to a 0-500 MHz baseband, digitized signal This 



sampled signal is the input to the SZA correlator, discussed in the next section. 



2.5.2 The SZA Correlator 



In the discussion of interferometry presented in ^1.3, details about the process of 
interfering the signals from a pair of antennae, and the instrument that performs this 
cross- correlation, were for simplicity ignored. I now discuss the correlator, which is 



the single instrument most crucial for performing interferometry. Fig. 2.22 shows a 
photo of the SZA correlator. 

The cross-correlation rxy{m) for digital time delay m of two discrete-time signals 
x{n) and y{m + n) (delayed by m samples with respect to x{n)), such as the digitized 
signals from a pair of antennae, is defined 



oo 



rxyifn) = a;[n] ?/[m -|- n]. (2.5) 



^^Recall that the Nyquist theorem states fg > 2/bw, where fs is the samphng frequency and /bw is 
the bandwidth of the signaL A samphng frequency of 1 GHz hmits the sampled signals bandwidth 
to 0.5 GHz. 
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Here, is a sample taken at time t = nr for sampling period r = 1//^, and m a whole 
number lag - a discrete time delay - applied to one antenna's signal. Since the SZA 
bands are sampled at = 1 GHz, r = 1 ns. In a real instrument, we cannot sum 



over the infinite number of samples indicated in Eq. 2.5 For the SZA, correlation 



is performed for 32 lags of each signal, after having been properly delayed by fringe 



tracking and downconverted to 0.5-1 GHz bands (as discussed in § 2.5.1). The overall 
effect of the correlator is to multiply the signals from two antennae, sliding one signal 
past the otherj^ Since the signals x{n) and y{n) scale as voltage, the correlator 
output scales as voltage-squared; cross-correlation is therefore proportional to the 
power received by an antenna. 

We are accustomed to dealing with visibilities as a function of frequency. Visi- 
bilities, when calibrated against an astronomical source of known fiux, give the fiux 



at each data point in ^jf-space, for each frequency in a spectrum (see {!.?>]. The 
spectrum is simply obtained by Fourier transforming the correlated data. 

The SZA correlator is an XF or lag-correlator, where X stands for the cross- 
correlation, and F stands for the Fast Fourier Transform (FFT). The SZA correlator 
performs cross-correlation in the time domain, but note that FX correlators are also 
common, since correlation has a simple Fourier relation 

r^y{m)=F{X[kYY[k]}, (2.6) 

where X[k\ and Y[k\ are the Fourier transforms of x{n) and y{n + m) for discrete 
frequencies /c, J-" is the Fourier transform operation, and * denotes the complex con- 
jugate, 



18 



The FFT is a computationally efficient implementation of the Discrete Fourier 



^^For comparison, note that convolution, a similar operation to correlation, is the flipping and 
sliding of one signal past the other. 

^^Note for comparison that convolution in the time domain is simply multiplication in the fre- 
quency domain (i.e. neither signal is complex-conjugated, unlike in correlation). 
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Figure 2.22 Photo of the SZA Correlator and Downconverter, courtesy David 
Hawkins. 

Transform (DFT), which in turn is a truncation of the Discrete Time Fourier Trans- 
form (DTFT). The DTFT is mathematically complete (e.g. it satisfies Parseval's 



Theorem; see Oppenheim et al. 1999 Bracewell 2000 for further information), and 



becomes the familiar (continuous time) Fourier Transform in the limit of infinitesi- 
mally small steps in time. 
The DTFT is: 

oo 



where 



x(n) = (l/27r) 



X(uj) e'""du. 



(2.8) 
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The DFT is a special case of the DTFT, where the number of samples is limited to 
a finite number N. This has the effect of discreting frequency at u = 2TTk/N, where 
k = N{f / fs) is the discretized wavenumber. In the DTFT, X{uj) wraps outside the 
range uo = [—tt,tt]; in the DFT, X{uj) is sampled at only the discrete values of k, and 
X{k) wraps outside the range k = [0,N — 1]. The DFT is: 

JV-1 

X{k) = Y,^in) e-''^'"^" (2.9) 

n = 

where 

JV-1 

= (1/iV) 6^^"'="/^ (2.10) 

fc=0 

The substitutions that take us from the DTFT to the DFT have the effect of limiting 
the frequency resolution to N/2 elements (in the Nyquist bandwidth of /s/2; there 
are N discrete fc's in the range / = [—/s/2, /s/2]). Notice that the wrapping in the 



DTFT and DFT (Eqs. 2.8 & 2.10) is the same as wrapping every fg in frequency, or 



every harmonic of the sampling frequency; this aliasing is simply that due to Nyquist 



sampling at fs, mentioned in ^2.5.1, Only the range k = [0, A^/2 — 1], corresponding to 
frequencies u = [0,7r], contains unique spectral information. Since the time-domain 
signal is real, its transform is Hermitian, meaning those frequencies above fs/2 are 



the complex conjugate of the positive frequencies. Table |2.2| lists the 17 channels 
produced by FFT-ing the 32 samples. 

In the time domain, the limited number of correlator lags m used in the DFT 
is equivalent to rectangular windowing - the one-dimensional equivalent to the top 



hat function - of the correlated signal rxyim) (Eq. 2.5). Outside of the range n = 
[0, — 1] (in discrete time), the signal is zero, and this has important implications 
when interpreting the DFT as a power spectrum (which are discussed in ^3.1.3 ). 



The mathematics of the cross-correlation and Fast Fourier Transform (FFT) are 
programmed into field-programmable gate arrays (FPGAs), which can be faster than 
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Channel 


Frequency 


Channel 


Frequency 


(#) 


(MHz) 


(#) 


(MHz) 





0.00 


cont'd 


cont'd 


1 


31.25 


9 


281.25 


2 


62.50 


10 


312.50 


3 


93.75 


11 


343.75 


4 


125.00 


12 


375.00 


5 


156.25 


13 


406.25 


6 


187.50 


14 


437.50 


7 


218.75 


15 


468.75 


8 


250.00 


16 


500.00 



Table 2.2 Each 500 MHz band has 17 channels, numbered 0-16. Channels and 
16 are both attenuated by the 500 MHz bandpass filters used in the downconverter 



2.5.1, and are not used for data. Note that these are discrete frequencies that arise 



from the DFT of a discrete-time signal (Eqs. 2.9), rather than continuous bands 



using general-purpose processors to perform digital correlation. The raw visibilities 
produced by this process are averaged for 0.5 second intervals and stored to disk as 



raw visibilities in correlator units The reader is encouraged to see Hawkins et al. 



(2004) for further details, where the SZA correlator is described in much greater 



detail. 



2.5.3 Thermal Considerations for the correlator trailer 

The correlator trailer draws ~ 1 MW of power, and thus requires an industrial- 
strength air conditioner (A/C) with two independent compressor units (though much 
of this power is drawn by the A/C itself, and its heat is dissipated outside). At 
the beginning of SZA commissioning, we used a state-of-the-art air conditioner with 
a commercially-available governor that can keep the temperature within ±1°C of a 
desired set point. It was soon discovered, due to the temperature sensitivity of the 
digitizers, that the 20 minute air conditioner cycling of ±0.25°C led to ±3° phase 



variations in the data (see Fig. 2.23), and a more robust solution had to be found. 



^The correlator electronics are set up such that the auto-correlation of each antenna's signal 



applying Eq. 2.5 with x{n) ~ yin) to find how a signal correlates with itself - is 1. Deviations from 
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Figure 2.23 Plot of phases on a single baseline (Baseline 0-1) and the corresponding 
downconverter temperatures during the same period. The blue lines are the tempera- 
tures measured in the downconverters for Antennae (dashed line) and 1 (solid line). 
The black line is the phase of the raw data on baseline 0-1, taken while staring at 
a strong point source, which ideally has flat, constant phase (see ^4.3.2 for details). 
Note that thermal variations in adjacent digitizers, due to cycling of the A/C, of 0.5°C 
(peak to peak) produced 6° phase variations. Modifications to the A/C removed these 
thermal oscillations on short time scales (see Fig 2.24). 



The resolution to this problem was largely worked out by the SZA engineer, Ben 
Reddall. He removed the governor and set the A/C to be continuously on. He then 
modified the unit so that louvers could push the air through the heat exchanger of 
one compressor at a time, so that the other compressor could be shut off when not 
necessary. This strong, continuous air flow in a closed cycle ensures good mixing 



this are used later, in calibration (Sffl, to diagnose problems in the correlator. 
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Temperature/Phase for the Downconverter, after upgrades 
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Figure 2.24 Plot of phases on a single baseline (Baseline 0-1) and the corresponding 
downconverter temperatures during the same period. The blue lines are the tempera- 
tures measured in the downconverters for Antennae (dashed line) and 1 (solid line). 
The black line is the phase of the raw, uncalibrated data on baseline 0-1, taken while 
staring at a strong point source, with the linear slope and mean value removed. Note 



that this is on a much shorted time scale than that of Fig. 2.23 but is long enough to 
have shown 3 full cycles of the A/C if they had persisted. Any residual phase error 
is both negligble and is removed by calibration (see §4^ 



and stability within the trailer; it also prevents the active compressor from icing up, 
as it essentially doubles the airflow through it. Various layers of insulation further 
contributed to the thermal stability of the trailer]^ The final steps implemented by 
Ben Reddall were the addition of water bottles (water has a high thermal capacity) 
to the air supply duct, and of three 1 kW space heaters used during the winter. 
The heaters help prevent the A/C, which must be run continuously, from freezing on 

^"Note that the original OVRO correlator was housed in a basement underground, and the newer 
CARMA correlator ~ for which the SZA correlator served as prototype ~ is in a large, heavily- 
insulated room with a concrete foundation, rather than in a wooden trailer located above pavement. 
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cold winter nights. Figure 2.24 demonstrates the continuing success of the correlator 
trailer thermal solution, three years after it was implemented. 

The temperature in the correlator trailer typically varies by 20° C over the day, but 
does so in a way that these slow, diurnal phase changes are common mode to the fibers 
and electronics. The slowly-changing phase drifts produced by these temperature 
changes can be tracked using our standard astronomical calibration (see Q and 
removed from the data. 



Chapter 3 
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Sources of Contamination 



This chapter addresses two sources of contamination found during the SZA com- 
missioning phase: radio interference from external sources and antenna cross-talk. 
While I have attempted to make this palatable for those who are not engineers or 
radio astronomers, the information is best suited to an audience interested in building 
and debugging a radio interferometer. 



3.1 Radio Interference 

Here we consider the problem of radio frequency interference (RFI), which we found in 
our early observations to be contaminating frequency channels or even entire bands. 
A "birdie," the colloquial name for RFI, is in general some contaminating, non- 
astronomical signal, which can increase the scatter in the radio visibilities at the 
birdie frequency, corrupting data in a given channel or even a whole band (discussed 



in ^3.1.3). The strongest birdies are generally from sources local to the antenna 
electronics, and could be the mixing product of the local oscillator (LO), the LO 
reference, and/or the YIG. Birdies can also come from external transmissions, gen- 
erally ground-based sources (e.g. military communications or satellite uplinks) or 
geostationary satellites. 

In this section, we examine ways of lessening the impact of birdies, since they 
could potentially contaminate maps made from observations with the SZA. For sim- 
plicity, we consider a birdie to be an ideal sinusoidal signal. We calculate how much 
fringe tracking helps to reduce birdies, since they are stationary with respect to the 
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antenna^!] and lose coherence in the cross-correlated data. 

Early observations with the SZA revealed the presence of a strong birdie, just 
outside our IF band, at the receiver output frequency 9.022 GHz. This birdie arose 
in the output of the receivers as the mixing product of the YIG 3'''^ harmonic (3 x 
8.972 GHz = 26.916 GHz) and the 35.938 GHz LO (35.938 GHz - 26.916 GHz = 
9.022 GHz). The 30-GHz receivers detected both the YIG 3"^^^ harmonic emitted 



within their own Eboxes and those from adjacent antennae (^2.3.3 for information 
on the YIG oscillator). This birdie showed up in Band 15, which is the lowest sky 
frequency band (26.938-27.438 GHz), despite the fact that this band runs from 8.5- 
9 GHz, and the birdie was 9.022 GHz. 

We demonstrate in this section that an additional 23 dB of attenuation, at the 
birdie's source, was required to bring the birdie to less than 10% of the rms noise 
level via fringe tracking. We then proceed to consider how a birdie can corrupt all the 
channels in a band, even if the birdie is a pure tone, singular in frequency. Finally, we 
demonstrate that a clever retuning of local oscillators can avoid band corruption by 
insuring all the birdie's power appears in one channel, which can readily be flagged 
and excised from the data. 

3.1.1 The Effect of Fringe Tracking 

Since birdies come from sources not associated with the field being observed, they 
can be treated as an untracked source. Since the SZA birdie was local to the antenna, 
it was effectivley moving across the sky in the opposite direction as the astronomical 
source being observed. The instantaneous correlations of the birdie signal affecting 



^As mentioned in the first paragraph of { 



3.1 



birdie sources are generally local to the antennae, 
from ground-based transmissions, or from geostationary satellites; all of these sources are stationary 
with respect fringe tracking. Other possible birdie sources include satellites in low earth orbits, 
which would move through the fringes on a given baseline much more rapidly than the astronomical 
sources we track (e.g. 45 minute orbits versus the sidereal day that a celestial object takes). We 
found no indication of contamination from any satellites with an orbital period on the order of a 
day, which - if they exist - could produce coherent, cross-correlated contamination. 
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a pair of antennae do not sum up coherently over time (see ^1.3), and the birdie 
is therefore attenuated by this loss of coherence. The factor by which the birdie is 
attenuated, F^, can be calculated as a sine function of the fringe frequency Pf (the 
rate at which an untracked source moves across the sky) and the integration time Tq. 
The attenuation factor is 

^^^sjni^^ (3.1) 

where i^/ is the fringe frequency 

Uf = uijJ^cos{S), (3.2) 

We is the earth's angular velocity, {u^ = 7.29115 x 10^^ rad s^"*^), 6 is the source's 
declination, and u is the Fourier conjugate (in -space, the standard notation used 
by radio astronomers) of the baseline along the east-west direction, normalized by 
the wavelength of observation: 

u = X^sinH + Y^cosH. (3.3) 

Here H is the hour angle of the source, where H = 0^ corresponds to the meridian, 
and each hour increment is 15°. Xx and Yx are respectively the x (north-south) and 
y (east-west) components of the baselines divided by the wavelength at which the 
observation is made. A longer baseline increases the fringe frequency, since it yields 
a tighter, more highly-resolved fringe pattern on the sky. A source at 5 = 0° (in 
the equatorial plane) moves most quickly across the sky, so the fringe frequency Uf is 
highest there. A point source at the pole, on the other hand, would not move through 
fringes at all. 



The l/iTUfTa envelope of the birdie attenuation factor (denominator of Eq. 3.1) 



results from the time averaging of the measured visibities, the smallest binning of 
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which is a 1/2 second frame (see ^2.5.2). We are interested in how much the birdie is 



attenuated during an observation, and how this compares to the Gaussian noise floor 



of the same observation. The noise floor decreases according to (see, e.g. Rohlfs & 



Wilson 1996 Thompson et al. 2001): 



AT. 



M 



sys 



a/Az/ X Ta 



(3.4) 



In the above equation, M is a factor - less than or equal to unity - that accounts for 
the quantization noise in the correlator system]^ is the integration time, Az/ is the 
bandwidth, and T4 is the antenna temperature. 



3.1.2 Can Fringe Tracking Beat Down the Birdie? 

We consider here if the SZA 9.022 GHz (IF band) birdie that plagued Band 15 of 



the SZA 30-GIIz system (see Table 2.1) in the first few months of operation could 
be sufficiently attenuated by fringe tracking. We assume our source is at declination 
5 = 0°, which is our best case scenario; this yields the highest fringe frequency, since 
the source moves most rapidly across the sky. If our source were at the north celestial 
pole, on the other hand, fringe tracking would not attenuate the birdie at all. Clearly, 
we cannot rely in all cases on fringe tracking to eliminate birdies. We will also for 
simplicity assume our source is close to the meridian (i.e. near transit). With if ~ 0, 



the Xx component is negligible (see Eq. 3.3), simplifying the calculation 



Let us set Y\ = 1000 (see Eq. 3.3). For the SZA, this would be among the longest 



of the inner array baselines, with a 10 meter, east-west baseline component while 
observing at 30 GHz. If we chose a purely north-south baseline (Fa = 0), a source 
near transit would give m = 0, and fringe tracking would not attenuate the birdie at all 



see Eq. 3.2). Again, this indicates that fringe tracking is not a panacea for reducing 



^We do not need to know M here. However, M 0.88 for the SZA correlator (Hawkins et al 



2004). 
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birdies, but we proceed to compute how much it helps in our best case scenario for 
the inner array. 



Calculating the fringe frequency (Eq. 3.1 ) for these conditions, Uf ^ 7.3 x 10 ^ Hz. 



We can see that Fi, hits its first null at r^, = 1/z// ~ 14 sec, and we might be tempted 
to choose this as our integration time so that no stationary sources could affect us. 
Since this frequency is different for each baseline, it is clear that this scheme will not 
work. 

We can proceed in calculating the birdie's attenuation without needing to know 



M and Tgyg (in Eq. 3.4) by measuring the ratio of the birdie signal power with fringe- 
tracking off (which is therefore unattenuated by integration) to the noise level in a 
known time interval. Given the functional behavior of the birdie's attenuation (Eq. 



3.1) we only measure the birdie's unattenuated amplitude Abfl and the noise's rms 



amplitude An{Tm) in a given time interval r^, with fringe tracking off. It should 



be noted that the noise integrates down as 1/a/t (as in Eq. 3.4) without regard to 
fringe tracking, so it is reasonable to compare the noise levels for given integration 
times, regardless of whether fringe tracking is on. The birdie and fringe amplitudes 
as functions of time are given by: 

AM = (3.5) 

and 

An{'Ta)^/Ta = (^m) V^Tm (3-6) 

SO that the ratio of the birdie to the noise level for integration time is: 

With fringe tracking off, the level of the 9.022 GHz birdie was measured to be 
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~ 250 times the noise in a 20 second, uncontaminated integration. Since it's more 
convenient to express integration times in seconds, we solve for Ah^ normalized to 
the noise in a 1 s integration: 



250 

A,o = 250 A„(20 s) = ^An(l s). (3i 



Plugging the ratio Abfi/An{l s) into Eq. 3.7, we find the ratio of the birdie to the 
noise level after time (in seconds): 

A(Ta) 250 

AM^TTi^'^^'^^' (3.9) 



Taking the local maxima of Fi, (defined in Eq. 3.1) as our worst case scenario for the 
birdie's amplitude, we ignore the sinusoidal component of the sine function and keep 
the l/^TTUfTa) envelope. This yields: 

Mra) 250 1 17.8 



^n(Ta) \/20 ^^h'f^/T~a V^^a 

Now, we can solve for the time when the envelope for the birdie amplitude falls 
below the noise fioor. For our fringe rate Vj^ we find that it takes about 16.5 hours 
for this birdie to integrate down to the same level as the noise. This of course violates 
the assumption that the source is near transit. Furthermore, this time scale is much 
longer than any amount of time we would sensibly bin as a single data point. 

We would like the birdie to be attenuated to < 10% of the noise level in each 
binned data point, so that maps are relatively unaffected by the birdie (with a ~ 
5-cr cluster detection, the birdie would affect the SZE flux on the 2% level, which 
is less than our calibration uncertainty). To get this birdie down to the 10% level 
in a 5 minute integration, the maximum time which we might reasonably bin into 
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a single data point, we need about 22 dB of attenuation on each antenna]^ When 
observing sources at higher dechnations, additional attenuation would be required to 
compensate for the lower fringe frequency uj. A source at 6 = 37° takes 25% longer 
to beat the birdie down to the noise level of a 5 = 0°, since its fringe frequency is 
80% of the 6 = 0° value. Conversely, an extra dB of attenuation would be required to 
bring the birdie down to the same level in the same time. Therefore, we would need 
23 dB of attenuation of the birdies at their sources in order for this measured birdie 
not to affect a typical SZA observation. 



In Eq. 3.10, we ignored the fact that Fi, varies sinusoidally (as in Eq. 3.1). In 
a typical SZA observation, many integrations are taken, and are often binned into 
intervals ranging from 20-300 s. For each binned data point, the birdie level would 
vary randomly (since Uf changes throughout the observation), with both positive and 
negative variations that are several times larger than the noise (e.g. 250 times larger 
for a 20 s observation, as measured before attenuation). Thus, as long as the fringe 
rate Uf ^ a. birdie would increase the scatter in the data taken at that frequency, 
rather than integrating to a coherent signal. 

3.1.3 Spectral Leakage 



As noted in ^2.5.2, the limited number of samples used in correlation has the effect 
of rectangular windowing of the cross-correlated data. In this section we examine 
how a sinusoidal signal - such as an idealized birdie - can contaminate an entire 
band of data, rather than only appearing at its frequency. This phenomenon is called 
"spectral leakage." 

Consider a birdie at frequency /o = k^fs/N, where the signal is sampled N discrete 
times at a sampling frequency fg, and fco is the birdie's wavenumber. Let x{n) be 



As noted in S 2.5.2 



cross-correlation scales as power. Therefore, each antenna-based birdie's 



power has to be attenuated by the same amount as the total desired birdie attenuation, 
rectangular window is the one- dimensional equivalent to the top hat function. 
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Figure 3.1 The rectangular window function and its transform. Left: The rectan- 
gular window function w{n) in the time domain. Right: Magnitude of the Fourier 
transform of w{n), \W{k)\. w{n) is unity for n = [0,31], corresponding to 32 discrete 
samples in time (i.e. the digitized waveform), and zero outside this range. |iy(fc)| is 
unity for k = 16 and k = —15, and zero at all other integer values of k. See discussion 
in text. 

the discrete-time representation of the continuous, sinusoidal birdie signal (i.e. the 
digitally-sampled birdie, cross-correlated with a similar birdie from another antenna, 
idealized to a sine wave). 



xin) 



A sin 



2'KnkQ 
N 



A 



^j2nnko/N _|_ ^-j2TTnko/N\ 



(3.11) 



The DTFT (Eq. 2.7) of x{n) is: 



X{k) 



A 



[6{k - ko) + 5{k + ko)] 



(3.12) 



The signal x{n) is limited to the range n = [0, — 1] because we only have 
samples, and is zero elsewhere. Therefore, the signal is rectangular windowed by the 
limited number of samples. The transform W{k) of a rectangular window function 
in the time domain, w{n), which is equal to 1 in the range n = [0, — 1] and is zero 



elsewhere (see See Fig. 3.1), is 



kn/N 



-j-Kk(N~l)/N 



sin(7rA;) 
sin(7rA;/A^) 



(3.13) 
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Since multiplication in the time domain is equivalent to convolution in the frequency 
domain, the DTFT of x{n) multiplied by the window w{n) is 

Yik) = f2 Wik)X{k) = ^ [W{k - ko) + W{k + ko)] . (3.14) 

n — — oo 

The convolution of W{k) with the birdie X{k) simply shifts the W{k) to be centered 



on ±kQ (hence W{k ± ko) in Eq. 3.14). In the DFT, which only samples the DTFT 
at discrete values of k, if ko is a multiple of fs/N, we have an integer value of ko, 
meaning we sample W{k) only at values where it is equal to one or zero. For k = ko, 
Y{k) = 1, while Y{k) = for A; ^ ko. 

For non-integer ko, W{k) is shifted so that it is sampled at values other than zero 
and one. A birdie not centered on a spectral channel thus leaks into all channels 
within the band, with the strongest contamination appearing in the channels nearest 
the birdie frequency]^ It is easiest to remove birdie contamination from an observation 
if it is centered on a frequency channel; in this case, all the birdie power appears in 



one frequency channel, which can be excised during data reduction (see ^4.2). 

Since spectral leakage arises from the DFT itself, this leakage will not show up in 
any band other than that in which the birdie lies. This is simply because each band is 



sampled and FFT'd individually (see Fig. 2.20). Furthermore, increasing the number 



of channels (which could be done by increasing the number of correlator lags) would 
localize the leakage; in the limit of infinite channels, a birdie would be measured at 



its true frequency. This is the case where the DFT limits to the full DTFT (Eqs. 2.7 



and 2.8), which has infinite spectral resolution within a band fs/2 in size. 



^Note that, in the XF/lag correlator scheme, the FFT is taken on the cross-correlated data, hence 
the rectangular windowing of the signal only happens once, and frequency response goes as sine, not 
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SZA 9.022 GHz Birdie in Band 15 



0.5 



0.4 



0.3 



a 



0.2 



0.1 




I \ I 

62.5 125 187.5 250 312.5 

Frequency witliin the Band (MHz) 



375 



437.5 



500 



Figure 3.2 9.022 GHz birdie in the IF, downconverted to 1.022 GHz, sampled in time 
and FFT'd. The blue line is the amplitude of the bandpass that would be determined 
from the SZA spectral data. The birdie does not alias cleanly to the center of any 
channel, and therefore leaks into every other channel in the band. The DTFT of the 
birdie (red curve) is the sum of two sine functions, centered at 478 and 522 MHz, (or 
at ±478 MHz, since the transform of a the windowed sinusoid is W{±ko))- Values of 
k = [0, 16] from the DFT are plotted here, corresponding to the positive frequencies of 



Channels through 16 in the SZA correlator (see Table 2.2). The band wraps outside 
k = [0,31] (or k = [—15,16]), which corresponds to frequencies of 0-1000 MHz (or 
-500-500 MHz). 
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Retuned '9.03125 GHz' Birdie in Band 15 
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Figure 3.3 Retuned "9.03125 GHz" birdie (placing it at an integer k frequency), 
downconverted to 1.03125 GHz, then sampled and FFT'd. This birdie aliases cleanly 



to the center of the Channel 16 (see Table 2.2), which can be excised from the data 



The other channels' centers sample the nulls of the rectangular window's transform 



(red curve), which is a summation of sines located at ±468.75 MHz. See Fig. 3.2 for 
more details. 

3.1.4 Example of Spectral Leakage Corrupting a Band 

The 9.022 GHz birdie, arising from the 3'''^ YIG harmonic (detailed at the beginning 



of this section, ^3.1), corrupted all of Band 15 during initial testing of the SZA (see 



Table 2.1). Just before digitization, this birdie is at 1.022 GHz; this is not a discrete 



multiple of a 500/16 MHz = 31.25 MHz. The birdie's presence in the final data 
product (before the problem was solved) indicates that the bandpass filters of the 
block downconverter could not completely eliminate strong signals from outside the 
band. 

Sampling at 1 GHz aliases the birdie to 522 MHz. Due to the Hermitian property 
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of a real signal's FFT, frequency is reflected about 500 MHz. The 522 MHz birdie 
therefore appears in the 0-500 MHz band, between Channels 15 and 16 (468.75 and 



500 MHz). The effect of spectral leakage, shown in Fig. 3.2 , is calculated by evaluating 



the DFT expression for X{k) (Eq. 2.9), where x{n) was a sinusoid at 522 MHz. 
While the birdie is strongest in the channels adjacent to its real frequency, it remains 
above 3% of its full power throughout the entire band, rendering the band useless for 
observations. 

A clever retuning of the last downconverter LO ( ^2.5.1 ) seen by Band 15 avoids 
contamination of the entire band, localizing the birdie's power to one channel. The 
details of this are as follows: The 9.022 GHz birdie first ends up at 978 MHz after 
mixing with the 10 GHz LO of the block downconverter; it is then mixed with a 



2.0 GHz LO (Downconverter Bank 2, in Fig. 2.20), which places Band 15 at in the 
0.5-1.0 GHz band that is input to the digiters. The digitizers therefore see the 
birdie at 1.022 GHz. We retuned the last LO to 2.00925 GHz, shifting the birdie to 



1.03125 GHz. Fig. 3^ shows the spectrum that results from this retuning. The result 
is that birdie's power only shows up in Channel 16, which is set in the data reduction 
pipeline to be automatically flagged ([ |4.2[ ). 

The final solution to the birdie question was to eliminate the birdies at their 



sources: the YIG oscillators. By wrapping the YIG in Eccosorb® (see Fig. 2.17), we 
found the signal was sufficiently attenuated; this was a temporary solution, since the 
form of Eccosorb foam used in this test has a tendency to flake and degrade in the 
field. A more permanent solution was implemented by simply shielding the YIGs in 
small aluminum boxes. While the retuning of the 2.0 GHz downconverter LO was 
retained, periodic measurements continue to show that the YIG harmonics do not 
escape their aluminum housings (though we still flag Channel 16 of Band 15, as this 
precaution only costs us one of our 16 x 15 spectral channels). 



3.2 Antenna Cross-talk 
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In the survey fields, which were all selected to pass close to zenith in order to mini- 
mize atmospheric noise, an apparent excess in nois^was discovered between adjacent 
antennae when observing at antenna elevations > 80°. It was recognized that all the 
antenna positions that allowed this "excess noise" on these short baselines corre- 
sponded to positions where the secondary mirror feedlegs on one antenna had direct 
lines of sight to the feedlegs on the adjacent antenna. 

It was further discovered that, with fringe tracking off and the antennae in a 
position known to produce "excess noise," the measured phase versus frequency on 
that baseline (which for blank, untracked sky should be random noise) was strongly 
wrapping every two channels; this observed "excess noise" was actually a signal due 
to "cross-talk" between elements of the antennae, and the wrapping was due to the 
pathlength between the elements involved in the antenna cross-talk. Since each chan- 
nel is 31.25 MHz, the measured wrapping corresponded to a delay of 16 nanoseconds 
(i.e. 2 X 31.25 MHz), and the position of the antennae was such that the astronom- 



ical delay (T in [1.3) for that pair of antennae was also 16 nanoseconds. This is 
because projected baseline at the cross-talk maximum position was roughly 480 cm, 
or cT = (3 X lO^'' cm s^^) x (16 ns). Thus the delay used to track astronomical fringes 
on the baseline formed by the two antennae also flattened the phase of the cross-talk 
signal, providing a coherent, correlated signal]^ 

Several ways of eliminating this antenna cross-talk were explored. While the 
precise cause is unknown, it is hypothesized that the feed horn from one antenna 
could launch a signal that scattered coherently into the other antenna, and the other 



^The rms noise of the data when observing blank sky would slowly rise and fall while tracking a 
field through certain antenna configurations. This excess in noise was not correlated with changes 
in the atmosphere. 

^Flattening the phase of a signal across a band is the spectral equivalent to sending the proper, 
fringe-tracking delay necessary to correlate data. 
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antenna would do the same in reverse. With fringe tracking on, the delay being sent 
to an antenna sometimes precisely flattened this phase wrapping, producing strong, 
correlated amplitudes in the data. 

In order to test that the scattering came from the feedlegs of the secondary, a 
wire mesh was added to the side of one antenna from a pair that exhibited cross- 
talk. This shifted the cross-talk to higher antenna elevations. We considered building 
collars for some or all of the antennae in the inner array, but this was determined to 
be too costly, and the high winds in Owen Valley would make it difficult to maintain. 
Furthermore, if the collars extended too far, this would increase the risk of collisions 
between close-spaced antennae, something the SZA was designed not to have. 



We then tried adding Eccosorb® to the feedlegs (see Figure 3.4). While the 
absorber did eliminate cross-talk, it also increased the noise in the system, as each 
feedleg was now radiating at ~ 300 K. For comparison, the cluster signals we observe 
peak at ~ 10 mK. This increase in system noise was determined to be unacceptable, 
as it increased the required integration times. 

The final solution to the cross-talk problem was to wrap the feedlegs in crinkled, 
heavy-duty aluminum foil (~50 mil), available at most grocery stores. The crinkles 
were carefully constructed to produce large, 3-10 cm facets of arbitrary orientations 
on the side of the feedleg facing the primary mirror. These facets have the effect of 
randomly changing the pathlength and angle between corresponding parts of each 
antenna. 
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Figure 3.4 Photo of the closest pair of antennae (Antennae 3 & 5, see Fig. 1.7), with 
the test Eccosorb® wrapped around the feedlegs. Uhimately, the Eccosorb® was 
replaced by crinkled aluminum foil. 



Chapter 4 
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Data Reduction 

In this chapter, I provide a brief overview of the steps necessary to cahbrate SZA 
data. For a more detailed treatment of the SZA data analysis routines, I refer the 



reader to the thesis of Stephen Muchovej (Muchovej 2008) 



The SZA data calibration routine was written by students in the SZA collabo- 
ration (primarily Mike Loh, Stephen Muchovej, Matthew Sharp, and Chris Greer) 
m MATLAB®Q code. Since MATLAB® is a scripting language optimized for matrix 
manipulation, it is well-suited to manipulate large data sets when performing uni- 
form mathematical operations. Some of the more computationally-intensive routines 
- in particular, those that required loops - were implemented in external, compiled 
C/C++ functions that are called from within MATLAB®. In the following sections, 
I describe each step in the data reduction pipehne in order of operation 

4.1 System Temperature Computation 

The first step in the data calibration routine is a calculation of the system temperature 



Tgys (see Eq. 2.2). As described in ^2.3.1, T^ys is defined similarly to Trx, but accounts 
for the noise of the entire system scaled to above the Earth's atmosphere. We therefore 
use Tgys to compute the rms noise level in an observation, 

Ts s 

arms OC — (4.1) 

^N{N-l)TAiy ^ ' 



http : //www . mathworks . com 
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Here is the number of antennae in the array, r is the integration time, and Az/ is 
the bandwidth of the system. For a given integration time, on a given instrument. 



Tgys characterizes the noise in an observation. 



During a typical cluster observation, we compute the noise temperature while 
observing the phase calibrator and while observing the target (usually a cluster) 
separately, since each is observed through a different column of atmosphere. See 



Fig. 4.1 for a plot of Tgys in an observation. 

Since Tgys is used to compute the theoretical rms noise in the observation at a later 
stage of the reduction pipehne (see ^4.3.4[ ), we interpolate our T^ys measurements past 
large deviations from ideal behaviorj^ This ideal noise is used in flagging data that 
show large deviations from the theoretical rms noise. A jump in Tgys of more than 
~ 1.5 K can indicate a sudden change in weather or an instrumental effect, so we 
flag and interpolate stable Tgys measurements past these jumps. Additionally, data 
are flagged if no Tgyg measurements are available for over ~ 30-35 minutes. This was 
determined during pipeline testing to be the longest reasonable time to go without a 
Tsys measurement, as we are effectively interpolating the behavior of the atmosphere 
between actual measurements. 



4.2 Automatic Flagging 

Some bad data are next automatically flagged, without any user intervention, on a 
number of conditions. These conditions can depend on the state of the array, in which 
case they are flagged by the control system (i.e. for these instances, the reduction 
routine simply propagates the flags set by the control system). Examples of the 



^As discussed in S 2.3.1 



determining Tgys requires a measurement of total output power (i.e. not 
correlated visibilities, but an actual measurement of power from each antenna) while observing two 
thermal loads with known temperatures. Therefore, Tgys cannot be determined while observing the 
source, and must be interpolated between measurements, which are taken every ^ 10 minutes. This 
interpolation provides an estimate of Tgys for each data point. 
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Figure 4.1 System temperatures of a typical SZA cluster observation, computed for 
each band. One representative antenna is shown, since the other seven all have similar 
system temperatures. The y-axis shows Tgys in K, and the x-axis is time in hours from 
the start of the observation. The dark blue point at the start of the observation is 
the Tgys of the bandpass calibrator (^4.3.1). The black points are the Tgyg for the 
target data (usually a cluster), and the orange points are T^ys for the phase calibrator 
(^4.3.2). The points differ due to the differing columns of atmosphere to each source. 
The minimum in each curve occurs when the source transits (reaches its highest 
elevation in the sky, and therefore the optical depth reaches its minimum for the 
observation) . 
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common flags set by the control system are: 

• Tracking (physical): If a telescope is not mechanically tracking a source for 
any reason (hardware malfunction, an antenna is off-line for repair, etc.), data 
from that antenna are flagged. 



• Fringe Tracking: If the software that implements fringe tracking (^1.3) is not 
functioning properly while observing an astronomical source, there can be no 
correlated data. Therefore, these data are flagged. 

• Correlator Bands not received: If the correlator software is not running for 
a given band, no correlated data are returned in that band; that band will be 
flagged. 

Other automatic flags can be due to properties of the array that can either be 
computed or have been measured. These are: 

• Antenna Shadowing: For given antenna positions and primary mirror sizes, 
in a homogeneous array like the SZA, it is straightforward to compute whether 
one antenna obstructs the view of another. We use a projected antenna spacing 
of 3.6 meters, which treats the effective primary diameter also as 3.6 m in this 
calculation (recall that the SZA primary diameter is 3.5 m). This buffer is 
used to avoid diffraction effects at the edges of the primaries, which are due 



to the sidelobes of the antenna beam pattern (see Fig. 2.2). If an antenna 
is shadowed, all data from that antenna are flagged during that period (i.e. 
we flag each baseline involving the shadowed antenna for the entire time it is 
shadowed) . 



Birdies: As discussed in ^3.1, radio interference was found to contaminate 



some of the channels in our observations. A list of bad channels was kept with 
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the rest of the data reduction code, in a Concurrent Versions System]^ (CVS) 
archive. 



• Cross-talk: As discussed in ^3^, coherent signals between closely-spaced an- 
tennae could corrupt baselines. Baselines measured to have significant amounts 
of cross-talk were noted in the CVS, and those corrupted baselines were ex- 
cluded. 

After the automatic fiagging, a baseline correction is applied to the data. By ap- 
plying small corrections to the antenna locations used by the control system, slightly 
different astronomical delays are computed for each antenna as the source traverses 
the sky. The set of correct antenna locations, which produce a constant phase versus 
time on each baseline, is called the "baseline solution." 

We typically measure the baselines every two weeks by observing many point 
sources across the sky over the course of a few hours, and permuting the old baseline 
solution to find the fiattest phases for the set of point sources. Common reasons for 
changes to the baseline solution, which are typically less than a millimeter, involved 
expansions and contraction of the ground due to rain, or freezing and thawing of 
moisture in the soil between the concrete antenna pads. No fiagging is performed 
when the baseline solution is applied, as this is simply a mathematical correction 
that is applied to the data. 



4.3 Interactive Data Calibration 

The next step in the pipeline is a series of calibrations that can be run interactively 
or by specifying scripts with desired parameters on which to fiag. While the following 
steps are now generally run in an automated fashion using scripts supplied by the 
user, we retain the fiexibility to fiag data by hand. This manual mode of fiagging was 

" http : / / www . nongnu . org/ cvs/ 
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particularly important when testing and debugging the pipeline, and was necessary 
to establish reasonable limits to use in script-driven data calibration. 

4.3.1 Bandpass calibration 

Any component - an amplifier, a cable, a filter, or a mismatch between compo- 
nents or sections of waveg uid^g can add a spectral shape to the signal that is 
non-astronomical in origin. By observing an astronomical source known to have a 
simple power law spectrum, instrumental phase and amplitude features can be cali- 
brated out of the data. Bandpass calibration is done by solving for the corrections 
necessary to fiatten both the amplitude and phase of each band's spectrum, on a per 
antenna basis. Since these instrumental features do not change on the timescales 
of the observation (and typically only do change when a component in the receiver 
or back-end electronics is changed), the bandpass calibration is performed once per 
observation, and can often be shared between adjacent observations. 

The first interactive fiagging step is therefore the bandpass calibration, in which 
the shape of the instrument's frequency response is removed from the channel-based 
data before averaging the bands. Recall from §2.5.1 & 2.5.2 that there are 15 channels 



in each of the 16 bands. Data are bandpass calibrated using an ^ 10 minute integra- 
tion on a strong point source (> 10 Jy). This short bandpass calibration observation 
is taken once per observation, at either the start or finish of an observation ("track"). 

Figures |4.2| & |4.3| show the bandpasses for each band of Antenna 5's amplitude 
and phase, respectively, before calibration. Figures 4.4 & 45] show these bandpasses 



after calibration. After the bandpass calibration has been applied, the channels of 
each band are binned, and we no longer work with the channel-based data. This 



binning reduces the size of the dataset by a factor of 17 (see Table 2.2). 



waveguide mismatch produces standing waves in frequency across a spectrum, due to reflec- 
tions. 
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Figure 4.2 Amplitude vs frequency channel of the bandpass calibrator, measured for 
each band of Antenna 5, before calibration. The y-axis is the amplitude in Jy, and 



the X-axis is the channel number. Fig. AA shows the calibrated amplitude of the 
bandpass. 
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Figure 4.3 Phase vs frequency channel for the bandpass cahbrator, measured for each 
band of Antenna 5, before cahbration. The y-axis is the phase in degrees, and the 



X-axis is the channel number. Fig. 4.5 shows the calibrated amphtude of the bandpass. 
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Figure 4.4 Amplitude (vs frequency channel) of the bandpass calibrator for each band 



of Antenna 5, after calibration. See Fig. 4^, which shows these bandpass amplitudes 
pre-calibration, for caption. 
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Figure 4.5 Phase (vs frequency channel) of the bandpass cahbrator for each band 



of Antenna 5, after cahbration. See Fig. 4.3, which shows these bandpass phases 
pre-cahbration, for caption. 
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4.3.2 Phase/ Amplitude calibration 

The next step in the cahbration corrects for drifts in the overall amplitude and phase 
of each band. Every ~ 15 minutes, a moderately strong (> 2 Jy) point source 
is observed for 3-5 minutes. An unresolved source at the array's pointing center 
(i.e. at the spatial location (0,0)) is mathematically a 2-dimensional Dirac delta 
function, 5(0, 0)|^ Therefore, its spatial Fourier transform is a purely real constant, 
which implies that each baseline ideally would measure the same amplitude, with zero 
phase. Deviations from zero phase, therefore, are tracked and corrected using this 
source as a "phase calibrator." 

Figures 4.6| -4.9 show the baseline-based phases and amplitudes of the phase cali- 



brator, measured by the example Baselines 0-1 and 5-6. By setting one antenna per 
band as the "reference antenna," which for that band is assumed to have zero phase, 
the calibrator phase for every other antenna can be computed. The reference antenna 
for each band is simply chosen to be the first antenna with no flagged calibrator data. 



In Fig. 4.11 , it is clear - from the fiat phase that is identically equal to zero - that 



Antenna was used as the reference for all bands except Bands 6 & 7, which had 



flagged phase calibrations mid-track (see red points in Fig. 4.6) 



The phase calibrator is used to track changes in the complex antenna-based gains 



throughout the track , which can happen for a number of reasons (see Figures 4.10 & 



4.11). The SZA receivers are well-characterized and have stable gains under normal 



^The integral of the point source's flux density over space is the total flux / of the point source 
(i.e. / = j S{x,y) dx dy). The general equation for a point source in M,w-space with a flux /q, 
normalized at frequency I'o = 30.938 GHz (center of the SZA 30-GHz band), and with a spectral 
index a, at image location {x,y) (in radians from the phase center), is 

/(m, u) = /o( — I [cos(2Tr(ux + vy)) + j sm(2n(ux + vy))]. (4-2) 

For {x,y) — (0,0), the point source is purely real, and has a constant flux at all locations in u,v- 
space. The simplicity of this solution motivates the use of point sources in tracking the system's 
response, as each baseline will ideally measure the same flux. Another way of stating this is that 
a point source is smaller than the beam formed by any baseline pair, so the flux probed by any 
baseline, in Jy/beam, equals the total flux in Jy. 
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Figure 4.6 Phase of the point source cahbrator versus time, as measured by a short 
basehne (Basehne 0-1). The y-axis is the phase in degrees, and the x-axis is time 
in hours since the start of the track. Band was entirely flagged for this track, 
due to broken digitizers and a lack of spares at the time. The magenta points are 
those flagged due to the automatic flagging, and the green points are those flagged 
by a user-specified script's limits, which catches 30° outliers from the underlying, 
interpolated phase. 
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Figure 4.7 Phase of the point source cahbrator versus time, as measured by a long 



basehne (Basehne 5-6). See Fig. 4.6 for further details. Note the slightly larger 



scatter in the phase on this baseline than on Baseline 0-1 (Fig. 4.6). This is due 
to atmospheric coherence being slightly poorer on long baselines (the antennae are 
looking through different columns of air, and atmospheric turbulence has a scale size 
on the order of tens of meters). The underlying slope in the phase calibrator is easily 
determined in this observation, as the atmospheric coherence for the example track 
was typical of a clear spring day. 
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Figure 4.8 Amplitude of the point source calibrator versus time, as measured by a 
short baseline (Baseline 0-1). The y-axis is the amplitude in Jy, and the x-axis is time 
in hours since the start of the track. Since a point source is unresolved (the beam is 
larger than the point source), the amplitude in Jy equals that in Jy/beam. This is 



the amplitude of the calibrator corresponding to the phase shown in Fig. 4.6 
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Figure 4.9 Amplitude of the point source calibrator versus time, as measured by a 
long baseline (Baseline 5-6). This is the amplitude of the calibrator corresponding to 
the phase shown in Fig. |4.7[ See Fig. 4.8| for further details. 
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operating conditions; changes in the complex gain of the actual instrument therefore 
usually arise when equipment is thermally unstable (failed refrigerator, poor regu- 
lation of the electronics box, rapid heating/cooling of the correlator trailer during 
sunrise/sunset, etc.). Sharp changes in phase, which lead to degradation of the mea- 
sured amplitude, are used to flag data. 

Scatter in the phase due to poor atmospheric coherence can lead to a lower mea- 
sured amplitude in the correlated data. Clouds and atmospheric turbulence during 
inclement weather are generally of a small enough scale that a baseline in the inner 
array measures the full amplitude of the calibrator, as the phase scatter is common- 
mode on short baselines (i.e. they see the same fluctuations); during such a period, 
long baselines measure a diminished signal, since the atmosphere is essentially adding 
random, uncorrelated phase fluctuations to each signal. The example observation pre- 
sented throughout this chapter was taken on a clear spring day with good atmospheric 
coherence. 



4.3.3 Orphaned data 



After phase calibration is complete, periods when data are not bracketed by Tgys or 
phase calibrator observations are flagged, as those data cannot be calibrated. We 
call these data points "orphans" and discard them. This includes data orphaned by 



flagged calibrator observations (such as the red points in Figures 4.6 and 4.7). Note 
that this step is antenna-based, since data are typically orphaned by antenna-specific 
problems (e.g. an antenna was shadowed, or was not tracking and missed a phase 



calibrator observation). See Figures 4.12 & 4.13 for plots that show the flagging of 
orphaned data. 
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Figure 4.10 Antenna-based gain of the point source calibrator versus time, for Antenna 
0. The y-axis is the gain (ideally 1), and the x-axis is time in hours since the start of 
the track. The scatter seen here, which is simply due to noise in the measurement, is 



smaller than our absolute calibration uncertainty (~ 5%, see Muchovej et al. (2007)). 
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Figure 4.11 Antenna-based phase of the point source cahbrator versus time. One 
antenna must be chosen as a reference in order to compute the antenna-based phase 
of the other antennae. In this plot, Antenna served as that reference for all but 



Bands 6 & 7, which had flagged calibration observations (see red points in Fig. 4.6) 
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Figure 4.12 Orphaned data from Antenna 0. Automatically flagged data in Band 
are shown as green points, and were due to a digitizer hardware problem (and 
a temporary lack of replacement digitizer boards). Data orphaned by gaps in the 



calibrator data in Bands 1, 6, & 7 are shown in red (see Fig. 4.11). Useful data are 



represented as blue points. For simplicity, target data are plotted with zero phase 
(rather than plotting the noisy distribution of raw target data phases, since we are 
only trying to determine which data are not bracketed by calibrator observations). 
The antenna-based phases of bracketing calibration observations are shown as black 
X's (which are identically zero in bands where Antenna was the reference; see 
Fig.liTTl). 
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Figure 4.13 Orphaned data from Antenna 5. Data orphaned by 



data in Band 2 are shown in red (see Fig. 4.7). See Figure 4.12 



Antenna 5 was not the reference antenna for any band in the 
calculation, so none of its phases are identically zero. 



gaps in the calibrator 
for details. Note that 
antenna-based phase 
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4.3.4 Noise and variance of the data 

The next step in the user-guided cahbration is to flag target data with a high or 
biased level of rms noise. The real rms noise of the target data is compared with 



the expected theoretical value, computed using the measured values of Tgys (see ^4.1 



Eq. 4.1). If the variance of the data does not agree with the theoretical value on the 
~ 35% level, the data are flagged. This threshold was determined by Matthew Sharp 
in a series of jackknife tests (differencing of essentially blank fields) on many CMB 
fields. Since Tgys data have been cleaned of outlier measurements, this flags periods 
with higher than expected noise. 

Additionally, any data points with more than a 4-0" deviation from the distribution 
are flagged. While this filter is useful in principle, data with such large error bars are 
typically caught in the previous steps. See Figures |4.14 & 4.15 for plots of the rms 



noise of target data for a long and a short example baseline. 

The expected rms noise is computed using cleaned Tgys data, where any jumps in 
Tsys have been interpolated past. However, if the jump was due to a real instrumental 
or atmospheric effect, the data will show a similar jump in measured rms noise. This 
is why large deviations from the expected theoretical noise are flagged. An example 
of an instrumental effect that leads to a higher than expected rms level (without 



affecting the measured Tgys) is antenna cross-talk, discussed in ^3.2 



4.3.5 Amplitude of the target data 

The last step of user-guided calibration flags for anomalies in the amplitude of the raw 



target data. For the example cluster observation (shown in Figs. 4.16 & 4.17), there 
are no strong (> 1 Jy) point sources in the field. Individual data points are dominated 
by noise, since the cluster signal is on the order of ten mJy (total, integrated over the 
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Figure 4.14 rms of target data, taken on short Baseline 0-1. The theoretical predic- 
tion, based on the measured Tgys, is shown in black. Flagged data - already caught 



by other steps (see Figures 4.6 & 4.12 in particular) - are plotted in red and magenta. 



Newly flagged data are in green. Useful, unfiagged data are in blue. 
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Figure 4.15 rms of target data, taken on long Baseline 5-6. Flagged data - already 
caught by other steps (see Figures 4.7 & 4.13) - are plotted in red and magenta. See 
Figure 4.14 for more details. 
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sky), and the rms noise per second is ~ 0.5 Jy/beam|j Data points with more than 
a 4-cr deviation from the mean - which could be due to high atmospheric noise - are 
typically flaggedj^ Also, any jumps in target amplitude are immediately suspect, as 
they are indications of undesired instrumental/atmospheric effects (since the cluster 
SZE signals and point source fluxes are not time-dependent on the scale of hours). 

Though we never attempted observations of clusters with (> 1 Jy) point sources 
in the field, the CMB and cluster survey data often do contain such sources. Their 
presence simply raises the mean value of the amplitude of the target data. Since we 
only flag on outliers from the distribution, note that this step does not prevent us 
from performing observations of fields with strong point sources]^ 



4.4 Dirty Maps 



After all steps of the calibration are complete, the machinery of Difmap (Shepherd 



1997), imported into MATLAB®, is used to produce rudimentary maps of the obser- 



vation (shown in Fig. 4.18). These maps are largely illustrative, as all quantitative 
results of the SZE cluster observations are determined by fitting the visibilities in 



UjV-space (see ^5.4). We do this because the noise properties of the visibilities in 



-space are well understood. 



In Fig. 4.18, the short baseline data (upper panel) show a central, dark blue 



"blob," which is the flux decrement due to the cluster CL J1226. 9-1-3332. To the 
left of the cluster, which is east on the sky, a ~ 4 Jy source is seen as a red "blob" 



^This is for typical binning of the visibilities - of course the signal-to-noise ratio (S/N) increases 
when data are binned, revealing a high-significance detection of the cluster's SZE signal, as well as 
the accompanying point sources in the field. 



^While data with larger, random rms noise would have a lower weight (see { 4.5 1, and thus would 



not bias datasets, including them increases the size of the calibrated dataset. Since noisy data points 
are the exception, it is easier to exclude them. 

®Since the SZA's dynamic range is ~ 100 the presence of a 1 Jy point source in a field would 
prevent us from obtaining good constraints on a ^ 10 mJy cluster SZE signal. 
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Figure 4.16 Amplitude of visibility data (in Jy, see |:1.3) for a short baseline (0-1) 



taken on the target source (typically a cluster). Blue points are unflagged data. Red 
and magenta points are flagged data caught by previous steps in the data calibration 



see Figures 4.6 & 4.12), while green points may indicate data flagged in this step as 



outliers or from previous steps (e.g. Bands 1, 2, & 14 contain some green points that 



are not outliers, flagged in Fig. 4.6). The SZE flux from the cluster and the fluxes 



of point sources in the cluster fleld are on the ~ mJy level; they are therefore not 
noticeable in the raw data plotted here (which are noise dominated). 
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Figure 4.17 Amplitude of visibility data in Jy (see ^1.3) for a long baseline (5-6) taken 



on the target source. See Fig. 4.16 Note that long baselines do not typically measure 



cluster scales, but the raw target data here and in Fig. 4.16 have similar amplitudes 
Also note that some green, flagged points (e.g. in Band 13) were flagged in Fig. 4.7 
not because they are outliers. 
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(attenuated by the primary beam to ~ 2.5 Jy, as it is 4.5' from the center). The 
resolution of the upper panel map is determined by the scale probed by the short 
baselines, and not the size of the source in the sky (if larger, the baseline would not 
be sensitive to the source). The lower panel shows a slight excess in flux at the point 
source position. The long baselines do not probe cluster scales. For well-calibrated 
data, only the very strongest sources show up in a single track, and a typical cluster 
observation requires more than a single track to constrain a cluster's SZE signal (as 
well as the weak ~ few mJy point sources in the field). We therefore expect a single 
track to be dominated by noise. 



4.5 Final Data Product 

Calibrated data are output, using a C++ function written by Erik Leitch, in the stan- 
dard Flexible Image Transport System (FITS)|^ format for astronomical data. Since 
these are radio interferometric data in n^f -space, we use the UVFITS specification 
of the FITS format. These files contain the real and imaginary components of the 



visibility data (in Jy, see ^1.3) from all 16 SZA bands, stored in binary format. They 



also contain the -space coordinates and weights (inverse variance, 1/cr^ms) GSich 



data point 



10 



Further useful information is also stored in the UVFITS file. An ASCII header 
provides observation specifics, such as the date of the observation, name of the ob- 
servatory, name of the target observed, the pointing center used for the target, and 
the frequency at which it was observed, while a binary table includes each antenna's 
physical location. 



3 See 



http : //f it s . gsf c . nasa . gov/ 



"'^'^Interferometric maps of the M,?;-data typically weight the data by inverse variance (1/cr. 
Information about the noise is stored in this form to avoid redundant computation steps. 



^ ) 

rms/ 
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Figure 4.18 Dirty maps of the short (upper panel, < 2 kA) and long (lower panel, 
> 2 kA). X and y axes are the map coordinates in degrees. The colors represent the 
signal in Jy/beam (the flux detected within the beam formed by each baseline). See 
text for more details. 
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Part II 



Modeling and Analysis of Clusters 
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Chapter 5 

Modeling the Cluster Signal 

5.1 Introduction to Cluster Models 



As detailed in ^1.3.2 the SZA interferometer probes a cluster's SZE signal in u,v- 



space (see Figures 1.5-1.11). In order to extract physical parameters from an in- 
terferometric cluster observation, we need to fill in information on the scales we do 
not probe. A simple model to extrapolate through the individually noisy, sparsely- 



sampled visibilitie^ (Eq. 1.6), is used to extract meaningful constraints on physical 



cluster parameters. Given the large gaps in the w,w-coverage, a simple, 2-parameter 
model was found to work best. 

When fitting data in u,v-space, it is helpful to consider two extreme cases: scales 
below the resolution of the instrument (i.e. smaller than the scales probed by the 
longest baseline), and scales larger than the shortest baseline probes. The first case 
defines what we consider to be a point source, where in u,v-spa.ce each antenna ideally 
measures constant, uniform power. 

The second case implies that we are not sensitive to backgrounds, such as the 
primary CMB (most of its power is on much larger angular scales than we probe) or 
large-scale galactic emission (from our galaxy). This makes the SZA a good spatial 
filter for clusters at high redshift, since it interferometrically probes Fourier modes of 
the spatial intensity pattern corresponding to the bulk of a cluster's signal. It also 
means we cannot constrain cluster features on scales larger than the shortest baseline 
can probe. A cluster model can predict an arbitrary signal at scales larger than we 



^There are typically ^10-100 thousand visibilities in an SZA cluster observation, yielding high 
signal-to-noise SZE observations when the data are binned. 
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can probe; a poorly-motivated model can have Fourier modes in its transform that 
agree with the measured SZE signal, but might predict too much or too little flux 
on scales not accessible to an interferometer. At the longest wavelength of the SZA, 
observing down to the shadowing limit, this upper limit is ~ 12.7'; above this scale, 
information is lost due to lack of -coverage. 

Modeling of any signal, including that due to the SZ effect in clusters, introduces 
priors on the resulting fit and its derived quantities. I discuss below the SZE models 
I tested, as well as the X-ray models used to complement them. I then discuss the 
Markov chain Monte Carlo (MCMC) fitting procedure we use to fit these models 
jointly to SZE+X-ray data. Finally, I show how we derive galaxy cluster properties 
from the resulting model parameters. 

Throughout this section, we use several common cluster overdensity radii that 
appear often in the literature]^ These overdensity radii are defined with respect 



to the critical density of the Universe Pc{z) at the cluster's redshift (see Eq. 5.35). 
We define as the radius within which the average (gas+dark matter) density is 
A pc{z) (e.g. r25oo is the radius containing an overdensity 2500pc(^)5 and v^qq is the 
radius containing 500pc{z)). Since the mean cluster density drops with radius, a 
smaller overdensity A corresponds to a larger radius. A good rule of thumb is that 
rsoo ~ 2 r25oo for most massive clusters. We show how we solve for these radii for real 



observations in ^5.5.6 



^The "overdensity" is simply the factor by which the average density in a given volume is higher 
than the critical density of the Universe at that redshift, which is the density it takes for an object 
to begin to collapse, rather than expand with the Universe. 
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line of sight 

Figure 5.1 Geometry for the line of sight integral of a spherically-symmetric model. 

5.2 /3-Model SZE Profiles 
5.2.1 The /3-model 



The /9-model, also called the King profile, has a long history (Cavaliere & Fusco 



Femiano 1976 1978) of performing remarkably well at parameterizing cluster radial 



gas density profiles. Part of its appeal is its simplicity. Here I show how to derive 
the line of sight Compton y parameter when a /3-model is used to describe the radial 
electron pressure profile, 



Pe(r) = Pe,o [1 + [r/rcf] 



21 -3/3/2 



(5.1) 



Here Pe,o is the central electron pressure, /3 is a single slope that describes how pressure 
decreases with radius, and is the "core radius," typically ~ 100 kpc for a massive 
cluster. For /3 = 2/3 - a typical value when this model is used to describe the intra- 
cluster medium's (ICM) density - is the radius at which the electron pressure falls 
to half its central value. 



Ill 



Figure 5.1 shows the geometry for the hne of sight integral of a spherically- 
symmetric ICM model. Here, r is the radius within the cluster, R is the radius 
on the sky (from the projected cluster center on the sky), and p = ^r'^ — B? is the 
distance along the line of sight. The radius on the sky, i?, is a constant in the inte- 



gration, equal to the minimum in f (see Figure 5.1). The Compton y parameter as a 
function of R is 



nip a 



/oo 
Pe,o [1 + (r/rj2] 
-oo 



21 -3/9/2 



Vr^ - R^ 



dr. 



(5.2) 



Substituting in p and dp/dr = r/y/r'^ — R'^, and taking advantage of the symmetry, 
we have 



poo 

y{R) = 2Pe,o^ / [1 + {p/Tc)' + {R/rc)'] dp. 



Vile C 

A rearrangement of terms yields 



,(i?)=2P,o^ [l + {R/ref]-'"' 



(5.3) 



p 



rl[l + {R/r, 



-3/3/2 



dp. (5.4) 



Defining the quantity pc = \Jr1 + R^, the integral in Eq. 5.4 becomes integrable in 
terms of the hypergeometric Eulerian function F. 



1 + 



p" 



r2[l + (i?/r,)2] 



-3/3/2 



■ /"OO 

dp= / [l + {p/Pc)T'^"dp 
Jo 

F(3/3/2 - 1/2) 



'Trpc 
V7r(r2 + ^2) 



F(3/3/2) 
F(3/3/2 - 1/2) 
r(3/3/2) 



(5.5) 



^,n + (^/,)Y/2r(W^^ 

vrr^i + ^it/rj j ^^^^^^^ 



Putting the result of Eq. 5.5 back into Eq. |5.4[ we finally have 



2 



TT re [1 + (R/r, 



N2]-3/3/2-hl/ 



r(3/3/2 - 1/2) 



F(3/3/2) 



(5.6) 
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Note that Eq. |5.6| has the same form as the integral of the hne of sight through the 



center of the cluster. We therefore can define the central Compton y parameter uq: 

ar r- r(3/3/2-l/2) , , 

This yields the expression for the SZE /3-model commonly fit to observations: 

= yo [1 + (i?/re)r''/'^'/'- (5.8) 



In terms of temperature decrement]^ AT/Tcmb = y (as in Eq. 1.1 ), we also define 
the central decrement 

ATo = /(x) I/O ^CMB, (5.9) 



where fix) contains the frequency dependence of the SZE (see Eq. 1.4) 



The /3-model thus provides a simple, semi-analytic formula for fitting SZE data. 
5.2.2 The Isothermal /3-model 



In we made no assumptions about the temperature, density, or metallicity 

distribution within the ICM; we only assumed the electron pressure can be described 
by a /3-profile. 

For the following, we make the common assumption of isothermality. We will 
also assume the X-ray-determined isothermal spectroscopic temperature Tx is the 



temperature of the ICM gas (Te(r) = Tx\ Using Eq. 5.1 and applying the ideal gas 
law, 

Pe = kBneT,, (5.10) 



^Below the null in the SZE spectrum, of course. Above ^218 GHz, we would refer to the 



■'temperature increment" instead. See Q.2 
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we arrive at this expression for the electron density: 

r^e(r)=r^,,o[l + (r/r,)2]-=^^/^ (5.II) 

A number of joint SZE+X-ray cluster studies have been performed using the 



Grego 


1999 


LaRoque et al. 


2006 


Bonamente 



et al. 2008). Throughout these studies, the isothermal /3-model has been applied 



in various ways; for simplicity we only discuss the most recent. In LaRoque et al. 



(2006) and Bonamente et al. (2008), the jointly-fit SZE and X-ray data were used 



to constrain an SZE-determined central electron density, ?T-e,o,sz- This comes from 



applying Eq. 5.10 to the SZE-constrained central pressure Pe,0; and using the X-ray- 
constrained temperature Tx to solve for the central density (i.e. They used Ue^^sz = 
Pefi/ksTx)- The analogous X-ray-determined central electron density, which does 
not require SZE data, is called ne,o,x in this context (and later is just called ne,o, as 
we never use Uefl^sz in this work). By not requiring riefi^sz = n^fi^x to hold when 



jointly fitting SZE+X-ray data, the central decrement (Eq. 5.9) can be allowed to 
fit freely the amplitude of the SZE signal. These different normalizations for the 
electron density profile were used to provide somewhat independent constraints on 
the hot gas fraction, using the X-ray-constrained estimate of the total cluster mass 
with the X-ray or SZE-constrained gas mass estimates (we show how these quantities 



are derived from the fit profiles in ^5.5). An important drawback of this method 
is that, since X-ray imaging data have a much higher spatial significance than SZE 
data, the jointly-fit shape parameters Vc and /3 are driven to those values preferred 



by the X-ray data (as we will demonstrate in ^6.6). Furthermore, it has been well- 
established by X-ray observations that isothermality is a poor approximation both 
within the cluster core and at large radii (see e.g. |Piffaretti et al.|2005 Vikhlinin et al. 



2005 Pratt et al. 2007, and references therein)). 
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Since current X-ray imaging data have higher signal to noise than current SZE 
observations, and X-ray imaging has been shown to recover well cluster gas density 



(Nagai et al. 2007b), we will instead only consider, for comparison, what the SZE fit 
yields for SZE-specific quantities, such as electron pressure and the integrated, intrin- 



sic Compton-y parameter, Yint (see ^ 5.6.1). We use the X-ray-determined density and 
temperature for all other quantities derived from cluster observations that have been 
fit with the isothermal /3-model, and do not attempt to constrain electron density 
with SZE data. 

As noted above, the assumption of isothermality does not remain valid over a 
broad radial extent of a cluster, even when the core is excluded from the X-ray data 



(as was done in LaRoque et al. (2006) and Bonamente et al. (2008), for example). 
Early joint SZE-|- X-ray data fitting only allowed the SZE data to constrain one unique 
parameter - yo - to within ~ 10% error bars, while X-ray data primarily constrained 
the cluster shape parameters {/3 and r^, used to fit both the SZE and X-ray data) and 
density normalization ne,o|^ Since the SZA can integrate to the OVRO/BIMA signal 
levels ~ 8 times more quickly, the resulting higher quality data of a typical 30-GHz 
SZA observation allows two parameters to be uniquely determined by SZE data. An 
additional benefit of the SZA over OVRO/BIMA is in the shorter baselines, due to 
the close-packed inner array (discussed in ^2.1 ). The short baselines of the SZA allow 
us to probe larger scales, where isothermality is not expected to hold. 

Note that we present the results of the isothermal /3-model jointfits for comparison 
only, and use its limitations to demonstrate how a new model, presented next, can 
move beyond the /3-model. 



""^Using X-ray data alone, temperature and density can both be determined. Assuming the ideal 
gas law, yo can be computed using only X-ray data. However, yo was allowed to be fit independently 
by SZE data in earlier works. In order to assess the performance of the /?-model as it was applied 
previously, I did not change how this fit was performed. 



5.3 A New SZE Pressure Profile 
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5.3.1 Motivations for a New Pressure Profile 



We motivate here the use of a new pressure profile, chosen not to have shape parame- 
ters hnked to those used in fitting the X-ray imaging data, which primarily determine 



ICM density (note the ^^(r) dependence in Eq. 1.8 for the X-ray surface brightness 



implies a strong density dependence)]^ By chosing a new pressure profile, we explic- 
itly allow the electron pressure and density profiles, Pe{r) and ne(r), to have different 
shapes, consistent with a non-isothermal electron temperature profile Te(r). In order 



to motivate this new pressure profile, first presented in Nagai et al. (2007al hereafter 



N07), we examine the empirical evidence afforded by recent, detailed X-ray studies 
of galaxy clusters. It is known that the isothermal assumption does not hold over a 



broad range of cluster radii (see, e.g.. Figure 5.3). In this section, we discuss some of 



the background work that provides the foundations of the new models. 

In looking at clusters over a broad range of masses, we hope to identify whether or 
not they exhibit self-similarity, which (if self-similarity holds) implies that they tend 
toward the same shape or slope when scaled to a fiducial radius such as the virial 
radius or rsoo- While several works have demonstrated that temperature can exhibit 



a nearly self-similar decline at large radii (r > 0.2r^ij., see e.g. Loken et al. (2002); 
Markevitch et al. ( 1998[ )), it is also clear in these works and in recent detailed X-ray 
measurements that cluster temperature profiles are not universal or self-similar over 



all cluster radii (particularly within the cores; see Vikhlinin et al. (2005, 2006) and 



Fig. 5.3). 



Recently, Vikhlinin et al. (2006 , hereafter V06) demonstrated a new density model, 
which they used to fit the surface brightness profiles of 11 nearby, relaxed clusters of 



1.^ 



^The X-ray surface brightness in Eq 
plasma emissivity A(T) a T^/^. Therefore 
than it is to temperature. 



behaves as Sx oc {1 + z) J Aee{Te, Z) d£, where the 
the surface brightness is much more sensitive to density 
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Figure 5.2 The above figure shows fits to the density profiles of 11 nearby, relaxed 
real clusters, as well as to the average density profiles of 16 clusters simulated using 
adiabatic (red) and cooling + star-formation (blue) physics. V06 motivate the use 
of a nine-parameter density model (Eq. |5.12[ ) to model these density profiles. Figure 
from NO 7. 



galaxies with deep Chandra X-ray exposures. The density profile V06 fit is 



ne{r) 



n 



eO 



(r/rj- 



+ 



n; 



e02 



[1 + (r/r,)2]3/3-"/2 ^ (^/^^)7]-/^ [1 + (r/r,2)2] 



213/92 ' 



(5.12) 



We refer to this equation as the "V06 density model." We note that the a component 



in the first term of Eq. 5.12 was introduced by Pratt & Arnaud (2002) to fit the inner 



slope of cuspy cluster density profiles. For a = 0, the first half of the first term 
is a /3-model. The second part of the first term - [1 + (^^/t^)'^]"^'''^ - accounts for 
any steepening around t^qq observed in these clusters. The additive second term - 
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^e02 [1 + (^Ac2 



,21-3/32 



- is simply another, smaller /3-model component, which is also 



present explicitly to fit the cluster core (and is degenerate with the a component of 
the first term). Of the nine free parameters in the V06 density model, four are present 
explicitly to fit the inner regions of the cluster (a, /32, ne02, & "f'c'i) 



V06 coupled Eq. 5.12 with a new radial temperature profile, which is 



T3D(r) = To 



(r/rcooi)"-°' + Tn,in/To 

(?^Acool)''"°°' + 1 



(r/rt)" 



:i + {rlrtYYl\ 



(5.13) 



This profile has eight free parameters. The term in the first set of square brackets 
describes any cool core present in the cluster, where the temperature approaches T^in 
as r — 0. The term in the second set of brackets describes the decline in temperature 
in the cluster outskirts, as r — )■ oo. 



Figure 5.2 compares the density profiles of the 11 real clusters presented in V06 to 



those fit to mock X-ray observations of the Kravtsov et al. (2005) sample of clusters; 



those clusters were simulated separately with adiabatic ('non-radiative') and with 
gas cooling and stellar formation ('cooling+SF' or CSF) feedback physics. The figure 
and analysis of the mock X-ray observations are presented in N07 (I simply use their 
results here). Note that outside the core (r > 0.15 rsoo), the density profiles of these 
relaxed clusters have slopes that are more similar than within r < 0.15 rsoo- This can 
be seen as the gas density profiles of the real clusters with Tx > 5 keV (magenta) and 
the CSF simulated clusters tend toward roughly the same slopes beyond r > 0.15 rsoo- 



Figure |5.3| compares the radially-averaged temperature profiles fit to Chandra X- 
ray spectroscopic data of these same 11 relaxed clusters from V06. It also shows the 



average temperature profiles of the Kravtsov et al. (2005) simulated clusters. Note 



in Fig. 5^ that the large spike in the temperature profiles of the simulated clusters, 
accompanied by a sudden drop for the CSF simulations, is generally accepted as evi- 

^V06 used this profile with £ = 3, so it is not a free parameter. 
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Figure 5.3 The above figure shows fits to the temperature profiles of 11 nearby, relaxed 
real clusters (in green, magenta, and cyan), as well as to the average temperature 
profiles of 16 clusters simulated using adiabatic (red) and cooling + star-formation 
(blue) physics. All profiles are scaled to their best-fit t^qq values. V06 used the eight- 



parameter temperature model (Eq. 5.13) to capture the details of the real cluster 
temperature profiles. Figure from N07. 

dence that one or more cluster feedback mechanisms is missing from the simulations. 
The CSF simulations in particular, which exhibit an exaggerated cusp in density 



(Fig. 5.2), suffer from this "over-cooling problem." 



5.3.2 Generalized NFW Model for ICM Pressure 



Applying the ideal gas law (Eq. 5.10) to the density and temperature fits presented in 
V06, N07 derived the ICM pressure profiles of the V06 relaxed clusters. Additionally, 
N07 analyzed the 3-D, radially-averaged pressure profiles of their simulated cluster 
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Figure 5.4 The above figure shows the X-ray-derived pressure profiles of 5 nearby, 
relaxed real clusters with Tx > 5 keV, as well as to the average 3-D, radially-averaged 
pressure profiles of 16 clusters simulated using adiabatic (red) and cooling + star- 
formation (blue) physics. The black lines are the best fit generalized NFW profiles 
for each type of cluster plotted. Figure from N07. 



sample (the Kravtsov et al. (2005) sample). N07 found that the pressure profiles of 



both simulated and real clusters can be accurately described by similar generalizations 



of the Navarro, Frenk, and White (Navarro et al. 1997, NFW) model used to fit dark 



matter halos of simulated clusters. Note that the cuspy core was excluded from the 
analysis in N07, so the over-cooling in the simulations does not affect the conclusions 
they draw about the derived pressure profile. 

The generalized NFW profile (abbreviated "genNFW" for convenience) applied 
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to the electron pressure has the form: 

Pe{r) = ^ jj-^. (5.14) 

where Pe,i is a scalar normalization of the pressure profile j^rp is a scale radius (typ- 
ically Vp ^ r5oo/l-8), and the parameters (a, 6, c) respectively describe the slopes at 
intermediate (r ^ Vp), outer (r > Vp), and inner (r ^ Vp) radii; we refer to (a, b, c) as 
the genNFW "slope parameters." Note that choosing a pressure profile similar to the 
dark matter halo's density profile is reasonable because the gas pressure distribution 
is primarily determined by the gravitationally dominant dark matter component. As 



we show in ^5.6.1, the integral of thermal pressure (ergs/cm ) can be directly re- 
lated to the (thermal) energy content of a cluster. As the dominant form of kinetic 
energy in the ICM, thermal energy tracks the underlying gravitational potential - 
and ultimately the dark matter halo - to the extent that the cluster is in pressure 
equilibrium (i.e. pressure balances the gravitational pull). It is therefore unsurprising 
that the pressure profiles of relaxed clusters should be self-similar and have a profile 
resembling that of dark matter. 

It is worth noting that, to derive pressure from the V06 profiles, seventeen pa- 
rameters (many degenerate) - eight for T(r) and nine for ne(r) - have been fit, while 
the X-ray derived pressures of all clusters in the V06 sample were well-fit by the five 
parameters of the genNFW profile (on the percent level). N07 provide best-fit, fixed 
slopes (a, b, c) for the genNFW profile, which they obtained by fitting both types of 
simulated clusters as well as the real clusters with spectroscopic temperatures above 
5 keV, presented in V06. 

The implications of this self-similarity in ICM pressure profile can be understood 
further by the analysis of cluster dynamical timescales. Since inhomogeneities in 

^Pe,i is not a central pressure, as this profile does not become flat in the core, unless c = 0. 
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pressure propagate as sound waves in the ICM, the sound-crossing time tgc is relevant 



in determining how rapidly a disturbance in pressure equilibrates. Ettori (2001) 
provides the sound-crossing time: 

Here R is the scale of a feature out of pressure equilibrium, and fcTgas is the feature's 
temperature in keV. As an example, a hot subcluster with a temperature of 10 keV 
and size of 0.1 Mpc merging subsonically with a larger cluster would reach pressure 
equilibrium with its immediate surroundings in only ~ 0.12 Gyrs. This is much less 
than the age of a typical cluster. For such a merging cluster not to be near pressure 
equilibrium, the merger would have to be very recent, or the subclump would have 
to be traveling close to the speed of sound. 

Some Common Simplifications of the Generalized NFW Model 

For different choices of slopes, the genNFW profile reduces to several interesting or 
familiar cases: 

• Dark Matter Halo Profiles : The slopes (a, 6, c) = (1, 3, 1) yield the NFW profile 



for dark matter halos. Similarly, the genNFW profile reduces to the Moore 



et al. (1999) profile by setting (a, 6, c) = (1.5,3,1), and to the Jing & Suto 



(2000) profile for galaxy clusters when (a, 6, c) = (1.1, 3, 1) 



The /3-model:: Setting the slopes to (a, 6, c) = (2, 3/3, 0) reproduces the /3-model, 
where h = 3(3. In this case P^^i = Pe,o and = r^. 

The N07 Pressure Profile : The slopes (a, 6, c) = (1.3,4.3,0.7) yield the profile 



we use in fitting real clusters of galaxies (see the discussion below, in ^5.3.3) 



The N07 Pressure Profile for Simulated Clusters : A slightly cuspier inner slope. 
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c = 1.1, yields the best-fit pressure profile clusters simulated with CSF physics. 
The full set of slopes is therefore (a, b, c) = (1.3, 4.3, 1.1). 



Slope of the Generalized NFW Pressure Profile (dP/dr) 

We provide here the derivative, with respect to r, of the genNFW pressure profile. 
This will be useful when estimating the cluster total mass, Mtot, using the equation 



of hydrostatic equilibrium (Eq. 5.32, ^5.5.4) 



dP{r) 
dr 



p . 



1 + 



-{a+b—c)/a 



+ C 



(5.16) 



5.3.3 The NOT Pressure Profile 

In order to reduce the number of free parameters to those that can be constrained 
by existing SZE observations of clusters, we fix the slopes of the genNFW profile to 
the best-fit slopes provided by NOT, and test their ability to extract cluster physical 
parameters in Chapter[6j By using these fixed slopes to fit the SZE, we take advantage 
of the self-similarity of cluster pressure profiles across a broad range of masses. N07 
found that the pressure profiles of all sixteen clusters simulated using cooling and star 
formation physics could be well fit with the slopes of the profile fixed to the same 
values (namely (a, b, c) = (1.3, 4.3, 1.1)), while the pressure profiles of several relaxed, 
real clusters, studied in detail using Chandra and presented in V06, could all be fit 
using a slightly different value for the inner slope, c = 0.7. Therefore, the pressure 



profile we test - Eq. 5.14 with slopes fixed at (a, 6, c) = (1.3, 4.3, 0.7) - becomes 



PAr) 



p . 



(5.17) 



The clusters simulated with non-radiative, adiabatic physics required quite different 
slopes, which we do not consider here, due to the simplified physics used in their 
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simulation. As demonstrated by N07, we note that the same outer slope used to 
fit the V06 (real) cluster and the CSF simulated clusters, b = 4.3, also fit well the 
clusters simulated using adiabatic physics. 

5.3.4 Simplified Vikhlinin Density Model — An X-ray Density 
Model to Complement the NOT Pressure Profile 

The density model we chose for fitting X-ray imaging data is a simplified, core-cut 
form of the V06 density model. We justify this core cut because we are primarily 
interested in recovering cluster parameters at r25oo and rsoo, where the contribution 
from the core is negligible. We refer to our simplification as the "Simplified Vikhlinin 



Model" or "SVM" hereafter. Many have demonstrated (see, e.g. LaRoque et al. 



2006 ) that a core-cut isothermal /3-model recovers Mgas accurately outside the inner 
100 kpc, out to r25oo- We chose this simple modification as a suitable description of 
intermediate- quality X-ray data outside the core, out to, and possibly beyond, rsoo- 
We note that V06 presented this profile as an intermediate step to obtaining the full 



9-parameter density model (Eq. |5.12 ), but never used the SVM to fit cluster data. 

We begin with the standard /3-model. Following V06, we add a scale radius at 
which the density starts to decline more rapidly, with a slope e. The SVM is therefore 

p(r) = po[l + (r/r,)2]-'^/' [l + (r/r,)3]-/^ (5. 18) 

The two extra degrees of freedom, in and e, attempt to overcome the limitations 



of the /3-model in fitting a cluster out to rsoo (for more detail, see Mroczkowski et al. 



2008D . Several authors (see e.g |Neumannl|2006| [Vikhlinin et al.||2006t |Maughan||2007[ ) 



have noted the /3-model is insufficient for describing clusters at r > rsoo, where their 
density profiles of steepen. Thus we have developed the N07+SVM pair of models - 



Eqs. 5.17 and 5.18 - primarily to fit clusters used in determing the SZE+X-ray bulk 
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Figure 5.5 This figure shows the density of the standard /3-model (red) and SVM (blue) 
using the steepest outer slope V06 allowed, e = 5. This was chosen to illustrate the 
largest disparity the radially- average outer cluster slope V06 considered "realistic." 
We chose typical parameters for the components common to the two profiles (namely, 
we chose = 100" and (3 = 0.7, and set the normalization po = I for simplicity). For 
the SVM, we show a typical massive cluster's scale radius for the slope to steepen, 
r, = 500". 



cluster scaling relations. 



The SVM simplifies to a /3-model when e = 0. Figure 5.5 shows the standard 
/9-model (red) and SVM (blue) using e = 5, the steepest outer slope V06 allowed in 
their density fits. Here, we chose typical parameters, found in cluster fits, for the 
components common to the two profiles (namely, we chose = 100" and /3 = 0.7, 
and set the normalization po = 1 for simplicity). For the SVM, we chose a scale radius 
of Tg = 500", corresponding to the radius at which the density of a typical cluster at 
low redshift begins to steepen. Note that since X-ray surface brightness scales as p^^, 
the signal will be < 0.1% that of the core for the regions where the profiles diverge 
(r > 200"). 
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5.3.5 Combining the NOT and SVM Profiles 



We apply the N07 pressure profile (Eq. 5.17) to fit the SZA data presented in Chap- 



ter combining it with the SVM (Eq. 5.18) to fit complementary Chandra X-ray 



imaging data. As discussed in ^1.4[ fitting X-ray surface brightness (Eq. 1.8) re- 
quires both density and temperature information. Having assumed the ideal gas law 



(Eq. 5.10), the electron temperature profile Te(r) occupies a unique position while 



jointly fitting the N07+SVM profiles to X-ray+SZE cluster observations: it can be 
derived from the pressure and density fits without relying on X-ray spectroscopy. For 
each iteration of the MCMC fitting routine, the temperature profile necessary to fit 
the X-ray surface brightness is derived for each set of test parameters of the density 
and pressure profiles. The weak temperature dependence of the X-ray imaging data 



(discussed in ^1.4) has the ability to exclude some of the pressure model's parameters 



that fit well the SZE data alone, as will be shown in Chapter [6] 



5.4 Markov chain Monte Carlo Analysis 

A Markov chain Monte Carlo (MCMC) method, implemented in a program called 
markov, is used to determine the distribution of models that fit the data, as described 



by Bonamente et al. (2004); LaRoque et al. (2006); Bonamente et al. (2006 2008). 



Iterations are accepted in the MCMC according to the Metropolis-Hastings algorithm, 
which ensures that a fully-converged MCMC reflects the true probability distribution 
of the resulting model fits. Assuming the model itself as a prior, the output chain of 
parameters - the "Markov chain" itself - can be used to determine the probability 
density of the model fits. 



The SVM has five free parameters to describe the gas density (see Eq. 5.18) 



and the N07 profile has two free parameters to describe the electron pressure (see 
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Eq. 5.17). Additional parameters such as the cluster centroidj^the X-ray background 
level, and the positions, fluxes, and spectral indices of field radio point sources are 
also included as necessary. 

The model fitting code was originally called jointfit, and used a downhill simplex 



method (see, for example, Bonamente et al. 2004 Joy et al. 2001 Grego et al. 2000 



200Tj [LaRoque et"aL]|2003| |Patel et al.||2000t [Reese et al.|[2000l [2002| and references 
therein), markov and jointfit have been used extensively to fit Chandra and ROSAT 
X-ray observations jointly with SZE observations made by OVRO and BIMA. Because 
the UVFITS files produced from SZA data contain sixteen bands, as opposed to the 
single-banded UVFITS files from OVRO/BIMAj^ mar/coi^'s UVFITS reader required 
an upgrade in order to read and properly account for the frequency of each band. 
Furthermore, in the transition from 32-bit to 64-bit desktop processors, the standard 
CFITSIO libraries did not work, and I upgraded the reader to the well-documented 



UVFITS reader Martin Shepherd wrote for Difmap (Shepherd 1997). 



For each iteration in the MCMC process, model test parameters are generated by 
stepping randomly from the last accepted link in the Markov chain. The size of the 
step is determined using a top hat distribution of adjustable width. These test param- 
eters are used to compute the model image over a regular grid, with sampling defined 
to be less than half the smallest scale the SZA can probe (i.e. Nyquist sampling of 
the scale probed by the longest baseline in the observation). This image is multiplied 
by the primary beam of the SZA, transformed via FFT to Fourier space (which is 
where the data are directly sampled by a radio interferometer), and interpolated to 
the Fourier-space coordinates of the SZE data. The likelihood calculation for the SZE 
data is then performed directly in Fourier space, where the noise properties of the 



^Since they are separate observations, and the coordinates used in fitting are simply in offset 
from the pointing center of the observations, the SZE and X-ray centroids are fit independently. 
They are, however, later checked for consistency. 

^OVRO had two bands, but each band was output as a separate UVFITS file. 
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interferometric data are well-characterized. 

The MCMC routine determines which trial sets of parameters are accepted into 
the Markov chain by computing the joint likelihood £ of the models' fits to the X-ray 
and SZE data. The SZE likelihood is given by 



where Ai?j and A/j are the difference between model and data for the real and 
imaginary components at each point i in the Fourier plane, and Wi is a measure 
of the Gaussian noise (1/cr^). Since the X-ray counts, treated in image space, are 
distributed according to Poisson statistics, the likelihood of the surface brightness 
model fit is given by 



where Mj is the model prediction (including cluster and background components), 
and Di is the number of counts detected in pixel i. Spatial regions of the X-ray 
image, containing for example the cluster core, gas clumps, or X-ray point sources, 
can be excluded from the X-ray analysis simply by excluding data in those regions 



The joint likelihood of the X-ray and SZE models is given hy C = Csze ■ ^Xvay- For 
the models that include X-ray spectroscopic temperature constraints, the likelihood 
is £ = £szE • £xray " ^xsPEC, whcrc £xsPEC comes from a table of values computed 
using a plasma emissivity model fit to the X-ray spectroscopic data, using the software 

^'^An alternative approach might involve, for example, interpolation past the region of the X-ray 
image that we wish to exclude. We choose not to make any such assumptions about the cluster in 
the excluded regions. 




(5.19) 




(5.20) 



from the X-ray likelihood calculation (Eq. 5.20[ ) p°| 
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package XSPEC (Arnaud 1996). The relationship between and £xspec is 



-CxsPEC = exp 



X 



(5.2i: 



Throughout a markov run, the probabihty distribution of fits to the data is ex- 
plored. Some number of iterations at the beginning of the Markov chain, however, 
will depend on the starting parameters markov chooses; We refer to this as the "burn- 
in period" - when the MCMC still depends on the starting conditions - and ignore 
iterations accepted during this period. After a sufficiently large number of total it- 
erations, the allowed range of parameters is fully explored. The number of required 
iterations depends on the number of free parameters, the size of the top-hat distri- 
bution out of which the test parameters are chosen, and the noise in the data. This 
is why we use convergence tests, described below, to determine the total number of 
required iterations. 

The likeliest parameters will appear multiple times in the Markov chain. Their 
distribution after the initial burn-in period is referred to as the "stationary distribu- 
tion," since further iterations do not yield better constraints on the model. We follow 



the methods outlined in Bonamente et al. (2004), first exploring large regions of pa- 



rameter space with test runs to ensure the solution is not stuck in a local minimum. 
Further runs are used to narrow in on the likeliest parameters by limiting the param- 



eter range to explore and tuning the top hat distribution width. As in Bonamente 



et al. (2004), we use the Raftery-Lewis test to determine convergence to a stationary 



distribution ( [Raftery fc Lewis||1992[ [Best et al.||1995| |Gilks et al.||1996t [Plummer et"aL 



2006), This is an efficient way to determine how many iterations the initial burn-in 



should contain, and how many further iterations are required to reach the stationary 
distribution. Once the MCMC routine has found the stationary distribution, we use 
the set of accepted model parameters to derive cluster parameters such as Mtot and 
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their error bars. 



5.5 Mass analysis 



5.5.1 Weighting factors for the cluster gas 

The first quantity we consider is the cluster gas density. X-ray spectroscopic data 
have the abihty to constrain the metalhcity Z of a gas (the mass fraction of elements 
heavier than helium) by fixing the abundance ratios to those observed in the solar 
atmosphere and fitting for the few prominent emission lines -notably those from iron 
- that can be found in X-ray spectra. Since the lighter elements (e.g. hydrogen, 
helium) are fully collisionally-ionized at cluster temperatures (> IkeV), no emission 
lines are observed from these, and their abundances are assumed to be those produced 
by Big Bang Nucleosynthesis (a.k.a. "cosmic abundance"), plus the amount produced 
by stellar nucleosynthesis, for the X-ray spectral fit to the metals. 



When comparing the SZE-derived temperature (see § 5.3.5) to the X-ray spec- 



troscopic temperature, we use the program XSPEC (Arnaud 1996) to fit a plasma 



emissivity model to the X-ray data. The plasma spectroscopic emission models we 



most commonly fit to these data are the Raymond-Smith (Raymond & Smith 1977) 



and MEKAL ( |Mewe et al^|1985t [Kaastra fc Mewe||1993| [Liedahl et aL]|1995D models, 
which can have free parameters of temperature, redshift, and metallicity (redshift is 
often fixed to that found using complementary optical observations, though typical 
X-ray data can constrain a cluster's redshift to better than ±1%). Fitting simulated 
cluster emission with X-ray spectra generated using the MEKAL model, we found 
that the metallicity recovered in XSPEC from fits of the Raymond-Smith model was 
within > 80% accuracy of the metallicity in the input MEKAL spectrum. 

Using the fit metallicity, we tested the X-ray surface brightness fitting procedure 
- which relies on the Raymond-Smith model when fitting X-ray emissivity - against 
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Element 


Mass Fraction 


Number of Species 


Mass 


H 




le" + lp+ 


IrUp 


He 


Y = 4nHe/ (riH + 4raHc) 


2e- + 1^2+ 


Arrip 



Table 5.1 Common parameters used to compute the weighting factors for a pure H 
and He plasma. 

mock X-ray observations of these simulated clusters. While these mock X-ray expo- 
sures were generated using the MEKAL model, the emissivity Aee(Te, Z) used to fit 



the mock observations (in Eq. 1.8) was computed using the Raymond-Smith model. 
The recovered gas density was accurate to better than 98%. We conclude that the 
continuum emission in an X-ray image is insensitive to the plasma emissivity model 
assumed in the fit, and continue to use the Rajmaond- Smith model to fit X-ray surface 
brightness. 

After determining the metallicity of the ICM (for either real or mock observations), 
several useful "weighting factors" can be determined by carefully accounting for the 
species of which the gas is composed. The weight per particle is determined by 
summing the masses of all the species (i.e. the electrons, protons, and atomic nuclei) 
in the plasma and dividing by one proton mass per particle. 

_ riH X Imp + riHe x Auip 

^ ~ [rinile- + lp+) + nueila^+ + 2e-)]mp' ^ ' ^ 

For example, a fully-ionized plasma consisting of 10% helium and 90% hydrogen by 
number yields 

0.9 + 0.1x4 

^ " 0.9(le- + lp+) + 0.1(la2+ + 2e-) " ^^'^^^ 

where we use the nomenclature e~ for an electron, p+ for a proton, and for a 
helium nucleus, and we have assumed a hydrogen atom is exactly one proton mass 



nip and a helium atom is Anip (see Table 5.1). 



Similarly, we can define fie, the weight per electron, by summing the masses of 
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the species and dividing by Imp for each electron. 



nn X Irrip + nne x 4mp 
[nu{le-) + raHc(2e-)] nip 



(5.24) 



For the example in Eq. 5.23, this means 



fJ-e 



0.9 + 0.1 X 4 
0.9(le-) + 0.1(2e- 



1.1^ 



(5.25) 



In practice, the detailed calculation of weighting factors requires accounting for all 
the species in the gas. Having fit spectroscopically a metallicity Z, we use the abun- 



dance calculation of Anders & Grevesse ( 1989 ) to determine the relevant weighting 



factors. We assume this abundance model for all results presented in this thesis 
(though we consider the effects of a non- universal, non-constant Y/X ratio in 



5.5.2 Relating the fit ne(r) and Pe{r) to ICM gas properties 

The quantities fie can now be used to relate the electron number density profile ne(r) 
to the gas (mass) density Pgas('"), which is 



r) = HempUeir). 



(5.26) 



Similarly, the number density of all species in the gas n{r) relates to the gas density 
through fi: 

(r) = ixmpTiir). (5-27) 



Assuming the electron temperature equals the gas temperature (Tg = Tgas)p] we 



combine Eqs. 5.26 & 5.27 and use the ideal gas law to find the electron pressure 



^^This is because ^ and /ie are mass weighting factors, so we have to relate pressure to density to 
take advantage of them. Hence we have /i//ie — n(.{r)/n{r) — Pe{r) / Pgi^sif) ■ 
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profile relates to the overall gas pressure via 



^gas(r) = ( ^ ) Pe(r). (5.28) 



5.5.3 Mass of the X-ray-Emitting ICM Gas (Mg,,) 

For a spherically-symmetric density profile, we now simply integrate Pgas('") to obtain 
the gas mass Mgas within a spherical volume with radius r: 

M^^{r) = 47r / p^^ir'y'^dr' . (5.29) 

JO 

5.5.4 The Total Mass of the Cluster (Mt„t) 

The total gravitational mass within a spherically-symmetric system in pressure equi- 
librium can be found by considering the balance between the inward gravitational 
force fin and the outward pressure font- This scenario of pressure support and neg- 
ligible net mass transport across the surface of the sphere is known as hydrostatic 
equilibrium (HSE). 

Consider the forces on a spherical shell of thickness dr, at radius r, with gas mass 
m — Airr^ Pgas('") dr. The inward force of gravity is 

_ GmMtotir) _ G Anr^ pg^,{r) Mtot{r)dr 

where G is Newton's gravitational constant and Mtot(r) is the total gravitating mass 
within radius r. The outward force of pressure is 

fout = 47rr2 dr (5.31) 
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Setting fin = font and solving for Mtot(r) yields 



Notice that this is a powerful way of using the observable gas distribution in a 
cluster of galaxies to estimate its total matter content (i.e. dark + baryonic). By 
assuming further the pressure is entirely therma ^ and that the ideal gas law holds, 
one only has to determine two of the three quantities {Pga.sij), Pga.s{r), and T{r)} in 
order to estimate a cluster's total mass. 

5.5.5 Hot Gas Mass Fraction (/gas) 

For each accepted iteration in the MCMC, we have a set of model parameters fit to 
the available SZE+X-ray data, from which we can compute Mgas and Mtot- Solving 



for Mgas and Mtot using Eqs. 5.29 & 5.32, we simply solve for the profile of the gas 



mass fraction individually for each set of MCMC parameters: 

(/gas(r)) = (Mgas(r)/Mtot(r)). (5.33) 

Doing this, the mean gas fractions we present in this thesis take advantage of the fact 
that MCMC explores the probability density of model fits to the data. Furthermore, 
the confidence intervals presented also directly refiect the range of models fit (rather 
than assuming the derived Mgas and Mtot are statistically independent). 

Turbulence and magnetic fields are possible sources of additional, non-thermal pressure support. 
The standard HSE total mass estimate could, for example, be systematically low by ~5-20% because 
it does not account for subsonic motions of the ICM, which can contribute to the total ICM pressure 



but will not affect the temperature and density (see e.g. Nagai (20061 for details). Additionally, 
there are indications that cosmic rays could contribute to radiative support (see e.g. Pfrommer et al. 
pOOTl)). 
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5.5.6 Overdensity Radius 

In order to compute global cluster properties such as Mgas, Mtot, and /gas, one needs to 
define a radius out to which all quantities will be calculated. Following |LaRo que et al. 



(2006); Bonamente et al. (2006), we compute global properties of clusters enclosed 
within the overdensity radius ta, within which the average density of the cluster is a 
specified fraction A of the critical density, via 



^7rp,{z)Arl = Mtot(rA), (5.34) 



where Pc{z) is the critical density of the Universe at redshift z, and is computed (with 
respect to Pc = SHq/StcG, the critical density at redshift z = 0): 

p,iz) = p,[nMii + zf + nA]. (5.35) 

Throughout this work, we evaluate cluster properties at density contrasts of 
A = 2500 and A = 500, corresponding to average densities of 2500 and 500 times the 
critical density at the redshift of the cluster. The overdensity radius r25oo has often 
been used in previous SZE+X-ray studies, since that was the maximum radius attain- 



able in many X-ray observations of intermediate redshift clusters (see, e.g. LaRoque 



et al. (2006)). The overdensity rsoo is now reachable in deep SZA and Chandra X-ray 
data, without extrapolating the model beyond the image fitting region or into regions 
with signal-to-noise ratio S/N < 1. 



5.6 SZE-Specific Quantities 
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5.6.1 The Integrated, Intrinsic Compton y Parameter {Y\nt) 

Since the SZE signal arises via Compton scattering, the integrated SZ effect, when 
scaled to the distance of the cluster, can be converted to Yint, the intrinsic, integrated 



Compton y parameter (see e.g. da Silva et al. (2004); [Nagai (2006)). The intrinsic 



line-of-sight (abbreviated "los" in the subscripts to follow) integrated Compton y 



parameter Yi^s is obtained by integrating y (Eq. 1.2) over the solid angle Vt subtended 



by the cluster, and scaling the area for the angular diameter distance of the cluster, 

>ios = D\ [ ydn = [ de [ P,dA. (5.36) 

Here A is the area of the cluster in the plane of the sky. 

In addition, SZE observations constrain the intrinsic volumetric (abbreviated "vol" 
in the subscripts to follow) integrated Compton y parameter Y^gi, obtained by inte- 
grating y within a spherical volume of radius r: 

noi = ^ fPe{r')r'Hr. (5.37) 



rrieC^ Jo 

Since l^nt ^ in the form of either l^oi and llos^ scales as pressure (ergs/cm^) 
integrated over volume (cm^), it is proportional to the thermal energy content of the 
ICM within that volume. As the dominant form of kinetic energy in the ICM, thermal 
energy tracks the underlying gravitational potential, and ultimately the dark matter 
halo of the cluster, to the extent that HSE holds. SZE flux thereby provides a robust, 
low scatter proxy, YJut, for total cluster mass, Mtot- 
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5.6.2 The Total Thermal Energy Content: 
An SZE-only Scaling Quantity 

A problem becomes apparent when one considers how the N07 pressure profile - or 
any profile fit to SZE data alone - could be used in an SZE survey to constrain 
cosmology. To maximize the utility of a cluster survey, it is crucial to understand 
how to relate a cluster's observable properties to its total mass. 

When using joint SZE+X-ray data, one has the luxury of solving for Mtot by 



assuming HSE (Eq. 5.32). Joint SZE+X-ray fits therefore allow one to solve for ta, 
the radius that encloses matter overdensity A. We can then compute Mtot('"A) and 
^intl'^A) (e.g. either yios('"A) or I^oK'^a), provided in Eqs. 5.36 & 5.37[ ), which, for 



example, could be used to establish the Yi^t{rA)-Mtot{rA) scaling relation]^ 

For an SZE-detected cluster with no complementary X-ray observation (or some 
independent estimate of Mtot, such as from lensing or galaxy velocity dispersion), the 
problem arises that we do not have a way to determine r a , and therefore do not know 
within which radius to compute Yi^t{r/^). There is no direct way to find from an 



SZE-only observation of a cluster, since we do not know Mtot{f) (see Eq. 5.34). Unlike 
the isothermal /3-model, the N07 pressure profile can be integrated over an infinite 
volume, yielding a finite estimate for the total thermal energy content of the cluster. 
The thermal energy content should scale directly with the total mass of a cluster - to 
the extent that thermal pressure supports the cluster against gravitational collapse - 



and has been studied before by others (e.g. Afshordi et al. 2007). The total thermal 



energy content i^tot of a cluster that can be described by the spherically-symmetric 



^■^One could also use the derived ta to obtain Aftot(?'A) from a previously-established scaling 
relation, though if one has ta, one necessarily has MtotC^'A) already. This is simply because 
Aftot(fA) = ^PcT^, where pc is the critical density of the Universe at that redshift. 
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pressure profile -Pgas('") is 



St, 



47rPgas (r)r^(ir. 



(5.38) 



Using Eq. |5.28 where the SZE-constrained electron pressure profile is parameter- 



ized in the form of the generalized NEW pressure profile (Eq. |5.14), a tool such as 



Mathematical^ can handily solve this integral^ Eor a > 0, 6 > 3, and c < 3, the 
integral converges, yielding 



-£'tot,genNFW — SvTr Pg, 



3 o / /^e 



/i 



r 



'3-c\ 



aVi 



(5.39) 



Here F is the Euler gamma function. For fixed slopes (a,6, c), the integral is only 
dependent upon profile parameters Pe,j and r^, and scales as Y^oi- The parameters 
recommended by Daisuke Nagai (personal communication) for his simulated clusters 
all had a G [0.5,2], c G [0.5, 1.5], and b G [4,5], which covers a range consistent with 
a finite thermal energy content. The N07 profile provides (a,6, c) = (1.3,4.3,0.7)^^ 
which also yields finite thermal energy. In Figure 5.6[ we show the NO 7 profile predicts 
a thermal energy content within a given radius that converges to a finite total value. 



The fact that the integral for the ICM thermal energy (Eq. 5.38) converges implies 
there is a single, physical quantity that can be constrained by SZE data alone, without 
having to know ta (which would otherwise have to be determined using independent, 
complementary observations). The total thermal energy is simply l^oK^ — >■ oo) ex- 



pressed in different units (c.f. Eqs. 5.37 and 5.38), and is not - by any means - a 
new concept. However, one could envision measuring how Etot scales with Mtot(?"5oo) 



//www. wolf ram. com 



-"^^Physically, this integral should not be carried out to infinity. However, the generalized NFW 
profile drops rapidly enough that any contribution from outside the virial radius is negligible. 

These are the original values for (a, 6, c), published in N07. These were updated after the 
completion of this work, however, and will be published in an erratum to N07. Fits utilizing the 



newer values are presented in Mroczkowski et al. (2008). 
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Ratio of Thermal Energy to Total 




Figure 5.6 The ratio of the thermal energy content within a given radius to the total 
as r — 7- oo, plotted for the fit to A1835 of the N07 pressure profile. The vertical line 
shows rsoo — 360" ~ 1.4 Mpc. 



for a small number of clusters that have high quality SZE and X-ray or lensing ob- 
servations. Then, as a larger SZE survey (such as those planned for ACT and SPT) 
discovers new clusters through the SZE, this scaling relation could be directly applied 
to obtain an Mtot('"5oo) estimate for each cluster in that sample, which is a useful step 



in constraining cosmology using clusters (see ^1.1). 
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Comparison to the /3-model for i?tot 



In contrast to -E'tot.genNFW? we note that the integral in Eq. 5.38 when applied to the 
/3-model diverges for /3 < 1. For /3 > 1, the total thermal energy of the cluster gas is 



r 



/ 3/3-3 ^ 



r(f ) ■ 



(5.40) 



Since /3 ~ 0.7 for X-ray observations, the isothermal /3- model generally predicts an 
infinite thermal energy content in clusters; this is a product of assuming isothermality 
where the assumption is no longer valid. As mentioned in §5.3.2[ the generalized 
NFW profile can be simplified to the /3-model by choosing {a,b,c) = (2,3/3,0). The 
condition that /3 > 1 is required to obtain a finite total thermal energy within a cluster 
is consistent with the condition that the genNFW parameter 6> 3 (i.e. /3 < 1 =^ b<3). 



Comparison of Total Thermal Energy to an X-ray Proxy 



Recent work by Kravtsov et al. (2006) has demonstrated an X-ray quantity that, like 



the SZE quantity Yint, provides a robust proxy for Mtot- This quantity is defined Yx = 
Mgas{f50o) Here Tx is measured within an annulus in radial range [0.15, 1] rsoo- 

This range was chosen because cluster temperature profiles are most self-similar over 



this range (Vikhlinin et al. 2005 Nagai et al. 2007b). This large core cut makes the 



measured Tx insensitive to effects within cluster cores. 

Using the ideal gas law (Eq. |5.10[ ) and noting the similarity between the integral 



to obtain the gas mass Mgas (Eq. 5.29) and that to obtain the thermal energy E 



(Eq. 5.38), it is clear that Yx scales as the thermal energy content within rsoo, assum- 



ing isothermality. The implication of this is that the SZE-determined thermal energy 



^■^Note that 



Maughan 



( 2007 ) demonstrated this quantity on a sample of clusters using an approxi- 
mate formula for rsoo , so the X-ray data did not have to be of sufficient quality to perform a detailed 
HSE mass estimate. An overestimate of r^oo yields an overestimate of Mgas ('"500)1 but also lowers 
the estimate of Tx- Similarly, an underestimate of r^QQ lowers the estimate of -Mgas (''500)1 but raises 
Tx- In this way, Yx is surprisingly robust to errors, both systematic and statistical in nature. See 



Kravtsov et al. (2006) for a more detailed analysis. 
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content could provide the basis for a powerful, SZE-only observable that could read- 
ily be applied to upcoming SZE cluster surveys, without the need for corroborating 
observations other than those required for redshift determination. 
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Chapter 6 

Applications of the Models 



In this chapter, I demonstrate the apphcation of the N07+SVM (Eqs. 5.17 & 5.18) 



to observations of real clusters. I initially tested the markov fitting code and new 



models on mock SZA observations of several of the Kravtsov et al. (2005) simulated 
clusters. These mock SZA observations were also jointly fit with mock Chandra 
observations. The tests provided sufficient confidence that the modeling code was 
correctly implemented, but the mock SZA observations suffered systematic biases 



due to the unrealistically cuspy cores of these simulated clusters (see Fig. 5.3 and 



discussion in 5.3.1 ). Any attempt to excise these cluster cores, which often dominated 
the SZE signal, from the input simulation data would involve several assumptions 
about the cluster core and physics missing from the simulations. Instead I chose to 
use the mock observations to verify the functionality of the modeling code, and then 
focused on tests involving real observations of real clusters. 

I selected three massive clusters, well studied at X-ray wavelengths, spanning a 
wide range of redshifts {z = 0.17-0.89) and cluster morphologies. These clusters are 
used test the joint analysis of Chandra and SZA observations, using the models and 



methods described in Chapter pj Specifically, I test the N07 pressure profile (Eq. 5.17) 



in conjunction with the SVM density profile (Eq. 5.18). I compare the results of these 
fits with results from a detailed. X-ray-only analysis, as well as with a joint SZE+X- 
ray analysis using the traditional isothermal /3-model. I assume a ACDM cosmology 
throughout the analysis presented in this chapter, with Qm = 0.3, = 0.7, and 
Qk = 0. The angular diameter distances to each of these clusters, as well as the 



redshifts and spectroscopic temperatures of the clusters, are hsted in Table 6.1 
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Cluster 


(Mpc) 


(keV) 


A1835 

CL1226 

A1914 


0.25 806.5 
0.89 1601.8 
0.17 600.6 


^•^^-0.37 

8.63l°f 



Table 6.1 Clusters chosen for testing the models. Angular diameter distances were 
computed assuming Qm = 0.3, = 0.7, & fi^ = 0. 



" Redsh i fts for A1914 and A1835 are from |Struble fc Rood] ( [19991 ) . Redshift for CL1226 is fro m|Ebeling 
et al.| (2001). All are in agreement with XSPEC fits to iron emission lines presented in LaRoque 



et al. (2006). 



Global X-ray spectroscopic temperatures were determined in the range r e [0.15, 1.0] rsoo- These 
temperatures were used in the isothermal /3-model analysis. 




Figure 6.1 X-ray image of A1835, showing it to be relaxed. The X-ray analysis 
of A1835 relies on a single Chandra ACIS-I exposure, with 85.7 ks of good time 
(unflagged exposure time). The pixels of the ACIS-I detector are binned to be 1.968" 
on a side. The X-ray image shown here is smoothed with a Gaussian that is 2 pixels 
in width (for display purposes only). See Table 6.3 for more details on the X-ray 
observation. The inner 100 kpc (core) and all detected X-ray point sources were 
excluded from the X-ray surface brightness and spectroscopic analyses. 
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Figure 6.2 X-ray image of CL1226, showing this high redshift cluster to be approx- 
imately circular on the sky, and thus apparently relaxed. This image combines two 
Chandra exposures. The inner 100 kpc (core) and all detected X-ray point sources 
were excluded from the X-ray analysis. 

6.1 Cluster Sample 



Abell 1835 (A1835) is an intermediate-redshift {z = 0.25) cluster (see Fig. 6.1, located 
at the mean redshift of a sample of clusters being used by the SZA collaboration to 
determine the SZE flux scaling relations (versus Mgas, Mtot, Tx, etc.). This cluster 
sample is X-ray fiux-limited, and is located in the redshift range z ~ 0.2-0.3. The 
sample has both X-ray observations and optical lensing mass estimates available, and 
will appear in the thesis of Ryan Hennessey as well as forthcoming SZA papers. A1835 
is relaxed, as evidenced by both its circular morphology in the X-ray images and its 



cool core (see Fig. 6.1, and Peterson et al. (2001) for more details). To demonstrate 



the applicability of the joint SZE+X-ray analysis to high redshift clusters, I analyzed 



CL J1226. 9-1-3332 (CL1226), an apparently relaxed cluster at 2; = 0.89 (see Maughan 



et al. (2007) - hereafter referred to as M07 -and Fig. 6.2). Note that M07 argues 
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Figure 6.3 X-ray image of A1914, showing it to be disturbed. This image combines 
two Chandra exposures. A large subclump, measured to be hot by M08, is located 
to the east (left) of the cluster center. This subclump, the inner 100 kpc (core), and 
all detected X-ray point sources were excluded from the X-ray analysis. 

despite appearing spherically symmetric, that the lack of a cool core in this cluster is 
likely due to a recent merger. Finally, to assess how applicable this pressure profile 
is to disturbed clusters, I also analyzed Abell 1914 (A1914), an intermediate redshift 
{z = 0.171) cluster with a hot subclump near the core (see Fig. 6.3); when the 



subclump is not excluded from the X-ray analysis, a large centroid shift is present 



in the fit to the surface brightness (see Maughan et al. (2008) - hereafter referred 



to as M08 - who use the size of the error in the fit X-ray centroid as a measure 
of cluster dynamical state). I exclude this subclump from the X-ray analysis of 
A1914 by defining a ~ 100 kpc region centered on the peak brightness of the clump 



and excluding this region from the X-ray likelihood calculation (see Eq. 5.20 and 
discussion of the X-ray likelihood in ^5.4[ ). 



Figures 6.4-6.7 show the c/ean-ec^ interferometric maps, produced from the SZA 



^ clean-ing interferometric data involves iteratively fitting point sources within a specified region 



and removing their flux. See e.g., Thompson et al. (2001) for more details 
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Cluster 


SZA Pointing Center (J2000) 


Short Baselines (0.34-1.5 kA) 


Long Baselines (3-7.5 kA) 


^int .red 




a 


S 


beam(" X 




beam(" X ")^ 


o-(mjy)': 


(hrs) 


A1914 




+37°49'35".7 


117.5x129.9 


0.30 


23.5x17.4 


0.35 


11.5 


A1835 


14''01'"02''.0 


+02°52'41".7 


116.6x152.1 


0.25 


17.5x23.5 


0.33 


18.6 


CL1226 (30 GHz) 


12''26'"58''.0 


+33°32'45".0 


117.4x125.4 


0.20 


16.0x21.2 


0.20 


22.0 


CL1226 (90 GHz) 






42.3x39. l<i 


0.42 


9.74x7.46=* 


0.32 


29.2 



Table 6.2 SZA Cluster Observations 



"Unflagged data after reduction. 

^Synthesized beam FWHM and position angle measured from North through East 
■^Achieved rms noise in corresponding maps 

''The short and long baselines of the 90-GHz observation probe 1-4.5 and 8-22 kA. 



Cluster 


Joint SZE+X-ray Analysis 


Maughan X-ray Analysis 


ObslD 




(ks)*^ 


// b 


(ks)^ 


// b 




A1914 


26.0 


34.4-423.1 


23.3 


0.0-462.5 


542+3593 


A1835 


85.7 


25.6-344.4 


85.7 


0.0-519.6 


6880 


CL1226 


64.4 


12.9-125.0 


50 


0.0-125.0 


3180+5014 


CL1226" 


N/A 


N/A 


75+68 


17.1-115 


0200340101 



Table 6.3 Details of X-ray Observations. The X-ray analysis presented here, as part 
of the joint SZE+X-ray modeling, was performed independently from that performed 
by Maughan. 



"Good (unflagged, cleaned) times for X-ray observations after respective calibration pipelines. 
''X-ray image fitting region. 

'^XMM-Newton observation of CL1226, presented in M07, was used only in the spectroscopic analysis. 
The exposure times are, respectively, those of the MOS and PN camera. 

observations analyzed here, of these three clusters. Details of the SZA and Chandra 
observations, including the X-ray fitting regions, the unflagged, on-source integration 
times, and the pointing centers used for the SZE observations, are presented in Tables 



6.21 and 6.3 Note that 90-GHz SZA data were included in the fits to CL1226, which 



is only partially resolved at 30 GHz due to its small angular extent (see ^6.4). 

In the following sections, I discuss the results of the joint analysis of the SZE and 
X-ray data. I show how the models fit the X-ray surface brightness and radio inter- 
ferometric SZE data, and use these fits to place constraints on cluster astrophysical 
properties. Finally, I compare these results with an independent, detailed. X-ray-only 
analysis performed by Ben Maughan, following the techniques presented in M07 and 
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M08 and described in §6.3[ The high quahty of the X-ray data on these three clusters 
allows this X-ray-only analysis, which is used in assessing the new models. 

6.2 Unresolved Radio Sources 



As discussed in ^ 4.3.2, an unresolved, compact radio source ("point source") can be 
represented mathematically by a 2-dimensional Dirac delta function in image space. 
In Fourier space, such unresolved sources - which can contaminate an SZE cluster 
observation - have an analytic transform: a constant, equal flux on all scales, as 
measured by all baselines. The flux from these sources is modeled analytically using 
a power law to describe the frequency-dependence of the source flux, with flux /o 
normalized at uq = 30.938 GHz, the center of the SZA 30-GHz band. I also account 
for attenuation by the primary beam, A^{x,y), at each band's frequency. Therefore, 
the measured flux /(z/) from a point source at map location {x, y) is fit using 

fiiy) = A,{x,y)foiiy/uor (6.1) 

where a is the spectral index. 

When modeling a compact source detected in an SZE observation, I fix its spec- 
tral index a so that it simultaneously fits the mean flux measured by the SZA at 



30.938 GHz and that measured in the NVSS ( [Condon et al.||1998[ ) or FIRST ( [White 
et al. 1997) surveys, which were performed at 1.4 GHzj^ The mean point source flux 
and the source's approximate coordinates are first identified from the SZA observa- 
tion (s) using the interfero metric imaging package Difmap. This position serves as the 



^This is a sufficient approximation since the signal to noise in each individual band is too low to 
leave flux a free parameter to be fit by SZA data alone, and even if a were under or overestimated by 
~ 1, the difference at each end of the band would only be biased by ~ 12%. Even for the strongest 
point sources in these fields 4 mJy, as discussed below) this would result in an error smaller than 
any l-cr fluctuation in a particular band. Furthermore, this method of point source modeling was 
tested by adding artificial point sources to both real and simulated observations, and was found to 
have a negligible impact on cluster parameters extracted from these observations. 
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Figure 6.4 Cleaned SZA 30-GHz map of A1835, showing the detection to have a high 
significance of ~ 22-0" at the peak. This image, made with Difmap using only the 
short baseline data (~ 0.35-1.5 kA), is for presentation purposes only; we fit our data 
directly in u,v-spa.ce, and do not use any interferometric mapping to determine cluster 
properties. Note that for a given -space coverage, any unresolved structure in an 
interferometric SZE map of a cluster essentially "looks" like the synthesized beam; 
this determines the effective resultion of the observation. The synthesized beam is 
depicted in gray in the lower left corner. Contours are overlaid at 2-a intervals. 
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Clean I map. Array; 5SS55555 
CI1226+33 at 30,938 GHz 2008 Moy 15 







200 -200 

Right Ascension (arcsec) 

Map center; RA; 12 26 53,001. Dec; +5-3 J2 44,992 (ZOOao; 
Map peak; -0,00354 Jy/beam 

Ct^ntours; 0.0003 Jy/beam k (-10 -S -6 -4 -2 Z ) 
Beam RVHM; 129 x 114 (arcsec) at -35.9' 



Figure 6.5 Cleaned SZA 30-GHz map of CL1226, showing the detection to have a 
significance of ~ lO-o" at the peak. See caption of Fig. 6.4 for further details. 
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200 -200 

Right Ascension (arcsec) 
Map center; RA; ^2 26 58.000. Dec; +33 32 44.399 (2000.0; 
Map peak; -0,00409 ■Jy/beam 

Contours; O.OOO^J Oy/beam x (-12 -10 -6 -6 -4 -2 2 
Contours; + & ) 

Beam RVHM; 42,3 x :39,1 {arcsec) at -19,9'' 




Figure 6.6 Cleaned SZA 90-GHz map of CL1226, showing the detection with this 
higher resolution instrument to have a significance of ~ 12-0" at the peak. This 
observation was included since 30-GHz data alone could not constrain the radial 
profile of this high-redshift, small angular extent cluster (see Figures 1.10 and 1.11). 
See Fig. 6.4 for further details, noting the short baselines of the 90-GHz instrument 
are ~ 1-4.5 kA 
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input location in markov, which is then refined in trial markov runs. Such MCMC 
trials would sufficiently constrain the location of any newly-discovered point source, 
at least for the purposes of constraining point source contamination in an SZE obser- 
vation. However, for all detected sources in the observations presented here, I found 
their trial positions agreed to within ~ 2" of point sources already found by NVSS or 
FIRST. Since the VLA has superior resolution to the SZA, which constrains sources 
to within ~ 7" positional uncertainty, I simply fixed each source's location to the 
catalogued NVSS or FIRST position. I leave the source fiux a free parameter in the 
Markov chains used to constrain cluster properties, so that the cluster SZE fiux and 
any compact sources are simultaneously fit. This marginalizes over uncertainties in 
the SZE and point source fiuxes. 

In the SZA cluster observations presented here, the A1835 field contains two 
detectable compact sources at 30 GHz: a 2.8 ± 0.3 mJy central source, and a 1.1 ± 
0.3 mJy source ~ 1' from the cluster center (using the fixed position and spectral 
index for each source, as discussed above). The central source was detected in both 
the NVSS and the FIRST survey, while the weaker was only detected by FIRST. 
The SZA observation of CL1226 contains one detectable compact source, identified in 
both FIRST and NVSS, with a fiux at 30 GHz of 4.0 ±0.8 mJy, located ~ 6' from the 
cluster center (see also Muchovej et al.|[2007 ). Flux from three compact sources, with 
positions constrained by NVSS and FIRST, was detected at 30 GHz in the A1914 
field. The fiuxes of the sources are 2.2±0.3 mJy, 1.3±0.3 mJy, and 0.6±0.2 mJy. The 
strongest was detected in both the NVSS and the FIRST survey, while the second 
strongest was only detected in the FIRST survey, and the weakest was only detected 



in NVSS. Figures 6.4-6.7 depict the SZA data with these point sources removed. 
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6.3 Independent X-ray Analysis 

For each of these clusters, I compare the results of the joint SZE+X-ray analysis to the 
results of independent X-ray analyses, performed by Ben Maughan and described in 
detail in M07 and M08, which respectively presented the observations of CL1226 and 
A1914. The ACIS-I observation of A1835 used in this work became public more re- 



cently than M08 was published, and is therefore presented here - and in Mroczkowski 



et al. (2008) - for the first time. The observation of A1835 was calibrated and ana- 
lyzed using the same methods and routines described in M08. Several key differences 
exist between both the M07/M08 data reduction and fitting methods, and those used 



in the joint SZE+X-ray analysis here (described in Bonamente et al. (2004, 2006), as 
well as in Appendix |A]) . I discuss a few salient differences here. 

In the X-ray-only analysis, Maughan used blank-sky fields to estimate the back- 
ground for both the imaging and spectral analysis. The X-ray images were binned 
into 1.97" pixels (as are the images used in the joint SZE+X-ray analysis; see Ap- 
pendix |A]). The imaging analysis (primarily used to obtain the gas emissivity profile) 
is performed in the 0.7-2 keV energy band in order to maximize cluster signal to 
noise. Both the Chandra^s efficiency and the (redshifted) cluster emission are highest 



in this range (see Fig. 1.12); since spectral information is not preserved in the X-ray 
images used in the surface brightness fit, binning noisier photons from outside the 
0.7-2 keV energy band can degrade the cluster signal, as the binning is unweighted 
(Appendix [A| discusses this further). 

In Ben Maughan's spectral analysis (used to derive both the global temperature Tx 
and the temperature profiles), spectra extracted from each region of interest were fit 
in the 0.6-9 keV banc^with an absorbed, redshifted APEC ( Smith et al.|2001 ) model. 



This absorption was fixed at the Galactic value. The global Tx, used in the isothermal 



^ A larger range of photon energies can be used for the temperature determination since photon 
energy information is preserved when performing the spectral fit. 
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/3-model fits, is determined within an annulus with radius r G [0.15, 1.0] rsoo- 

An important consideration when using a blank-sky background method is that 
the count rate at soft energies can be significantly different in the blank-sky fields than 
the target field, due to differences between the level of the soft Galactic foreground 
emission in the target field and that in the blank-sky field. Ben Maughan accounted 
for this in the imaging analysis by normalizing the background image to the count 
rate in the target image in regions far from the cluster center. In his spectral analysis, 
this was modeled by an additional thermal component that was fit to a soft residual 
spectrum (the difference between spectra extracted in source free regions of the target 
and background datasets; see Vikhlinin et al. ( 2005[ )). The exception to this was 



the XMM-Newton data used in addition to the Chandra data for CL1226. A local 
background was found by M07 to be more reliable for the CL1226 spectral analysis, 
thus required no correction for the soft Galactic foreground. 

The M07/M08 X-ray analysis methods exploit the full V06 density and tempera- 



ture models (Eqs. 5.12 & 5.13) to fit the emissivity and temperature profiles of each 
cluster. The results of these fits are used to derive the total hydrostatic mass profiles 
of each system. Uncertainties for the independent. X-ray-only analysis are derived 
by Ben Maughan using a Monte Carlo randomization process. These fits involved 
typically ~1000 realizations of the temperature and surface brightness profiles, fit 
to data randomized according to the measured noise. For further details about this 
fitting procedure, see M07 and M08. 

6.4 SZE Cluster Visibility Fits 



As mentioned in both 91.31 and ^4.5, interferometric SZE data are in the form of 



visibilities, V{u,v). Equation 1.7, which is an integral over all space {fjdxdy), 



relates the visibilities (in flux) to the spatial intensity pattern (in units of flux per 
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solid angle on the sky, or flux density). Each visibility is thereby a measure of the 
flux within the Fourier mode probed by a given baseline. 



Using relations provided in, e.g., Carlstrom et al. (2002), we can rescale the visibil- 



ities to a new quantity Y{u, v), which is the integral of line-of-sight Compton y within 
the Fourier modes probed by an interferometer]^ In the image plane, we can convert 
intensity to line-of-sight Compton y using the CMB intensity Jq, the derivative of the 
blackbody function, and the spectrum of the SZE distortion f{x) as a function of 



dimensionless frequency x (Eq. 1.4). The derivative of the blackbody function mul- 



tiplied by f{x) yields the spectral dependence of the intensity shift in primary CMB 
due to the SZE, which is 



X 4 



The shift in primary CMB intensity due to the SZE, AIsze, is 



(6.2) 



AIszE = gix)Ioy, 



(6.3) 



where the primary CMB intensity Iq is 



2{kBTa 



{hey 



(6.4) 



Here ks is Boltzman's constant, T^mb is the primary CMB temperature, h is Planck's 
constant, a nd c is the speed of light. 



Figures 6.8 -6.10 show representative, high-significance fits of the N07 and isother- 



mal /3-model to each cluster's radially-averaged SZE visibility data (i.e. binned ac- 
cording to u,v distance). This quantity is then rescaled by d\/[g{x)Io] to obtain the 



^Note that the Fourier transform of the hne-of-sight Compton y, Y{u,v), scales hke Yint but is 
not to be confused with the actual Yjnt - in the form of either Ijog or yyoi(see Eqs. 5.36 & 5.37) - 



the recovery of which is discussed in S6.7 
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Figure 6.8 Radially-averaged SZE model fits to A1835 in ■u,f-space, from the jointly- 
fit N07+SVM and isothermal /3-model. The upper panels show the real (left) and 
imaginary (right) components of the visibilities, radially-averaged to be a function of 
{u, v) radius, and rescaled to units of intrinsic, line-of-sight integrated Compton y; I 
plot here the frequency-independent quantity Y{u,v) d\ = V{u,v) d\/g{x)lQ, where 
each band is scaled appropriately before binning, and the angular diameter distance is 
computed using the assumed ACDM cosmology. The lower panels show the reduced 
of the fits to the data for the chosen binning. The black points with error bars 
(1-0") are the binned Y{u,v)d\ data, with the point source models first subtracted 
from the cluster visibilities. The blue, solid line is a high likelihood N07 model fit, 
while the red, dashed line is a similarly-chosen fit of the /3-model. For the available 
data points, both SZE models fit equally well (see lower panels, which shows the 
for each model is indistinguishable). However, note that as the w^f -radius approaches 
zero kA - where there are no data to constrain the models - the /3-model predicts a 
much higher integrated Compton y than the N07 model. For cluster data centered 
on the phase center of the observation (with both cluster and primary beam sharing 
this center), the mean imaginary component would equal zero. Since these model 
fits include the primary beam, which is not necessarily centered on the cluster, the 
small but non-zero imaginary component in the upper right panel is expected (note 
the smaller units for the y-axis of the right hand plot) (See, e.g. Reese et aL]|2002[ for 
comparison) . 
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Figure 6.9 Radially-averaged SZE model fits for CL1226 in u,v-space. See Figure [6^8 
for additional caption details. Note that the 30-GHz -space coverage alone does 
not probe a sufficient range of cluster scales (see Figures 1.5 and |l.ll"| ) to determine 
where the SZE signal falls to zero. This results in poor constraints on the cluster's 
radial profile when using 30-GHz data alone, and is a result of this high-redshift source 
being relatively compact on the sky (compared to A1835 and A1914; note that r25oo 
is on the order of 1 arcminute, as shown in Table 6.4, ^6.7). I therefore included 



90-GHz SZA data (green) in the joint SZE-I- X-ray fit, since the 90-GHz instrument 



was designed to complement the -coverage provided at 30-GHz (see Figures 1.10 & 



1.11). Again, the small but non-zero imaginary component is expected, but note the 



smaller units for the plot of the imaginary component of the radially-binned Y{u,v). 



radially-averaged Y{u,v) intrinsic to the clusterjj This rescaling removes both the 
redshift and frequency dependence from the cluster visibility data, after having first 
subtracted all point source models (as modeled simultaneously with cluster in the 
MCMC fitting code) from the visibility data in u,v-s]iace. 



Note that for Figure 6.9, 1 combine both 30 and 90-GHz data to increase both 
the dynamic range and UjV-coverage on the high-redshift cluster CL1226. The u,v- 
coverage provided by the 30-GHz observation alone was insuffient to constrain the 



'^These fits were performed jointly on the X-ray imaging and SZE interferometric data. However, 
the inclusion of X-ray data does not significantly affect the high-likelihood fits of the N07 profile; 
rather, it excludes some of the low Hkelihood N07 pressure fits (see discussion in ^6.7 1. 
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Figure 6.10 Radially-averaged SZE model fits for A1914 in UjV-space. See Figure 6^ 
for details. Note that the slightly larger imaginary component in the upper right 
panel could be due to cluster asymmetry, since A1914 is disturbed and elliptical (see 
Fig. 6.7). Note that the reduced of this fit is no larger than those for the other 



two clusters (Figures 6.8 & 6.9) 



cluster's profile. Figure 6^ shows the 30-GHz data points at low -radii agree to 
within their 1-a error bars, meaning the SZE signal from the cluster was not resolved 
by those baselines, nor does the 30-GHz data for CL1226 accurately determine where 
the signal falls to zero. The inclusion of 90-GHz data provides the necessary u,v- 
coverage and dynamic range to resolve the cluster profile. As detailed in §1.3.2 , the 
90-GHz instrument was designed to complement the t/jf-coverage provided at 30-GHz 



see Figures 1.10 & 1.11) 



In Fig. 6.10[ it can be seen that the average imaginary component of the radially- 
averaged Y{u,v) of A1914 is slightly larger than that seen for the other two clusters 



(Figures 6.8 & 6.9). As mentioned in the caption, I attribute this to the disturbed. 



asymmetric nature of A1914, which is also apparent in the SZE image (see Fig. 6.7). 



158 




Figure 6.11 X-ray surface brightness profile fits to the clusters. The vertical dashed 
line denotes the 100 kpc core cut. The blue, solid line is the surface brightness com- 
puted using a high-likelihood fit of the N07+SVM profiles (analogous to the SZE fits 



plotted in Figures 6.8-6.10), while red, dot-dashed line is the surface brightness fit 



of a /3-model. Both model lines include the X-ray background that was fit simulta- 
neously with the cluster model (i.e. the plotted lines are the superpositions of each 
set of the background and cluster models, which were fit simultaneously to the X-ray 
imaging data). The black squares are the annularly-binned X-ray data, where the 
widths of the bins are denoted by horizontal bars. The vertical error bars are the l-cx 
errors on the binned measurements. Arrows indicate r25oo and rsoo derived from the 



N07+SVM profiles (see ^6.7). 



6.5 X-ray Surface Brightness Fits 



The X-ray data reduction used in the joint SZE+X-ray analysis follows that from 
Bonamente et al. (2004 2006), and is discussed in detail in Appendix |Aj Each image 
was fit, excluding the inner 100 kpc from the X-ray imaging likelihood calculation 



see ^5.4), out to ~ rsoo, using the Markov Chain Monte Carlo technique, described 



in ^5.4 , which jointly fits the X-ray and SZE data|^ Table 6.3 lists the precise fitting 
regions for each X-ray observation. 



The surface brightness (Eq. 1.8) was modeled separately with both the isothermal 



/3-model, using the spectroscopically-determined, global Tx (measured within r G 
[0.15, 1.0] rsoo), and the SVM, with temperature derived using the ideal gas law from 



the N07 pressure profile fit to the SZE data (as discussed in ^5.3.5). Since the SVM 



was developed explicitly to be jointly fit with some form of the generalized NEW 



^Recall that the SZE data are used here to obtain temperature information in the joint N07+SVM 
fitting procedure. 
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pressure profile, requiring both SZE and X-ray data, the X-ray surface brightness 
fit for the SVM contains information determined by the N07 pressure profile, which 
is primarily determined by SZE data. However, since X-ray surface brightness is a 



much stronger function of density than of temperature (Eq. 1.8), the SVM primarily 
determines the shape of the surface brightness fit. 



Fig. 6.11 shows representative, high-likelihood (low x ) fits to the surface bright- 
ness of each cluster, for both the SVM and isothermal /3-model. For plotting purposes, 
the X-ray data are radially- averaged around the cluster centroid, which is determined 
by fitting the two-dimensional X-ray imaging data with the spherically-symmetric 
SVM and isothermal /3-model profiles. 

6.6 Fit Cluster Gas Profiles 

The profiles parameterized in the models, which I apply in fitting the SZE visibilities 
and X-ray surface brightness data, are Pe{r), ne{r), and Te(r). In this section, I 
discuss the recovered model fits of the parameters, and compare them with those fit 



in the detailed, independent X-ray analysis ([ |6.3 ). For the N07+SVM fits presented 
here, metallicity is fixed to the best-fit value determined from X-ray spectroscopic 
data. 



Figures |6.12f|6.14| show the three-dimensional ICM radial profiles fit in the joint 
analysis of the SZA visibility + Chandra imaging data for A1835, CL1226, and A1914 
(from left to right). In each figure, I compare the results, shown with their respective 
68% confidence intervals, of each of the models and fitting procedures tested. The 
resulting constraints on Pe{r), ne{r), and Te(r), using the MCMC routine to determine 
the probability density of each of the fits, are as follows: 



Pressure: Fig. 6.12 shows the pressure profile fit to the SZE data (blue for the 



N07 pressure profile fit jointly with the SVM profile to describe X-ray imaging) 
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Figure 6.12 Pe{r) for each set of models fit to each cluster. The pressure from the 
jointly-fit N07+SVM is plotted in blue with vertical hatching. Pressure constrained 
by the SZE fit of the isothermal /3-model is plotted using red, dot-dashed lines; note 
that the isothermal /3-model's shape is constrained by X-ray imaging data, and the 
only unique parameter to the SZE data in this fit is the central decrement (SZE 



normalization, see ^5.2.1). Pressure derived from the density and temperature fits 



of the V06 profiles in the independent X-ray analysis is shown in black with grey 



shaded regions. See text in ^6.6 for details. The vertical, black dashed line shows 



?"25oo derived from the N07-I-SVM fits, while the magenta dashed line is for rsoo (see 



Tables 6.4 and 6.5). 



red for the isothermal /3-model), and compares it with pressure derived from 
the V06 density and temperature profiles, which is constrained by the indepen- 
dent X-ray analysis (black), which utilizes X-ray spectroscopic information not 
used in the N07+SVM fittingj^ For all three clusters, the N07 pressure pro- 
file fits agree within their 68% confidence intervals with those predicted by the 
independent X-ray analysis. 

The pressure predicted by the isothermal /3-model is in good agreement with 
the other model fits at r < r25oo; but is systematically higher in the cluster 
outskirts (> r25oo)- This is likely to be due to the fact that cluster temperature 



declines with radius (see Fig. 6.14), since the density estimates agree within 



f ^ ^^500 (discussed in next subsection; see Fig. 6.13). Since the isothermal 



^The ideal gas law is assumed when deriving the electron pressure from the V06 density and 
temperature fits from the independent, X-ray-only analysis. 




Figure 6.13 Pgas{r) for each set of models fit to eacli cluster. Colors and line styles are 



the same as in Fig |6.12[ See text in §6.6| for details. Note that the vertical, black and 
magenta dashed lines show r25oo and rsoo, respectively, derived from the N07+SVM 
fits. 

temperature Tx is emission-weighted, the temperature is biased to favor the 
inner regions of a cluster (even when the core is excluded), where the density 
(and therefore the emission) is higher. The isothermal /3-model predicts a large 
SZE signal in Fourier modes that are longer (lower u,v-Ta.du) than the SZA 



interferometer can probe (see Figures 6.8-6.10). In contrast, the extra degree 
of freedom in the N07 model allows it to capture the true shape of the cluster's 
pressure profile, since it fits pressure directly without relying on X-ray imaging 
to constrain the pressure profile's shape. The X-ray imaging data are of much 
higher significance than the SZE visibility data, so they determine the shape of 
the isothermal /3-model. In contrast, the parameters of the N07 model are only 
linked to the X-ray fit through the temperature derived from the N07+SVM 



profiles (see ^5.3.5). 



Density: Fig. 6.13 shows the density fits to the X-ray imaging data from each 
cluster (where the SVM and isothermal /3-model were jointly fit with the SZE 
data). Since the X-ray imaging data provide high significance measurements of 



Sx (see Fig. 6.11), and since they are most sensitive to density (see Eq. 1.8) 
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the X-ray data dominate the fit density. Since all three clusters have excellent 
X-ray data, density is tightly constrained for all three types of model fits. The 
additional flexibility in the SVM and V06 density profiles, which were built to 
fit departures from the /3-model density profile in the cluster outskirts, were not 
necessary for any of these clusters. 

For all three clusters, the fits of the SVM density generally agree to within 
1-2% of the full V06 density profile over the radial range r G [100 kpc,r25oo]- 
Differences in the fit densities arise at larger radii, due to the differing X-ray 
background fitting procedures! 

The strong agreement between the densities derived from the SVM and the 
isothermal /3-model at large radii is due to the fact that they use the same X-ray 
data and X-ray background levels, and none of the clusters had density profiles 
that diverged significantly from a /3-model for the fitting regions considered here 
(r G [100 kpc, rsoo]). Note that the density fit from the isothermal /9-model is 
overlaid in the plots by the SVM density fit (due to the strong agreement). This 
supports the argument the SZE data do not infiuence the shape parameters of 



the isothermal /3- model (rc and j3, Eq. 5.1). 



Temperature: By not relying on X-ray spectroscopic temperature informa- 
tion when fitting the N07+SVM profiles to the SZE-I- X-ray data, Te(r) can be 
derived. This derived temperature profile provides a means by which we can di- 
agnose how well the X-ray+SZE derived and X-ray spectroscopic temperatures 



independently agree. As shown in Fig. 6.14, the N07+SVM derived temperature 



profiles broadly agree with the V06 temperature profile fit to the X-ray spec- 



^This could explain Maughan's higher fit V06 density at r > r25oo (see Fig. |6.13 ), which would 
arise if the X-ray background was underestimated; as one approaches a signal-to-noise ratio of 1, any 
systematic errors in the X-ray background (noise level) are fit by the cluster density profile. Although 
it appears negligible in the logarithmic plot, any discrepancy integrates to produce non-negligible 
contribution to, e.g., the Compton-y computed for this cluster from the X-ray only analysis (see 



Fig. 6.16) 




Figure 6.14 Te(r) for each set of models fit to each cluster. Colors and line styles are 



the same as in Fig 6.12 Note that the isothermal /^-model's constant Te(r) = Tx, 



constrained by the X-ray data, is plotted using red, dot-dashed lines and dark red 



shading. See text in ^6.6 for details. Note that the vertical, black and magenta dashed 



lines show r25oo and rsoo, respectively, derived from the N07+SVM fits. 

troscopic data. For the relaxed cluster A1835, the median derived N07+SVM 
temperature profile is well within 1-a of the fit V06 Te(r) for nearly all radii out- 
side the 100 kpc core. For CL1226, which is likely spherically-symmetric but has 
not had enough time to form a cool core (as argued in M07), the temperature 
agrees within 1-a at r G [r25oo, ''^soo], but does not agree in the core. 

For the disturbed cluster A1914, the N07-I-SVM derived Te(r) agrees with 
Maughan's fit Te(r) over a larger range than the isothermal Tx does, but fails 
to capture the overall slope of Te(r). Since the N07 pressure profile was derived 
using relaxed clusters, it may not be surprising that the derived Te(r) could be 
biased; however, note that Pe{r) does agree for the N07 and Maughan derived 



pressure (see Fig. 6.12), so the bias in the derived Te(r) could arise from the 



assumption of spherical symmetry being applied to the density profile (which 



in the X-ray image is clearly clumpy; see Fig. 6.3). 

The assumption of isothermality shows good agreement (typically within 1-a) 
within r < r25oo with Maughan's fit of the V06 temperature profile. As expected. 
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the isothermal Tx does not agree with the more sophisticated analysis in the 
cluster outskirts (r > r25oo)- This implies that isothermality is a relatively poor 
assumption at large radii, as the temperature universally declines at large radii 
in fits that have more degrees of freedom. 

The radial electron temperature profiles Te(r) derived from the N07+SVM joint- 
fit model ( ^5.3.5 ) reproduce the spectroscopically-determined Te(r) from deep 



Chandra observations of clusters quite well. This approach provides a unique 
and potentially powerful probe of Te(r) for high-redshift clusters, for which 
X-ray spectroscopic temperatures are difficult and often expensive to obtain. 

6.7 Derived Cluster Properties 



In Tables 6.4 and 6.5, I report the global properties of individual clusters derived 



from the N07+SVM model fits to the SZE+X-ray data. I calculate all quantities 



(following the methods outlined in ^5.5 & 5.6) at overdensity radii r25oo and rsoo, and 



compare them to results from both the jointfit isothermal /3-model analysis and to 
the Maughan X-ray-only analysis. 

At both r25oo and rsoo, the measurements of Y\os derived from the joint N07+SVM 
and the X-ray-only analysis are consistent at the 1-a level, for all three clusters. The 
isothermal /3-model analysis, however, overestimates Fios by ~ 20%-40% at r2500) 
and by ~ 30%-115% at rsoo- This is due to the large contribution to Y\os (at every 
projected radius) from the cluster outskirts, where the /3-model significantly overes- 
timates the pressure. 

In contrast, the determinations of l^oi at r25oo are generally consistent among the 
three analyses, except in the case of CL1226, due to the lower fit pressure. The excel- 
lent agreement at this radius between y^oi derived using fits of either the isothermal 
/3-model or the N07 profile arises directly from the constraints on SZE fiux provided 
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Table 6.4 Y|os, ^ob ^gas, M^ot, and /gas for each model, computed within each model's 
estimate of r25oo- 



by Y{u,v) (see ^6.4). At rsoo, however, the median Y^oi values from the isothermal 
/3-model are ~20%-60% higher than either the N07+SVM or M08 resuhs, due to the 
fact that isothermality is a poor description of the cluster outskirts. With more data 
points probing cluster scales ~ T250O) the good agreement between the isothermal 
/3-model's determination of l^oi and that from the N07+SVM fit is expected at r25oo- 
However, with the X-ray data determining the shape parameters of the /3-model (in 
particular, the value of (3 that can describe density is too shallow to describe pressure 
in the cluster outskirts), the over-constrained isothermal /3- model fails to capture 
accurately Fvoi(r5oo)- 



The N07 pressure profile (Eq. 5.17) - fit to SZE data alone or fit jointly to SZE-I-X- 
ray data - has just two free parameters to describe the SZE data: Pe,j and Vp. In 
the joint fit N07+SVM, these parameters are only linked to the X-ray imaging data 



through the derived temperature (see discussion in [:5.3.5). Figure 6.15 shows the de- 



generacy between Pe,i and Vp when fitting the SZA observations of A1835, comparing 
the SZE-only and the joint SZE+X-ray fits of the N07 profile. This degeneracy is 
similar to that between Vc and /3, when fitting the /3-model to SZE data alone (see 
Grego et al.|200l" for example). These two quantities are not individually constrained 



166 



Cluster 
Model fit 



^500 
(") 



>"500 

(Mpc) 



-5 



5a 



(lOi^M©) 



Mtot 
(IO^Mq) 



Abell 1835 



N07+SVM 


3651^2 


1 


42+0 


13 

09 


Maughan (this work) 




1 


421° 


07 
05 


isothermal /3-model 




1 


4llg 


03 
03 



CL J1226+3332.9 

N07+SVM 



Maughan et al. 



(20071 



isothermal /3-model 
Abell 1914 

N07+SVM 



Maughan et al. 



(20081 



isothermal fl-model 



125 
115 
127; 

425 
448 
461 



+ 12 





971° 


09 
07 


?3 





891° 


02 
02 


+ 10 
-10 





991° 


08 
08 



+37 
-27 


1 


241° 


11 
08 


+24 
-21 


1 


291° 


07 
06 


+ 13 


1 


34lg 


04 
03 



20.871^.^^0 
21.37l?;« 
34.53 



58 
1.78 
1.68 



11 6+1-^ 

I3.9li;i 

5+2.7 
-2.3 



18.3'' 



12.69l^g7 
13.47li;77 
29.08lij^ 



1 7 7Q+'3''^2 
J-'-'»_2.45 



-1.61 
-0.99 
-1.09 
-1.02 



+ 1.50 



17.41: 

21.29: 
9.53 



10.59l°:f3 

11.9111:^1 



ii-oit?:58 

10.78ll:°| 
17.07l°:|I 



1 Q cc + 1.29 
13-66^0.90 

13 Q4+0.64 
lo.M'i_o.52 



+0.27 
0.27 



13.29 



8 44+°-'^'' 
^■^^-0.57 

Q or, + 0.68 

"■■J^-0.60 



10.111^:^7 

10.24l0:« 

ii.05ll5:f3 



10.681?:^^ 

10.6811:^4 
^+0.58 



10.30" 



■0.54 



7.17 



+2.22 
1.48 



5 4Q+0-46 
^•^^-0.47 



7.49 

6.38 
7.49 
8.14 



+ 1.93 
1.60 



1.83 
1.13 
1.29 
1.00 
+0.69 
0.59 



0.128 



+0.017 
-0.020 



r, iqo+0.009 
U.100_o 012 
-,+0.005 
-0.005 



0.129"* 



f, 11O+0.021 
U.ii»_o.022 
(, ,^,+0.008 
U.i01_o 008 
1+0.020 
-0.016 



0.111"" 



0.159 



+0.021 



-0.024 

n iqo+0.015 
O.i3»_o.oig 

"+0.007 
-0.007 



0.136"" 



Table 6.5 l^os, K 
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vol, Mgas, Mtot, and /; 



gas 



for each model, computed within each model's 
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Figure 6.15 lios computed within 6' (~ 1.4 Mpc, which is ~ rsoo for this cluster) for 
fits to SZA observation of A1835. The bold, black contours contain 68% and 95% 
of the accepted iterations to the jointly-fit Chandra + SZA data, while the thinner, 
blue contours are those for fits to SZA data alone. Note that the vertical, black and 
magenta dashed lines show r25oo and rsoo, respectively, derived from the N07+SVM 
fits. 



by our SZE observations, but they are tightly correlated. The preferred region in the 
Pe,i — 'fp plane encloses approximately constant Fios- As a result, the 68% confidence 
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region allows less than ±25% variation in Yi^s, despite the much larger variation in- 
dividually in Pe,i and r^. Since Fios at large radii scales with the cluster SZE flux 



probed by an interferometer (e.g. Y{u,v) discussed in ^6.4), this parameter is most 
directly constrained when fitting models to the SZE data. 

The gas mass estimates - computed by integrating the gas density fit with either 
the jointly-fit SVM or the isothermal /3-model (as discussed in ^5.5.3 ) - agree with 



the gas mass estimates derived from the Maughan X-ray fits (Tables pA and 6.5). 
This agreement is not surprising, given that the gas mass is determined in all cases 
from density fits to the X-ray surface brightness data. It demonstrates, however, that 
the 100 kpc core makes a negligible contribution to the total gas mass even at r25oo, 
and that excluding the core from the joint analysis does not therefore introduce any 
significant bias in our estimate of Mgas. Incidentally, it also shows that the additional 
components in the SVM and the full V06 density models were not necessary to fit 
these clusters. 



Tables 6.4 and 6.5 also present estimates of the total masses, computed using 



each model's estimate of the overdensity radius (r^,, Eq. 5.34) for each Monte Carlo 
realization of the fit parameters. For two of the clusters, I find that the error bars 
are significantly larger for Mtot determined from the N07+SVM fits than for the 
isothermal /3-model or M08 fits. This is a consequence of the fact that the /3-model 
analysis, with fewer free parameters and the assumption of isothermality, typically 
places strong but poorly-motivated priors on the total mass, while the M08 fits make 
use of the spatially resolved temperature profile afforded by the deep X-ray observa- 
tions. CL1226 is the exception, with both the N07+SVM and M07 constraints on 
Mtot being tighter than those provided by the isothermal /3-model; the N07+SVM 
fits rely on high-significance SZE constraints provided by combining 30+90 GHz SZA 
data, while the M07 constraints are narrowed by using both XMM-Newton and Chan- 
dra spectroscopy. I find that the N07+SVM and M08 total mass estimates broadly 




Figure 6.16 y{R) - Compton y (integrated along the line of sight, Eq. 1.2 ) as a function 



of sky radius R. Colors and line styles are the same as in Fig |6.12[ See text in §6.6 
for details. Note that the vertical, black and magenta dashed lines show r25oo and 
rsoo, respectively, derived from the N07+SVM fits. 

agree at both r25oo and t^qq, leading to good overall agreement between gas fractions 
computed using the N07+SVM profiles and those from the Maughan X-ray fits. As 



discussed in ^5.5.5, I compute the gas mass fraction /gas in a way that takes advantage 



of the fact that MCMC explores the probability density distribution of model fits to 
the data. Both the Maughan V06 and N07+SVM estimates of /gas within r25oo are all 



consistent with the constant fraction found in Allen et al. (2004, 2007). The isother- 



mal /3-model, on the other hand, is too constrained to agree with the non-isothermal 
fits at both r25oo and rsoo; its estimate of Mtot is moreover sensitive to the annulus 



within which Tx is determined. This trend can also be seen in Figure 6.19 which 
shows Mtotif') for each cluster. 



In addition to the data presented in Tables 6.4 and 6.5, which only show parame- 



ters computed at each model's estimate of the overdensity radii r25oo and rsoo, I plot 
the derived cluster profiles and discuss them here: 



Compton y: Figure 6.16 shows y{R) (integrated along the line of sight through 
the cluster, Eq. 1.2), for sky radius R, from the isothermal /3-model and the 
N07 model, and compares it to y{R) derived from Maughan's fit density and 




Figure 6.17 yios - Compton y integrated over sky radius R. Colors and line styles are 



the same as in Fig |6.12[ See text in §6.6| for details. Note that the vertical, black and 
magenta dashed lines show r25oo and rsoo, respectively, derived from the N07+SVM 
fits. 

temperature profiles. This quantity is simply the line-of-sight integral of electron 



pressure (see Eq. 1.2). It is worth noting that the excess pressure predicted by 
the /3-model can translate to an over-prediction of y{R) even at the cluster core. 
Because this excess is nearly constant over sky radius, it is not constrained by an 
interferometer, which is not sensitive to scales larger than a baseline probes (see 



discussion in ^1.3). y{R) from the NOT fits to A1835 and A1914 agrees well with 
that derived from the independent X-ray analysis, but disagrees to the extent 
that the N07 pressure drops off more rapidly for CL1226 than the X-ray-derived 
pressure. Simply put, the V06 density could be biased by assumptions about 
the X-ray background, leading to an overestimate of the density at large radii, 
which would contribute to the X-ray estimate of y{R) (see discussion of the 



density fits in ^6.6). While this may be disconcerting, the reader should note 
this is one of the drawbacks when using quantities integrated along an infinite 
sight line through the cluster; if the profile does not decline rapidly enough, the 
overestimate at large radii yields a non-negligible contribution. 

Flos ~ the line-of-sight Integrated Compton y, integrated over an area 




Figure 6.18 Mgas - the gas mass integrated within a spherical volume defined by 
cluster radius r. Colors and line styles are the same as in Fig 6.12 See text in §6.6 
for details. 



of the sky: Figure 6.17 shows the SZE scaling quantity l^os integrated along 



the line of sight through the cluster and within a region of the sky (see ^5.6.1 



Eq. 5.36). Given the isothermal /^-model's y{R), and the fact that the integral 



for Yios diverges for /3 < 1, it is unsurprising that there is a large systematic 
difference, at large radii, between the models' estimates. For A1835 and A1914, 
the estimates of llos(-R) from the SVM+N07 fits agree with those calculated 
from Maughan's fits, which is expected given the agreement in y{R). 



M, 



gas 



Gas Mass: Figure 6.18 shows the gas mass obtained by integrating 



the gas density within a spherical volume (using Eq. 5.29). Since the densities 



obtained using each method of fitting were similar (see Fig. 6.13), Mgas(r) ob- 



tained from each method is also similar. Note, however, that the independent 
X-ray analysis also modeled the cluster core, which contributes negligibly to the 
overall gas mass at large radii (such as r25oo and rsoo)- 

Mtot~ Total Mass: Figure |6.19| shows the total mass profile estimated by 



applying hydrostatic equilibrium to each model (Eq. |5.32). The total mass 



from the N07+SVM agrees with that from the independent X-ray analysis for 
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Figure 6.19 Mtot - the total mass estimated assuming hydrostatic equihbrium at 
cluster radius r. Colors and line styles are the same as in Fig 6.12 See text in §6.6 
for details. 
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Figure 6.20 1-D histograms of the N07+SVM jointfit estimates, for A1835, of Mtot 
and Flos normalized by their respective median values, Mtotis the cyan region with 
a dashed outline, and Yjos is the vertically hatched region with a solid black outline. 
Both Mtot and Fios are computed within a fixed radius of 6' = 360". The derived Fios, 
which scales with integrated SZE flux, has a more tightly constrained and centrally 
peaked distribution than that of Mtot, as Mtot is sensitive to the change in slope 
of the pressure profile (see Eqs. 5.16 & 5.32). Since the N07 model was developed 
primarily to recover Yjnt from SZE observations, it is useful to note how well this 
model performs in this capacity. 
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most cluster radii, for all three clusters. In contrast, the isothermal /3-model 
disagrees for A1835 and A1914 over a large range. Due to the steepness of the 



derived N07+SVM temperature profile of A1914 (see Fig. 6.14), its estimate of 
Mtot flattens before rsoo- 



Since the density fits are all well-constrained (see Fig. 6.13), the error bars 
from the isothermal /3-model and the V06 fits performed by Ben Maughan arise 
almost entirely from the temperature constraints (including constraints on the 
derivative of temperature), while the error bars on Mtot from the N07+SVM 



fits arise from those on the pressure profile and its derivative. Fig. |6.20| shows a 
comparison of the ability of the N07+SVM, when fit to X-ray imaging and SZE 
data, to constrain Mtot and l^os- The sparsely sampled u,v-data provide poorer 
constraints on dP/dr than they do on Yios, which scales with SZE flux. The 
relatively poor constraints on dP/dr translate directly to error bars on Mtot- 

Unlike the N07-I-SVM, the isothermal /5-model lacks the flexibility (by defini- 
tion) to model a gas temperature profile that varies over the cluster's radius. In 



Fig. |6.19 one can see that the total mass estimate from the isothermal /3-mode 



never agrees with the independent, detailed X-ray-only analysis at both r25oo 
and rsoo- Where the isothermal /3-model does agree depends mainly on the re- 
gion within which Tx was measured. Since the HSE estimate of Mtot depends 
on the slope of the pressure profile (as discussed in ^:5.32)j^the fact that the 



isothermal /3-model assumes a constant temperature {dT/dr = 0), but recovers 



density well (e.g. Fig. 6.13), indicates the isothermal /3-model will provide a 



systematically-biased estimate of Mtot at some radii. 



/gas ~ Hot Gas Fraction: Figure 6.21 shows the gas fraction, computed using 



each accepted MCMC iteration's prediction for Mgas and Mtot (as discussed in 

Assuming the ideal ga 
nkB{dT/dr) + kBT{dn/dr). 



^ Assuming the ideal gas law, as we do when performing the HSE estimate of A/tot, dP/dr 




Figure 6.21 /gas - the gas mass fraction, computed as Mgas/Mtot at cluster radius r 
for each accepted iteration in the Markov chain. Colors and line styles are the same 



as in Fig 6.12 See text in ^6.6 for details 



^5.5.5). The N07+SVM estimates of /gas at intermediate radii (r ~ ^"2500) are 



consistent with previous, X-ray-only results (such as Allen et al. 2004, 2007). 



The flatness of the Mtot profile derived from N07+SVM fits to A1914 leads to 

assuming 



the prediction that /gas.soo > ^h/^M ~ 0.165 (Hinshaw et al. 2008 



h = 0.7), which seems unlikely. This suggests that the N07+SVM estimate of 
Mtot for A1914 is too low at r ~ rsoo, due to the steeply declining temperature 
profile discussed above. 



6.8 Conclusions 

I have applied a new model for the ICM pressure profile - motivated both by theory 
and detailed cluster observations - to fit the SZE signal from three galaxy clusters. I 
have also developed and tested a complementary density model - a simple extension 
to the /3-model - that can be used jointly with this new SZE model to fit X-ray 
observations. I show the new pressure profile accurately captures the bulk properties 
of relaxed clusters outside the core, and out to rsoo- 

I also argue that this new model should supplant the isothermal /3-model when 
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attempting to use SZE data - interferometric or otherwise - to determine cluster SZE 
scaling relations. The additional degree of freedom, versus the isothermal /3-model, 
in the N07 pressure profile allows it to more accurately describe a cluster that has 
been only sparsely sampled in u,v-space. 

Finally, the derived temperature from the N07+SVM could prove to be a useful 
tool, either in independently confirming X-ray spectroscopic temperature measure- 
ments, or in measuring the temperatures of high-redshift clusters for which suffi- 
ciently deep X-ray exposures are unavailable, and are difficult to obtain. I present in 
Chapter [7] a way in which X-ray spectroscopic data can be used to provide additional 
constrains in the context of the N07+SVM profile. 



Chapter 7 
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Extensions to the models 

Throughout the cluster analysis presented in Chapter |6| I tested the new pressure 
and density models by discarding all X-ray spectroscopic temperature information 
from the joint SZE+X-ray fit, and using this information as an independent test. 
I used the SZE constraints on pressure to derive temperature, assuming the ideal 
gas law and a fixed ACDM cosmology with Qu = 0.3, Qa = 0.7, fik = 0. In this 
section, I consider what additional tests can be performed by including spectroscopic 
temperature information. 

For these proposed tests, we will choose a sample of relaxed clusters comprising 
many of the 6 > —10° clusters with publicly-available Chandra observations. By 
selecting apparently relaxed clusters, the impact of assuming spherical symmetry can 
be reduced. 



7.1 Using X-ray Spectroscopic Data 



Medium exposure observations of high redshift clusters typically provide too few pho- 
tons to constrain their temperature profiles in detail. The most robustly-determined 
spectroscopic temperature for any cluster is a single, global Tx- Measure Tx within a 
core cut annulus. This choice attempts simply to probe the more self-similar portions 
of a cluster, avoiding the systematic discrepancies between "cool-core" and "non-cool- 



core" clusters (see, for example, Kravtsov et al. (2006); Maughan (2007), where a core 



cut was used to make Yx and Lx more robust proxies for Mtot)- 

The observable, emission-weighted spectroscopic temperature Tx can be predicted 
if a cluster's density and temperature distributions are somehow known (e.g. for a sim- 
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ulated cluster, or when using the N07+SVM profile without including spectroscopic 



information). Mazzotta et al. (2004) did precisely this, and verified this method is 



reliable for temperatures T > 3.5 keV. Vikhlinin (2006) provides a simple, integral 



form for the "spectroscopic-like" temperature Tgi, which for spherically-symmetric 
profiles reduces to 



{7.1] 



where the temperature weighting factor a = 0.75. 

Using the ideal gas law, and assuming the electron temperature is equal to the 
overall ICM temperature (Te(r) = T(r)), we can derive Tsi from fits of the N07+SVM 



profiles to X-ray and SZE imaging alone. Eq. 7.1 becomes 



r«i = k 



-1 Jr^i 



B 



(7.2) 



By ensuring Tgi is computed over the same volume within which Tx was measured, 
one can include the likelihood that Tx = T^i in the MCMC fitting process (described 



in ^5.4). This allows for the tests proposed in the following sections, which rely on 



the combination of SZE data with X-ray imaging and spectroscopic dataj^ 



7.2 Refining Constraints on M^^^ and / 



gas 



While the ICM gas density is well-constrained by X-ray imaging, large uncertainties in 
cluster pressure and temperature remain when using either X-ray spectroscopically- 



measured or X-ray -|-SZE-derived temperature alone (see ^6.6). These uncertainties 
dominate the uncertainty in M^ot estimates (assuming systematics in the HSE mass 
determination can be accounted for). However, one can use the combined data to 



^ Note that using Tx ~ Tgi does not force the temperature profile T{r) to be isothermal; it is sim- 
ply stating that, with a known T(r) and n{r), one can compute the emission- weighted temperature 
Tx that would be measured within a given region of an X-ray observation. 
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Figure 7.1 Constraints on Mtot{r) and fgas{r) using the N07+SVM combined with 
the spectroscopically-measured temperature, using Tx = Tgi, to fit A1835. The 
spectroscopically-measured Tx used here is the same as that used in the isothermal 
/3-model analysis presented in Chapter [6] The N07+SVM results that include X-ray 
spectroscopy are plotted using red, dot-dashed lines and red shading. The N07+SVM 
results without spectroscopic constraints (i.e. the resulted detailed in Chapter [6]) are 
plotted in blue with vertical blue hatching. Results derived from the density and 
temperature fits of the V06 profiles, from the independent X-ray analysis, are shown 
in black with grey, shaded regions. Both panels show that including spectroscopic 
information in fits of N07-I-SVM profile tightens constraints and improves the already 
remarkable agreement between it and the independent X-ray analysis. 



constrain more tightly the observable properties of a cluster, using the direct SZE 
pressure measurements to complement the X-ray-derived pressure. The combination 
of high-significance X-ray and SZE data on relaxed, spherically-symmetric clusters 
will allow for detailed astrophysical measurements of these systems. 

For these studies of cluster astrophysics, we continue to assume a fixed (Ia for each 



cluster. A test of how well this performs is presented in Figure 7.1 



7.3 Sensitivity to Angular Diameter Distance ((Ia) 

Rather than constraining Mtot more precisely, we can relax the assumption of cos- 
mology by including spectroscopic information in the joint fit. We consider here the 
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Cluster 


z dx" 
(Mpc) 


Tx'" dA^ Tsi 
(keV) (Mpc) (keV) 


A1835 
CL1226 


0.25 806.5 
0.89 1601.8 


^•^^-0.37 ""'-152 ^-^^-0.54 

9.30li:i ^^Ortttt 9AAt\fs 



Table 7.1 Clusters chosen for testing the ability of the upgraded N07+SVM profile to 
constrain cosmology. 



"Computed angular diameter distance, assuming Qm ~ 0.3, fl^ — 0.7, & 17^ = 0. 

''Global X-ray spectroscopic temperatures were determined in the range r G [0.15, 1.0] r^oQ. 

■^Angular diameter distance fit to the SZE + X-ray imaging and spectroscopic data. 

cosmological sensitivity of our fits to the data. 



The Compton y parameter (Eq. 1.3), which scales as electron pressure integrated 
along sight line i, has an inverse linear dependence on dA, since di = dAdO: 



P,di 



Pe OC dj^. 



(7.3) 



Because X-ray imaging is sensitive to the surface brightness Sx (Eq- W^)-, X-ray 



'1/2. 



constrained density is dependent on 

Sx = ^^^^^^Y j ^eAee(Te, Z) di = 

Combining these dependencies, we see 

Tsl OC Pe/Ue OC d^^^"^ . 



He OC Ua 



d,~'l\ 



(7.4) 



(7.5) 



The spectroscopically-measured temperature Tx does not have any dependence on 
dA, so setting Tx = Ts\ results in a sensitivity to cosmology. This is the same 

cosmological sensitivity exploited in previous, joint X-ray-|-SZE dA determinations 



such as Bonamente et al. (2004), but relies on a more sophicated pair of models for 



the cluster gas. 
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The inverse square root cosmological sensitivity of a joint SZE+X-ray fit to a clus- 
ter implies that, for example, a single cluster observation with ±10% uncertainty in 
the measurements of Tgi and Tx (the measurements of which should be independent) 
can constrain to within ±37.6%. Here, I have added the errors on Tx and Tgi in 
quadrature, since constraints provided by X-ray spectroscopic and the temperature 
derived from X-ray±SZE imaging should be independent. A determination of the 



angular diameter distance benefits from a large sample of clusters (e.g. Molnar et al. 



2002 Bonamente et al. 2006), and we can expect a sample of ~ 40 relaxed clusters 



to yield ~ 5% statistical uncertainty on dA (where I have ignored systematics due 
to cluster asphericity and systematic effects in X-ray temperature measurements; see 



for example Ameglio et al. (2006)). 

I tested this method on the relaxed clusters A1835 and CL1226, using the global 
Tx, measured within r e [0.15, 1.0] rsoo, from Chapter [6j The results of this test are 
consistent with ACDMand are shown in Table 17.11 

It has been suggested (Alexey Vikhlinin, private communication) that a compari- 
son between the X-ray proxy Yx = M^^^^Tx and the SZE quantity l^oi could be used 
to constrain dA more robustly than previous X-ray±SZE attempts, since both Yx and 
Fvoi are integrated quantities that can be robustly determined. However, since the 



N07±SVM profiles share a density profile, simple inspection of Eq. 7.1 reveals that 
setting Tx = Ts\ is in fact the correct implementation of comparison between Yx and 
Yvoi- The X-ray spectroscopic Tx and the derived Tgi contain the same weighting by 
density, while Yx and l^oiare not in general equivalent. Since l^oi is the integral of 
pressure, it scales as Mg^sT^^, where Tmw is the mass-weighted temperature. In gen- 
eral, the emission-weighted temperature Tx is not equivalent to the mass-weighted 
temperature, as Tx favors denser regions more heavily than T^^ does, and clusters 
are not isothermal. 



Recently, Allen et al. (2007) used X-ray observations of clusters to constrain dA 
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by assuming the gas fraction /gas('"250o) is constant with redshift. They find /gas = 
Mgas/Mtot oc (^1^1 which is far more sensitive to dA than the hne-of-sight SZE+X-ray 
joint fit when no assumptions are made about the evolution of /gas- 

For SZE+X-ray joint fits without spectroscopy, Mgas has same dependence as it 
does when constrained by X-ray data alone, which is 



5/2 



(7.6) 



Volume is proportional to the cube of angular diameter distance {V oc 6/^4), and 
density is proportional to the inverse square-root of dj^ ine oc rf^^^^^, as implied 



by Eq. 7.4). For the Mtot cosmological dependence, we have a slightly different 



cosmological dependence than that from X-ray alone: 



tot I 



Gp. 



gasV 



dP, 



gas 



dr 



oc 



d 



^d- 
1/2 "A 



73/2 



(7.7) 



using r = 9 dA and Eqs. |7.3| & |7.4[ Since the SZE-determined pressure exhibits a 
different dependence on cosmology than the X-ray-determined pressure, we arrive at 



/gas cx: dA- 



(7i 



We conclude that an SZE+X-ray determination of dA using an assumed evolution 



of /gas would provide precisely the same sensitivity as Allen et al. (|2007|) exploit, since 

1/2 

we recover an additional dy^ dependence through Eq. 



7.5 



However, the additional. 



independent constraints afforded by the inclusion of SZE data could provide a more 
precise determination than possible with X-ray data alone. 
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7.4 Final Words 

In this thesis, I have presented an exciting new tool for observing the Sunyaev- 
Zel'dovich effect from galaxy clusters. I have shown how data from this - or any 
interferometric SZE instrument - can be combined with X-ray imaging to place rea- 
sonable constraints on galaxy clusters, including at high redshift where the SZE can 
contribute most significantly to the understanding of large scale structure]^ Finally, I 
have shown how the new models can be combined with X-ray spectroscopy to improve 
constraints and perform key measurements of the expansion of the Universe. 



^These methods are, of course, more broadly applicable, and could be applied to non- 
interferometric SZE instruments. 
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Appendix A 

Chandra X-ray Data Analysis 

The X-ray data used in this analysis were obtained with the Chandra ACIS-I 
detector, which provides spatially resolved X-ray spectroscopy and imaging with an 
angular resolution of 0.492" and with energy resolution of ~ 100-200 eV. Data analysis 
was performed with the CIAOj^ software (version 3.2) and the CALDB calibration 
information (version 3.1) provided by the Chandra calibration team. 

Both images and spectra of the low-redshift clusters were limited to the 0.7- 
7.0 keV energy band in order to exclude the low-energy and high-energy data that 
are more strongly affected by background and by calibration uncertainties. For the 
high-redshift cluster, CL1226 {z = 0.89), the image was limited to the 0.7-2.0 keV 



band, where Chandra^s efficiency peaks (see Fig. 1.12). This range of energies was 
chosen because, with a spectroscopic temperature of ~9.8 keV, only ~ 50% of the 
redshifted cluster emission is measured at photon energies > 2 keV; However, the X- 
ray background in a 0.7-7.0 keV image is ~ 4 times higher than that in a 0.7-2.0 keV 
image. Combining this with the fact that the ACIS-I efficiency peaks below 2 keV 



(see Fig. 1.12), the highest S/N is obtained in the detector energy range 0.7-2.0 keV. 

The X-ray images were binned in 1.968" pixels; this sets the limiting angular 
resolution of our X-ray data, as the Chandra point spread function in the center of 
the X-ray image is smaller than our adopted pixel size. The X-ray background was 
measured for each cluster exposure, using peripheral regions of the adjacent detector 
ACIS-I chips that are source free. Additional details of the Chandra X-ray data 



analysis are presented in Bonamente et al. (2004, 2006). 



^ Chandra Interactive Analysis of Observations, 



http : //cxc . harvard . edu/ ciao/ 



